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.

Fluid-Solid Interaction in Particle-Laden Flows

Liu Qingquan, Singh V.P

Prediction of particle distribution in isotropic turbulence by large-eddy simulation

The recent application of large-eddy simulation (LES) to particle-laden turbulence requires that the LES with a subgrid scale (SGS) model could accurately predict particle distributions. Usually, a SGS particle model is used to recover the small-scale structures of velocity fields. In this study, we propose a rescaling technique to recover the effects of small-scale motions on the preferential concentration of inertial particles. The technique is used to simulate particle distribution in isotropic turbulence by LES and produce consistent results with direct numerical simulation (DNS). Key words: particle distribution, particle-laden turbulence, large-eddy simulation, subgrid scale model.

Agglomeration front dynamics: Drying in sessile nano-particle laden droplets

The drying of sessile, nano-silica laden water droplet is studied under ambient conditions, in the absence of any convection. The drying process can be divided into two distinct regimes. During regime 1, the outer edge of the droplet remains pinned and particles agglomerate at the droplet periphery similar to the traditional coffee ring. However in regime 2, with further evaporation, both the liquid contact line and the agglomeration front starts moving radially inwards from the initial contact edge. The contact between the liquid and the agglomerate is maintained throughout regime 2 and the vaporisation driven liquid edge recession essentially drives the inward growth of the particle deposition. Fast kinetics of particle aggregation results in rapid growth of this agglomeration front as seen from the experiments. A theoretical formulation involving a simplistic model of the agglomeration front growth based on particle mass balance has been proposed. (C) 2014 Elsevier Ltd. All rights reserved,

Subgrid-scale contributions to Lagrangian time correlations in isotropic turbulence

The application of large-eddy simulation (LES) to particle-laden turbulence raises such a fundamental question as whether the LES with a subgrid scale (SGS) model can correctly predict Lagrangian time correlations (LTCs). Most of the currently existing SGS models are constructed based on the energy budget equations. Therefore, they are able to correctly predict energy spectra, but they may not ensure the correct prediction on the LTCs. Previous researches investigated the effect of the SGS modeling on the Eulerian time correlations. This paper is devoted to study the LTCs in LES. A direct numerical simulation (DNS) and the LES with a spectral eddy viscosity model are performed for isotropic turbulence and the LTCs are calculated using the passive vector method. Both a priori and a posteriori tests are carried out. It is observed that the subgrid-scale contributions to the LTCs cannot be simply ignored and the LES overpredicts the LTCs than the DNS. It is concluded from the straining hypothesis that an accurate prediction of enstrophy spectra is most critical to the prediction of the LTCs.

Large-eddy simulation of turbulent collision of heavy particles in isotropic turbulence

The small-scale motions relevant to the collision of heavy particles represent a general challenge to the conventional large-eddy simulation (LES) of turbulent particle-laden flows. As a first step toward addressing this challenge, we examine the capability of the LES method with an eddy viscosity subgrid scale (SGS) model to predict the collision-related statistics such as the particle radial distribution function at contact, the radial relative velocity at contact, and the collision rate for a wide range of particle Stokes numbers. Data from direct numerical simulation (DNS) are used as a benchmark to evaluate the LES using both a priori and a posteriori tests. It is shown that, without the SGS motions, LES cannot accurately predict the particle-pair statistics for heavy particles with small and intermediate Stokes numbers, and a large relative error in collision rate up to 60% may arise when the particle Stokes number is near St_K=0.5. The errors from the filtering operation and the SGS model are evaluated separately using the filtered-DNS (FDNS) and LES flow fields. The errors increase with the filter width and have nonmonotonic variations with the particle Stokes numbers. It is concluded that the error due to filtering dominates the overall error in LES for most particle Stokes numbers. It is found that the overall collision rate can be reasonably predicted by both FDNS and LES for St_K>3. Our analysis suggests that, for St_K<3, a particle SGS model must include the effects of SGS motions on the turbulent collision of heavy particles. The spectral analysis of the concentration fields of the particles with different Stokes numbers further demonstrates the important effects of the small-scale motions on the preferential concentration of the particles with small Stokes numbers.

