Modeling and numerical simulation of flow field in three types of standard new design cyclone separators

In this paper, numerical simulation of flow field in three types of standard new design cyclone separators namely 1D2Dn, 1D3Dn and 2D2Dn are investigated. In these standard cyclones, the length of cylindrical top part of the body is equal to 1, 1 and 2 times of the body diameter, respectively; and the length of the cylindrical bottom part is 2, 3 and 2 times of the body diameter. The new design cyclone is based on the idea of improving cyclone collection efficiency and pressure drop by increasing the vortex length. The Eulerian-Lagrangian computational procedure is used to predict particles tracking in the cyclones. The velocity fluctuations are simulated using the Discrete Random Walk (DRW). Results show that among the three standard new design cyclones, cyclone 2D2Dn has the highest efficiency followed by 1D3Dn one with about only 2% lower efficiency. Cyclone 1D2Dn possesses the lowest efficiency among all. Similarly, the highest pressure drop occurs in cyclone 2D2Dn. Cyclones 1D3Dn and 1D2Dn followed 2D2Dn one with a marked difference of about 20%. In result section, the details of the flow field including velocity, pressure contours, turbulence, velocity vectors and particle trajectory will be presented.     

Numerical study of flow field in new design cyclones with different wall temperature profiles: Comparison with conventional ones

In this paper, numerical study of flow field in the new design cyclones with five different wall temperature profiles are investigated. The new design cyclone is based on the idea of improving cyclone collection efficiency and pressure drop by increasing the vortex length. In this paper, the five wall temperature profiles are as follows: (A) cooling with uniform distribution, (B) without temperature change, (C) heating with uniform distribution, (D) incremental linear heating, (E) reduction linear heating. Results are compared in new design and conventional cyclones. The Reynolds averaged Navier–Stokes equations with Reynolds stress turbulence model (RSM) are solved. The Eulerian-Lagrangian computational procedure is used to predict particles tracking in the cyclones. The velocity fluctuations are simulated using the Discrete Random Walk (DRW).Results show that generally, heating the bottom zone of the cyclones can improve the collection efficiency and reduce the pressure drop while heating the top zone of the cyclones marginally affects the flow field. Moreover, cooling the cyclones reduces the efficiency and causes a higher pressure drop. Among five different wall temperature profiles, C and E profiles can increase the efficiency about 8% and profile C reduces the pressure drop by about 9%. The mentioned values in different conditions including particle diameter, flow rate, etc. can be different.     

The impact of cylinder vortex stabilizer on fluctuating turbulence characteristics of a cyclone separator based on Large Eddy Simulation

The main purpose of present study is to comprehensively clarify the impact of cylinder vortex stabilizer on fluctuating turbulence structure of a Stairmand cyclone separator on basis of Large Eddy Simulation. The cylinder vortex stabilizer is easy and could be applied to any existing cyclone model without any major replacement. This novel modification in cyclone body is considered to alleviate the negative effect of entrainment of particles from the ash hopper and swing of the vortex end in swirling flow. The numerical simulations were conducted based on Stairmand cyclone separator and three new models with variation of vortex stabilizer length and diameter. The results showed that the cylinder vortex stabilizer could enhance flow instability and improve fluctuating turbulence structure to some extent. It is confirmed that cylinder vortex stabilizer could significantly reduce the tangential velocity in the inner quasi-forced vortex region of the cyclones. Comparing with Stairmand cyclone, the swirling first and second peak frequency of cyclone model with vortex stabilizer (Length L/D: 6.5, diameter d/D: 0.12) have been confirmed to get considerable reduction of 11.54% and 10.86%, respectively. This modified cyclone model is comparatively better for enhancement of flow stability, providing about 18.4% maximum reduction of normalized flow angle, 24.8% of rotational kinetic energy in dust collector and 14.2% in the main body of cyclone.     

Near-wake flow characteristics of two side-by-side circular cylinders close to a wall

