Stratification by solid particle suspensions

Turbulence was generated by an oscillating grid above a bed of sediment of spherical glass beads. As expected, part of the sediment was lifted up by the grid action and a suspension layer of depth D formed above the grid. This depth was found remaining independent of grid action but varying with the sediment layer depth when the grid action was kept constant. Volume concentration measurements show the existence of only weak concentration gradients over the layer depth with a rapid fall off in concentration at the outer edge. The theoretical analysis based on a concentration flux model is in good qualitative agreement with observations.

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.

Turbulent collision of inertial particles: point-particle based, hybrid simulations and beyond

Point-particle based direct numerical simulation (PPDNS) has been a productive research tool for studying both single-particle and particle-pair statistics of inertial particles suspended in a turbulent carrier flow. Here we focus on its use in addressing particle-pair statistics relevant to the quantification of turbulent collision rate of inertial particles. PPDNS is particularly useful as the interaction of particles with small-scale (dissipative) turbulent motion of the carrier flow is mostly relevant. Furthermore, since the particle size may be much smaller than the Kolmogorov length of the background fluid turbulence, a large number of particles are needed to accumulate meaningful pair statistics. Starting from the relative simple Lagrangian tracking of so-called ghost particles, PPDNS has significantly advanced our theoretical understanding of the kinematic formulation of the turbulent geometric collision kernel by providing essential data on dynamic collision kernel, radial relative velocity, and radial distribution function. A recent extension of PPDNS is a hybrid direct numerical simulation (HDNS) approach in which the effect of local hydrodynamic interactions of particles is considered, allowing quantitative assessment of the enhancement of collision efficiency by fluid turbulence. Limitations and open issues in PPDNS and HDNS are discussed. Finally, on-going studies of turbulent collision of inertial particles using large-eddy simulations and particle- resolved simulations are briefly discussed.

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4:2244-2251; Int. J. Multiphase Flow 2000; 26:1583-1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes.

The influence of Saffman lift force on nanoparticle concentration distribution

The lift force on a spherical nanoparticle near a wall in micro/nanofluidics has not received sufficient attention so far. In this letter the concentration of 200 nm particles is measured at 0.25–2.0 m to a wall in a microchannel with pressure-driven de-ionized water flow pressure gradient 0–2000 kPa/m . The measured data show the influence of the lift force on the nanoparticle concentration distribution. By introducing the Saffman lift force into the Nernst–Planck equation near a wall, we find that the lift force is dominant at the range of 2particle radius, z is the distance from the wall).

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.

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.

MEASUREMENT OF SOLID SLURRY FLOW VIA CORRELATION OF ELECTROMAGNETIC FLOW METER, ELECTRICAL RESISTANCE TOMOGRAPHY AND MECHANISTIC MODELLING

The study presented here was carried out to obtain the actual solids flow rate by the combination of electrical resistance tomography and electromagnetic flow meter. A new in-situ measurement method based on measurements of the Electromagnetic Flow Meters (EFM) and Electrical Resistance Tomography (ERT) to study the flow rates of individual phases in a vertical flow was proposed. The study was based on laboratory experiments that were carried out with a 50 mm vertical flow rig for a number of sand concentrations and different mixture velocities. A range of sand slurries with median particle size from 212 mu m to 355 mu m was tested. The solid concentration by volume covered was 5% and 15%, and the corresponding density of 5% was 1078 kg/m(3) and of 15% was 1238 kg/m(3). The flow velocity was between 1.5 m/s and 3.0 m/s. A total of 6 experimental tests were conducted. The equivalent liquid model was adopted to validate in-situ volumetric solids fraction and calculate the slip velocity. The results show that the ERT technique can be used in conjunction with an electromagnetic flow meter as a way of measurement of slurry flow rate in a vertical pipe flow. However it should be emphasized that the EFM results must be treated with reservation when the flow pattern at the EFM mounting position is a non-homogenous flow. The flow rate obtained by the EFM should be corrected considering the slip velocity and the flow pattern.

The effect of food concentration on the life history of three types of Brachionus calyciflorus females

The effect of food concentration on the life history of three types of Brachionus calciflorus females (amictic, unfertilized mictic and fertilized mictic female) was studied with replicated individual cultures at 25 degrees and at four food concentrations (1.5, 3.0, 6.0 and 9.0 x 10(6) cells mL(-1)) of Scenedesmus obliquus. There were highly significant effects of both food concentration and female type, independently and in interaction on the duration of juvenile period of the rotifer, but neither a;ere the effects on the duration of post-reproductive period and mean life-span. The duration of juvenile period of unfertilized mictic female at the food concentration of 9.0 x 10(6) cells mL(-1) was the longest among all the food concentration-female type combinations. Both food concentration and female type influenced significantly the duration of reproductive period and the number of eggs produced by each type of female per life cycle, respectively. There was, however, no significant interaction between food level and female type. Among the three types of females, the number of eggs produced by an unfertilized mictic female was the largest. and that of a fertilized mictic female was the smallest.

Modeling and optimizing of high-concentration erbium-doped fiber amplifiers with consideration of ion-clusters