Flow field measurements of pulverized coal combustion using optical diagnostic techniques

This paper demonstrates the application of laser Doppler velocimetry (LDV) and particle image velocimetry (PIV) techniques to a particle-laden reacting flow of pulverized coal. A laboratory-scale open-type annular burner is utilized to generate velocity profiles of coal particles and micrometric alumina particles. Pair-wise two-component LDV measurements and high-speed stereo PIV measurements provide three-dimensional velocity components of the flow field. A detailed comparison of velocities for alumina and coal particle seeding revealed differences attributed to the wide size distribution of coal particles. In addition, the non-spherical shape and high flame luminosity associated with coal particle combustion introduces noise to the Mie scatter images. The comparison of mean and RMS velocities measured by LDV and PIV techniques showed that PIV measurements are affected by the wide size distribution of coal particles, whereas LDV measurements become biased toward the velocity of small particles, as signals from large particles are rejected. This small-particle bias is also reflected in the spectral characteristics for both techniques, which are in good agreement within the range of frequencies accessible. PIV measurements showed an expected lack of response of large coal particles to the turbulence fluctuations. The overall good agreement between LDV and PIV measurements demonstrates the applicability of the high-speed PIV technique to a particle-laden, high luminosity coal flame while highlighting some of its limitations. © 2013 Springer-Verlag Berlin Heidelberg.

Revisiting the carrier-phase source-term formulation in mixed Lagrangian-Eulerian particle models

The formulation of the carrier-phase momentum and enthalpy source terms in mixed Lagrangian-Eulerian models of particle-laden flows is frequently reported inaccurately. Under certain circumstances, this can lead to erroneous implementations, which violate physical laws. A particle- rather than carrier-based approach is suggested for a consistent treatment of these terms.

The transport of cosmic rays in self-excited magnetic turbulence

The process of diffusive shock acceleration relies on the efficacy with which hydromagnetic waves can scatter charged particles in the precursor of a shock. The growth of self-generated waves is driven by both resonant and non-resonant processes. We perform high-resolution magnetohydrodynamic simulations of the non-resonant cosmic ray driven instability, in which the unstable waves are excited beyond the linear regime. In a snapshot of the resultant field, particle transport simulations are carried out. The use of a static snapshot of the field is reasonable given that the Larmor period for particles is typically very short relative to the instability growth time. The diffusion rate is found to be close to, or below, the Bohm limit for a range of energies. This provides the first explicit demonstration that self-excited turbulence reduces the diffusion coefficient and has important implications for cosmic-ray transport and acceleration in supernova remnants.

Insight into instabilities in burning droplets

The complex multiscale physics of nano-particle laden functional droplets in a reacting environment is of fundamental and applied significance for a wide variety of applications ranging from thermal sprays to pharmaceutics to modern day combustors using new brands of bio-fuels. Formation of homogenous nucleated bubbles at the superheat limit inside vaporizing droplets (with or without nanoparticles) represents an unstable system. Here we show that self-induced boiling in burning functional pendant droplets can produce severe volumetric shape oscillations. Internal pressure build-up due to ebullition activity ejects bubbles from the droplet domain causing undulations on the droplet surface and oscillations in bulk. Through experiments, we establish that the degree of droplet deformation depends on the frequency and intensity of these bubble expulsion events. In a distinct regime of single isolated bubble residing in the droplet, however, pre-ejection transient time is identified by Darrieus-Landau evaporative instability, where bubble-droplet system behaves as a synchronized driver-driven system with bulk bubble-shape oscillations being imposed on the droplet. The agglomeration of nanophase additives modulates the flow structures within the droplet and also influences the bubble inception and growth leading to different levels of instabilities. (C) 2014 AIP Publishing LLC.

Classical plume theory: 1937-2010 and beyond