Flow characteristics in the near wake of two identical side-by-side circular cylinders located close to a fully developed turbulent boundary layer are investigated experimentally using the particle image velocimetry (PIV) and the pressure sensor. The Reynolds number based on the cylinder diameter (D) is 1,696, the boundary layer thickness is 6.6 D, the cylinder center-to-center spacing (T) is varied from T/D = 1 to 1.906, and the gap spacing between the lower cylinder and the wall (G) is varied from G/D = 0 to 1.811. To study the effects of changing the gap ratios of T/D and G/D on the wake flow, various wake characteristics such as averaged streamlines, Reynolds stress and vorticity contours as well as other key flow features including the length scales and the Strouhal number are investigated for different ratios of T/D and G/D. According to these wake characteristics, five basic flow patterns have been identified.     

The influence of projectile hardness on ballistic performance

The ballistic performance of 17 penetrator materials, representing 5 distinct steel alloys treated to various hardnesses along with one tungsten alloy, has been investigated. Residual lengths and velocities, as well as the ballistic limit velocities, were determined experimentally for each of the alloy types for length-to-diameter (L/D) ratio 10 projectiles against finite-thick armor steel targets. The target thickness normalized by the projectile diameter (T/D) was 3.55. For some of the projectile types, a harder target, with the same thickness, was also used. It was found that the ballistic limit velocity decreases significantly when the projectile hardness exceeds that of the target. Numerical simulations are used to investigate some of the observed trends. It is shown that the residual projectile length is sensitive to projectile hardness; the numerical simulations reproduce this experimental observation. However, the observed trend in residual velocity as a function of projectile hardness is not reproduced in the numerical simulations unless a material model is invoked. It is assumed that the plastic work per unit volume is approximately a constant, that is, there is a trade off between strength and ductility. Using this model, the numerical simulations reproduce the experimentally observed trend.     

Vortex ring formation for low Re numbers

The dynamics of formation and evolution of vortex rings with low Reynolds numbers created in a piston-cylinder arrangement are studied. The ratio of the piston displacement L _m to the nozzle diameter D ₀ determines the vortex size and evolution. Experiments with different conditions are presented: translation velocity of the piston and stroke ratio L _m /D ₀ for 150 < Re < 260. Measurements of the 2D velocity field were obtained with a PIV technique. The vortex circulation was computed considering a vortex identification scheme (Q criterion). The results show that there is a critical value of L _m /D ₀ above which the circulation inside the vortex cannot increase and remains constant. For the Reynolds numbers studied, we found that the limit stroke ratio is 4 ≤ L _m /D ₀ ≤ 6. As Re decreases, the vortices become “thicker”; therefore, they are able to accumulate more vorticity and increase their circulation.     

Polling systems with batch service

Motivated by applications in production and computer-communication systems, we study an N-queue polling system, consisting of an inner part and an outer part, and where products receive service in batches. Type-i products arrive at the outer system according to a renewal process and accumulate into a type-i batch. As soon as D _i products have accumulated, the batch is forwarded to the inner system where the batch is processed. The service requirement of a type-i batch is independent of its size D _i . For this model, we study the problem of determining the combination of batch sizes ${\vec{D}^{({\rm opt})} }$ that minimizes a weighted sum of the mean waiting times. This model does not allow for an exact analysis. Therefore, we propose a simple closed-form approximation for ${\vec{D}^{({\rm opt})}}$ , and present a numerical approach, based on the recently proposed mean waiting-time approximation in Boon et?al. (Perform Eval 68, 290–306, 2011). Extensive numerical experimentation shows that the numerical approach is slightly more accurate than the closed-form solution, while the latter provides explicit insights into the dependence of the optimal batch sizes on the system parameters and into the behavior of the system. As a by-product, we observe near-insensitivity properties of ${\vec{D}^{({\rm opt})}}$ , e.g. to higher moments of the interarrival and switch-over time distributions.     

Deagglomeration of airborne nanoparticles in a decelerating supersonic round jet