Starting from the modeling of isolated ions and ion-clusters, a closed form rate and power evolution equations for high-concentration erbium-doped fiber amplifiers are constructed. Based on the equations, the effects of the fraction of ion-clusters in total ions and the number of ions per cluster on the performance of high-concentration erbium-doped fiber amplifiers are analyzed numerically. The results show that the presence of the ion-clusters deteriorates amplifier performance, such as the signal power, signal gain, the threshold pump power for zero gain, saturated signal gain, and the maximum gain efficiency, etc. The optimum fiber length or other parameters should be modified with the ion-clusters being taken into account for the amplifiers to achieve a better performance. (c) 2007 Elsevier B.V. All rights reserved.

Ignition of single coal particle in a hot furnace under normal- and micro-gravity condition

An experimental study on ignition and combustion of single particles was conducted at normal gravity (1-g) and microgravity (l-g) for three high volatile coals with initial diameter of 1.5 and 2.0 mm, respectively. The non-intrusive twin-color pyrometry method was used to retrieve the surface temperature of the coal particle through processing the images taken by a color CCD camera. At the same time, a mathematical model considering thermal conduction inside the coal particle was developed to simulate the ignition process. Both experiments and modeling found that ignition occurred homogeneously at the beginning and then heterogeneously for the testing coal particles burning at l-g. Experimental results confirmed that ignition temperature decreased with increasing volatile content and increasing particle size. However, contradicted to previous studies, this study found that for a given coal with certain particle size, ignition temperature was about 50–80 K lower at l-g than that at 1-g. The model predictions agreed well with the l-g experimental data on ignition temperature. The criterion that the temperature gradient in the space away from the particle surface equaled to zero was validated to determine the commence of homogeneous ignition. Thermal conduction inside the particle could have a noticeable effect for determining the ignition temperature. With the consideration of thermal conduction, the critical size for the phase transient from homogeneous to heterogeneous is about 700 lm at ambient temperature 1500 K and oxygen concentration 0.23. 2009 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Comparison of aggregation behaviors between branched and linear block polyethers: MesoDyn simulation study

The comparison of aggregation behaviors between the branched block polyether T1107 (polyether A) and linear polyether (EO)(60)(PO)(40)(EO)(60) (polyether B) in aqueous solution are investigated by the MesoDyn simulation. Polyether A forms micelles at lower concentration and has a smaller aggregation number than B. Both the polyethers show the time-dependent micellar growth behaviors. The spherical micelles appear and then change to rod-like micelles with time evolution in the 10 vol% solution of polyether A. The micellar cluster appears and changes to pseudo-spherical micelles with time evolution in the 20 vol% solution of polyether A. However, the spherical micelles appear and change to micellar cluster with time evolution in the 20 vol% polyether B solution. The shear can induce the micellar transition of both block polyethers. When the shear rate is 1x10(5) s(-1), the shear can induce the sphere-to-rod transition of both polyethers at the concentration of 10 and 20 vol%. When the shear rate is lower than 1x10(5) s(-1), the huge micelles and micellar clusters can be formed in the 10 and 20 vol% polyether A systems under the shear, while the huge micelles are formed and then disaggregated with the time evolution in the 20 vol% polyether B system.

Shell model study of single-particle and collective structure in neutron-rich Cr isotopes

The structure of neutron-rich Cr isotopes is systematically investigated by using the spherical shell model. The calculations reproduce well the known energy levels for the even-even Cr52-62 and odd-mass Cr53-59 nuclei, and predict a lowering of excitation energies around neutron number N = 40. The calculated B(E2; 2(1)(+) -> 0(1)(+)) systematics shows a pronounced collectivity around N = 40; a similar characteristic behavior has been suggested for Zn and Ge isotopes. Causes for the sudden drop of the 9/2(1)(+) energy in Cr-59 and the appearance of very low 0(2)(+) states around N = 40 are discussed. We also predict a new band with strong collectivity built on the 0(2)(+) state in the N = 40 isotope Cr-64.

The electrochemical behavior of a system with a limited number of molecules

In this paper, based on Einstein relationship between diffusion and random walk, the electrochemical behavior of a system with a limited number of molecules was simulated and explored theoretically. The transition of the current vs time responses from discrete to continuous was clearly obtained as the number of redox molecules increased from 10 to 10(6).

Independence of the brittle-ductile transition from the rubber particle size for impact-modified poly(vinyl chloride)

A series of acrylic impact modifiers (AIMS) with different particle sizes ranging from 55.2 to 927.0 nm were synthesized by seeded emulsion polymerization, and the effect of the particle size on the brittle-ductile transition of impact-modified poly(vinyl chloride) (PVC) was investigated. For each AIM, a series of PVC/AIM blends with compositions of 6, 8, 10, 12, and 15 phr AIM in 100 phr PVC were prepared, and the Izod impact strengths of these blends were tested at 23 degrees C. For AIMs with particle sizes of 55.2, 59.8, 125.2, 243.2, and 341.1 nm, the blends fractured in the brittle mode when the concentration of AIM was lower than 10 phr, whereas the blends showed ductile fracture when the AIM concentration reached 10 phr. It was concluded that the brittle-ductile transition of the PVC/AIM blends was independent of the particle size in the range of 55.2-341.1 nm. When the particle size was greater than 341.1 nm, however, the brittle-ductile transition shifted to a higher AIM concentration with an increase in the particle size. Furthermore, the critical interparticle distance was found not to be the criterion of the brittle-ductile transition for the PVC/AIM blends.