Developing a theoretical description of turbulent plumes, the likes of which may be seen rising above industrial chimneys, is a daunting thought. Plumes are ubiquitous on a wide range of scales in both the natural and the man-made environments. Examples that immediately come to mind are the vapour plumes above industrial smoke stacks or the ash plumes forming particle-laden clouds above an erupting volcano. However, plumes also occur where they are less visually apparent, such as the rising stream of warmair above a domestic radiator, of oil from a subsea blowout or, at a larger scale, of air above the so-called urban heat island. In many instances, not only the plume itself is of interest but also its influence on the environment as a whole through the process of entrainment. Zeldovich (1937, The asymptotic laws of freely-ascending convective flows. Zh. Eksp. Teor. Fiz., 7, 1463-1465 (in Russian)), Batchelor (1954, Heat convection and buoyancy effects in fluids. Q. J. R. Meteor. Soc., 80, 339-358) and Morton et al. (1956, Turbulent gravitational convection from maintained and instantaneous sources. Proc. R. Soc. Lond. A, 234, 1-23) laid the foundations for classical plume theory, a theoretical description that is elegant in its simplicity and yet encapsulates the complex turbulent engulfment of ambient fluid into the plume. Testament to the insight and approach developed in these early models of plumes is that the essential theory remains unchanged and is widely applied today. We describe the foundations of plume theory and link the theoretical developments with the measurements made in experiments necessary to close these models before discussing some recent developments in plume theory, including an approach which generalizes results obtained separately for the Boussinesq and the non-Boussinesq plume cases. The theory presented - despite its simplicity - has been very successful at describing and explaining the behaviour of plumes across the wide range of scales they are observed. We present solutions to the coupled set of ordinary differential equations (the plume conservation equations) that Morton et al. (1956) derived from the Navier-Stokes equations which govern fluid motion. In order to describe and contrast the bulk behaviour of rising plumes from general area sources, we present closed-form solutions to the plume conservation equations that were achieved by solving for the variation with height of Morton's non-dimensional flux parameter Γ - this single flux parameter gives a unique representation of the behaviour of steady plumes and enables a characterization of the different types of plume. We discuss advantages of solutions in this form before describing extensions to plume theory and suggesting directions for new research. © 2010 The Author. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.

Polycyclic aromatic hydrocarbon emission loads and profiles from the burning of locally grown wood species in Victoria, Australia - Some pilot project results

A testing facility for combustion of biomass and sampling of emissions has been established at Deakin University. In this pilot project using this facility, four kinds of locally grown wood species were burned and the particle emissions sampled and analysed for Polycyclic Aromatic Hydrocarbons (PAHs). The selected wood species covering pine, red gum, yellow box and sugar gum, are the most popular domestic fuel wood in Australia. Particulate matter emissions from burning of each load of wood were sampled from the flue using a standard stack emission sampling train. The particle-laden filters were extracted and the .extract analysed to determine PAH concentrations by Gas Chromatographyl Mass Spectrometer (Gc/MS). The sampling was conducted under two different burning conditions with the air inlet of the combustion chamber fully open and with it half open. A suite of 15 PAHs, ranging from naphthalene (C IOHB) to dibenzolahlanthracene (C12H14), were selected for analysis. PAH profiles for the four wood species, under the different burning regimes, have been generated. Some preliminary emission factors for the different wood species have been derived as microgram of summed PAHs (rPAHs) emittedlkilogram of wood burned. Total Particulate Matter (TPM) emission factors were also obtained from gravimetric measurement of the sample filter before and after the combustion. Based on these emission factors, pine displayed the highest level of rPAHs emitted from the combustion of the four wood species, with sugar gum showing the lowest level of rPAHs emission. Emission factors associated with the slow burning condition clearly showed higher l:PAH levels compared to the faster burning condition. During the faster burning condition, red gum and pine show a higher percentage of rPAH to TPM than sugar gum or yellow box. Under the slower burning. the l:PAHlTPM ratio in every case was greater.

Study of polycyclic aromatic hydrocarbons emissions from burning of local wood species in Australia