The size and the agglomeration state of nanoparticles (NPs) play a very important role in nanoparticle applications due to the effect of the dispersion level of NPs on their toxicity, pharmaceutical activity, and catalytic activity. In this study, a novel two-stage nanoparticle disperser system was proposed. We focused on the effects of the aspect ratio (L/D, length (L) and diameter (D)) of the round jet (a straight pipe) on the nanoparticle size distribution (PSD) and we attempted to explain how the jet wakes induced the NP deagglomeration based on the numerical results. The median diameters, mode diameters, and geometric mean diameters (GMD) of the PSDs that were obtained using the second-stage were significantly smaller than those achieved using the first-stage disperser alone. Overall, the GMD, mode diameters, and median diameters could all reach sub-70 nm, which was approaching close to the NP primary diameter (30–40 nm). The particle sizes decreased with the increasing dimensionless nozzle diameter (D* = D/D_p × 10^?3, D_p represents the particle diameter). The increase in L/D for a constant D hardly affected the PSDs but led them to be closer to a standard lognormal distribution. Regarding the NP traveling path and its considerable inertia as it penetrated the Mach disc, these particular field distributions of the jet-wake might explain why the suspended sub-micro NP agglomerates could be highly deagglomerated. Namely, the nano NP agglomerates would suffer the severe shear stress around the Mach disc.     

Design,preparation and properties of online-joints of C/SiC–C/SiC with pins

This work is aimed at providing a new joining technology for C/SiC composites and investigating the influence of drilling holes, hole distribution (including ratios of edge distance to diameter (E/D), width to diameter (W/D) and hole distance to diameter (H/D)) and the number of applied pins on the mechanical properties of C/SiC substrates and joints. The mechanical testing results show that drilling holes and hole distribution greatly affects the mechanical properties of C/SiC substrates but when adopting an optimized design principle (E/D ⩾ 3, W/D ⩾ 3 and H/D ⩾ 3) the effect could be neglected. 1D C/SiC pins with higher shearing strength (107.2 MPa) are more suitable to join the substrates. With the increase of pins (1, 2 and 4), the bearing loads of the joints increase almost linearly, and the reliability of joints is also improved in that the fracture mode changes from the interlayer damage to the substrate rupture. Besides, the joining process generates uniform and dense joining layer (composition of ZrC and SiC) and a strong bonding without obvious interface.     

Effect of a protruding rod-supported disk on the drag of a nose-controlled cylinder

The drag C _x of a cylinder of diameter D with a front protruding disk supported on a rod of length l has been studied as a function of the relative distance l/D under the conditions of high (supersonic) flight velocities. It is established that the optimum (minimum) drug C _x exists, the value of which agrees with the results of numerical simulations.     

Numerical simulation of square cyclones in small sizes

The idea of using square cyclones was first introduced in early 1990s because of some problems of big conventional (round) cyclones in Circulating Fluidized Bed (CFB) industries, such as huge volume and long start-stop time of cyclones. Now there is this question, in spite of the main reason for making square cyclones, how square cyclones behave in small sizes.In this paper, two small cyclones with the same hydraulic diameter, which one is square and the other one is round, are numerically compared. The Reynolds averaged Navier–Stokes equations with Reynolds Stress Turbulence Model (RSTM) are solved by use of the finite volume method based on the SIMPLE pressure correction algorithm in the computational domain. The Eulerian–Lagrangian computational procedure is used to predict particles tracking in the cyclones. The velocity fluctuations are simulated using the Discrete Random Walk (DRW).The results show that the pressure drop in small square cyclone is less than the pressure drop in small round one. Also at each flow rate collection efficiency of small square cyclone is less than round one, but by increasing flow rate this difference decreases.     

Experimental and numerical failure analysis of carbon/epoxy laminated composite joints under different conditions

The aim of this study is to investigate the effects of preload moment, moisture and interference-fit on bearing strength and failure mode in pin-jointed and bolted carbon–epoxy plates which were subjected to a traction force. Two different geometrical parameters, end distance to pin diameter ratio (E/D) and width to pin diameter ratio (W/D), were considered experimentally. E/D and W/D ratios were selected ranging from 1 to 4 and from 2 to 4, respectively. The test results showed that the ultimate failure loads were directly affected by the geometrical parameters, preload moments and interference-fit. Furthermore, for the specimens in wet condition, failure loads increased with preload moments while failure loads stayed almost the same in non-preloaded specimens. Finite element models of the specimens were also developed by using ANSYS software. Contact and stress analyses of the models were carried out for both free and interference-fit. Tsai-Wu criterion was used to determine bearing strength corresponding to first failure load. It was concluded that the numerical results are in good agreement with experimental results.     

Accurate quasi static capacitance for abrupt homojunction under forward and reverse polarization

In this work, we present a new approach to derive the capacitance–voltage characteristic for an abrupt homojunction with uniform doping (N _A in p-region and N _D in n-region) under forward and reverse polarization. Under thermal equilibrium conditions, we show that it is possible to obtain analytically the exact capacitance–voltage characteristic without any simulation for the symmetric case (N _A ?=?N _D ). We also propose a model of the total capacitance under forward and reverse polarization for the asymmetric case ( $N_{\bf A} \ne N_{\bf D}$ ).     

Shape effect on weakly nonlinear elliptical composites

The electric field responses of two types of weakly nonlinear dielectric composites consisting of elliptic cylindrical inclusions, one with identical shape and the another with distributed shapes, randomly embedded in the linear host media in the dilute limit are investigated. The dielectric property of the inclusions is that the relation between the displacement (D) and electric (E) fields obey the form D = ?E + χ∣E∣^βE where β is a nonlinear integer exponent and ? ? χ∣E∣^β. By using the decoupling approximation, the effective nonlinear susceptibility (χ_e) is determined and analyzed for varying the aspect ratios and the shape distribution parameters for the composites with identical and distributed inclusion shapes, respectively. In addition, the exact analytic result of χ_e for the elliptical composites with distributed inclusion shapes for the case of β = 2 is derived in this article.     

Frequency and Size Dependences of Superfluidity in?Low-Dimensional 4He Fluids

We have studied superfluidity of ⁴He fluids in two- and one-dimensional states. In the 2D state of the ⁴He films on flat substrate, superfluidity was observed in the normal fluid state above the 2D Kosterlitz-Thouless temperature at high measurement frequencies. The superfluidity in 2D also depends on the system size, e.g. pore diameter of porous glasses and grain size of powder. The 1D state was realized for ⁴He fluid nanotubes formed in 1D nanopores 1.8?C4.7 nm in diameter and about 300 nm in length. Superfluidity in the 1D state was observed by the torsional oscillator experiment. The results are qualitatively well reproduced by the Monte Carlo calculation for a classical XY spin system modeled on the present ⁴He nanotubes. Although the 1D state is in the normal fluid state at any finite low temperatures due to the 1D phonon fluctuation, the superfluid frequency shift ??f of the oscillator can be observed. Above a temperature, ??f decreases due to another kind of 1D thermal fluctuation of which excitations destroy the phase coherence in the nanotubes. The excitations depend on the tube length as well as the tube diameter.     

ESR and X-ray Investigations of Deuterium Atoms and Molecules in Impurity-Helium Solids

Impurity-helium solids created by injecting deuterium atoms and molecules into superfluid ⁴He have been studied via electron spin resonance (ESR) and x-ray diffraction methods. We measured the g-factor, the hyperfine constant and the spin-lattice relaxation time of D atoms in D-D₂-He solids. These measurements show that D atoms are mainly stabilized in D₂ clusters. Using an x-ray method we found the size of D₂ clusters to be ～90Å in diameter and the densities of D₂ molecules in the samples to be of order 2.5?10²¹ cm^?3 . The highest average concentration of D atoms achieved in D-D₂-He solids was ～1.5?10¹⁸ cm^?3 . The local concentrations of D atoms within D₂ clusters is found to be large (～2?10¹⁹ cm^?3).     

The use of smoothed particle hydrodynamics for simulating crystal growth from solution