Four kinds of Australian local wood species were burned in a domestic wood heater installed in a laboratory. The selected wood species include pine, red gum, yellow box and sugar gum, that are the most popular domestic fuel wood in Australia. Particulate matter emissions from burning of each load of wood were sampled isokinetically on filter media from the flue by standard stack emission sampling train. The particle laden filters then went for Gas Chromatography/ Mass Spectrometer (GC/MS) analysis to determine polycyclic aromatic hydrocarbons (PAHs) concentrations. The sampling was conducted under two different burning conditions – air inlet of the combustion chamber fully open and half open. Approximately 15 types of PAHs were detected. Emission factors were derived as microgram of PAHs /kg of wood burned. Total particulate emission factors were also obtained from gravimetric measurement before and after the sampling. PAH emission profiles for four species were generated from the results. Comparisons of emission factors have been conducted among different species of wood, as well as under different burning conditions, ie. fast burning and slowing burning. According to the derived emission factors, pine displayed the highest level of PAHs among the four species, followed by red gum and yellow box, whereas sugar gum showed the lowest level of PAHs. Emission factors were compared between each type of wood under two different burning conditions, the slow burning condition, which was air inlet half open, clearly showed higher PAH levels compared to the fast burning condition. Total PAH fractions on particulate matter were calculated and compared among wood types under two burning conditions. During the fast burning condition, red gum and pine have the higher percentage of PAH to total particulate matter emission than sugar gum and yellow box. When changed to slow burning, the PAH fraction on particulate matter are all increased with sugar gum having the largest increase.

Thinning of a vertical free draining aqueous film incorporating colloidal particles

The drainage under gravity of a vertical foam film formed on a wire frame has been investigated. Dual-wavelength optical interferometry was used so that unambiguous fringe order assignments could be made, enabling absolute film thicknesses to be calculated with confidence. Films were stabilized by nonionic polypropylene glycol surfactant. Halfmicrometer silica particles with varying degrees of hydrophobicity were added to the film-forming liquid to investigate their effect on film drainage rate and stability. Hydrophilic particles had little or no effect, while hydrophobic particles slowed the drainage of the film and caused a minor increase in film lifetime, from ∼10 to ∼30 s. In both the hydrophilic and hydrophobic cases the films ruptured when they reached a thickness of ∼2 particle diameters. Particles of intermediate hydrophobicity had the most significant effect, increasing film lifetime by an order of magnitude over that for hydrophilic particles. The intermediate particles allowed films to thin down to a thickness less than the particle diameter, indicating that particles bridge across the entire film. This did not occur with more hydrophobic particles even though they were embedded in each of the two film surfaces. These results correlate well with previous literature on particle-laden foams. The film thickness and drainage measurements allow drainage mechanisms for the different particles to be identified, thus providing a mechanistic explanation for the observation by several previous authors that foams formed in the presence of particles, for example during mineral processing, have the greatest stability when the particles are of intermediate hydrophobicity.

The simultaneous role of an alveolus as flow mixer and flow feeder for the deposition of inhaled submicron particles

In an effort to understand the fate of inhaled submicron particles in the small sacs, or alveoli, comprising the gas-exchange region of the lung, we calculated the flow in three-dimensional (3D) rhythmically expanding models of alveolated ducts. Since convection toward the alveolar walls is a precursor to particle deposition, it was the goal of this paper to investigate the streamline maps' dependence upon alveoli location along the acinar tree. On the alveolar midplane, the recirculating flow pattern exhibited closed streamlines with a stagnation saddle point. Off the midplane we found no closed streamlines but nested, funnel-like, spiral, structures (reminiscent of Russian nesting dolls) that were directed towards the expanding walls in inspiration, and away from the contracting walls in expiration. These nested, funnel-like, structures were surrounded by air that flowed into the cavity from the central channel over inspiration and flowed from the cavity to the central channel over expiration. We also found that fluid particle tracks exhibited similar nested funnel-like spiral structures. We conclude that these unique alveolar flow structures may be of importance in enhancing deposition. In addition, due to inertia, the nested, funnel-like, structures change shape and position slightly during a breathing cycle, resulting in flow mixing. Also, each inspiration feeds a fresh supply of particle-laden air from the central channel to the region surrounding the mixing region. Thus, this combination of flow mixer and flow feeder makes each individual alveolus an effective mixing unit, which is likely to play an important role in determining the overall efficiency of convective mixing in the acinus.