Most numerical simulations for crystal growth processes have been performed using the finite volume and finite element techniques. Both of these methods require adaptive meshes to track the growth and dissolution interfaces. In addition, in the growth of ternary systems the solid and liquid phases must be solved simultaneously, which requires iterations at the interfacial boundaries between solid and liquid regions. Such iterations slow down simulations tremendously, and also lead to numerical instabilities. In order to address these issues, an alternative mesh-free, Lagrangian method, known as smoothed particle hydrodynamics (SPH) has been investigated. For the simulation, the liquid phase diffusion (LPD) growth of Si_xGe_1−x has been selected, since both experimental and finite volume simulation results were available for comparison. This comparison between finite volume and SPH simulations shows that although SPH has the potential to accurately model the crystal growth process, the number of SPH particles required for accurate predictions is upwards of 60,000 particles for the Reynolds number of the LPD system. This high number of particles translates to a computational time of approximately six times longer than the equivalent finite volume simulations. However, future improvements made to the relatively young SPH method may overcome such computational difficulties.     

Two-phase pressure drop in return bends: Experimental results for R-410A

This study presents 238 pressure drop data points measured for two-phase flow of R-410A in horizontal return bends. The tube diameter (D) varies from 7.90 to 10.85 mm and the curvature ratio (2R/D) from 3.68 to 4.05. The mass velocity ranges from 179 to 1695 kg m⁻² s⁻¹ and the saturation temperatures from 4.6 °C to 20.7 °C. Preliminary tests show that the recovery length necessary for a correct pressure drop measurement downstream of the return bend is less than 20D, for the experimental conditions covered in this study. The singular pressure drop is determined by subtracting the regular pressure drop in straight tube from the total pressure drop. The experimental data are compared against four available correlations found in the literature. The present experimental database for the return bend pressure drop is presented in the Appendix A.     

Vortex Diffusivity and Core Diameter of 2D Superfluid in 4He Films on Gold and H2 Substrates

The two-dimensional (2D) ⁴He fluid films show the Kosterlitz-Thouless (KT) transition where pairing and unpairing of the 2D vortices play an important role. However, the vortex properties (the diffusion constant D, the core diameter a ₀) have not been precisely obtained for various conditions. Here, we accurately determined the parameter D/a ₀ ² by the high frequency dependence of the superfluid onset up to 180 MHz for the submonolayer ⁴He fluid films adsorbed on gold and H₂ (3.3 layers) preplated on gold, respectively. The superfluid onset coverage changes from 1.6 (gold) to 0.5 layers (H₂), which clearly indicates the large difference of the adsorption potential. The parameter D/a ₀ ² , on the other hand, has the same value for the coverages with the same KT temperature T _KT. This suggests that the vortex diffusions on both substrates have the largest value D～?/m in the quantum limit. The core diameter a ₀ was estimated to be the same magnitude as the de Broglie wavelength at T _KT between 0.1 and 0.9 K.     

Vortex ring propagation and interactions studies

This paper experimentally investigates the vortex ring propagation and interactions with thin cylindrical and flat surfaces. Dye-based visualization technique is adopted for the interaction studies. Vortex rings are generated from a circular nozzle of 19 mm diameter with the stroke length ratios (length of the fluid slug to nozzle diameter, L_N/D_N) of 1 to 5, and ejection velocities in the range of 0.05 to 0.2 m/s. Vortex interaction studies are carried out with two different bodies; firstly, with the circular cylinders having the diameters of 0.2, 0.6, 1.5 and 2.5 mm, and secondly with a flat solid surface. Results indicate that the trails in the vortex ring start following at L_N/D_N = 4. The influence of the initial velocity is found to be insignificant on the vortex ring diameter, however, found to depend on stroke length ratio. Vortex-cylinder interaction studies indicated that vortex velocity decreases with increase in cylinder diameter after the interaction. Reconnections of vortex rings are observed in lower cylinder diameter cases. In case of vortex ring interaction with the flat surface, stretching of the vortex core is observed leading to a considerable increase in the vortex ring diameter.