Transient heat transfer behaviors in cylindrical porous beds at relatively large Reynolds numbers

The transient behavior of heat transfer in a cylindrical porous bed was examined experimentally under various factors (flow rate, diameter of spherical solid particle, temperature of flowing fluid and physical properties of porous bed). In these factors, it was understood that especially flow rate and the diameter of the particled have important role in evaluating the transient behavior of heat transfer in the porous bed. That is, as the flow rate and the diameter of the particle under a constant diameterD of the cylindrical bed are increased, mean heat transfer coefficient between flowing fluid and the solid particles is increased and the time period to reach a thermally steady state is decreased. The useful experimental correlation equations of mean heat transfer coefficient and the time period to reach the steady state were derived with the functional relationships of Nusselt numberNu _d =f(d/D, Reynolds numberRe _d ) and Fourier numberFo =f (modified Prandtl numberPr*, d/D, Re _d ).     

Some remarks on “Drag coefficients of a slowly moving sphere in Non-Newtonian fluids”

Viscoelastic solutions were ejected vertically downwards into air and various Newtonian fluids. The measured swell increased significantly when ejected into a liquid rather than air. The observed increase is considered a result of both bouyancy and drag forces on the solution. The following dimensions expression relating the ratio of the swell diameter in liquid and air D_L/D_A to the elastic shear compliance of the ejected solution J_e was experimentally observed.(D_L/D_A)⁶-1=30(Δ?/?_s)^{?$\text{1}\text{2}$}([g²η²_N?_s]$\text{1}\text{3}$J_e)^{$\text{3}\text{5}$}, where Δ? is the density difference between the extruded and Newtonian fluid, ?_s is the solution density, g is the gravitational constant, and η_N is the Newtonian fluid viscosity. Thus with this expression a simple extrudate swell technique exists to estimate the elastic shear compliance of a viscoelastic solution.     

Flow fields and packing states in the discharge flow of noncircular particles—A SIPHPM simulation

Although the discharge flow of spherical materials has been extensively explored, the effect of particle shape on discharge is still poorly understood. The present work explores the two-dimensional discharge flow fields of noncircular particles using the soft-sphere-imbedded pseudo-hard particle model method. Rectangular particles having different aspect ratios (R_a = 1, 1.5, 2–5) and regular polygonal particles having different numbers of sides (N_s = 3–8, 10) are discharged through hopper beds having different orifice widths (D_i = 40, 70.77, 99.13, 125.74, 151.13 mm). The discharge rates of differently shaped particles in different beds are consistent with Beverloo’s relation. Moreover, the flow fields are computed and evaluated to study the effects of R_a, N_s, and D_i on particle discharge. The characteristics of particle–particle connections in the discharge process are evaluated according to the temporal evolution and spatial distribution of the contact points. Additionally, the effect of the initial packing on the discharge profile is investigated. The findings help clarify the discharge of noncircular particles.     

Experimental investigation on flow asymmetry in solid entrance region of a square circulating fluidized bed

To study the influence of back feeding particles on gas-solid flow in the riser, this paper investigated the flow asymmetry in the solid entrance region of a fluidized bed by particle concentration/velocity measurements in a cold square circulating fluidized beds （CFB）. The pressure drop distribution along the riser and the saturation carrying capacity of gas for Geldart-B type particles were first analyzed. Under the condition of u0 = 4 m/s and Gs = 21 kg/（m^2 s）, the back feeding particles were found to penetrate the lean gas-solid flow near the entrance （rear） wall before reaching the opposite （front） wall, thus leading to a relatively denser region near the front wall in the bottom bed. Higher solid circulation rate （u0 =4 m/s, Gs = 33 kg/（m^2 s）） resulted in a higher particle concentration in the riser. However the back feeding particles with higher momentum increased the asymmetry of the particle concentration/velocity profile in the solid entrance region. Lower air velocity （u0 =3.2 m/s） and Gs =21 kg/（m2 s）, beyond the saturation carrying capacity of gas, induced an S-shaped axial solid distribution with a denser bottom zone. This limited the penetration of the back feeding particles and forced the flnidizing air to flow in the central region, thus leading to a higher solid holdup near the rear wall. Under the conditions of uo = 4 m/s and Gs = 21 kg/（m^2 s）, addition of coarse particles （dp= 1145 μm） into the bed made the radial distribution of solids more symmetrical.     

Large strain deformation of polycrystalline metals at low homologous temperatures

Large strain compression data (true strains to about ?3.0) are presented for polycrystalline α U and α Fe at room temperature. The results, together with other published data at low homologous temperatures (≈0.2 T_m), where T_m is the absolute melting temperature, suggest that a steady-state flow stress σ_s is approached after extensive strain-hardening, α U exhibits a very high strain-hardening rate, with σ_s ≈ 2900 MPa (420 ksi) indicating that cold-working is a very potent method of strengthening this metal. All the data evaluated can be fit by the stress-strain relation σ = σ_s? exp (?(Nε)^p)(σ_s? σ_y), where σ_y is the yield stess, p is a constant equal to a for the metals analyzed, N is a constant associated with the strain-hardening characteristics of a material, σ is true stress, and ε is true strain.     

Particle resolved direct numerical simulation of a liquid–solid fluidized bed: Comparison with experimental data

Particle-resolved direct numerical simulations of a 3-D liquid–solid fluidized bed experimentally investigated by Aguilar-Corona (2008) have been performed at different fluidization velocities (corresponding to a range of bed solid volume fraction between 0.1 and 0.4), using Implicit Tensorial Penalty Method. Particle Reynolds number and Stokes number are O(100) and O(10), respectively. In this paper, we compare the statistical quantities computed from numerical results with the experimental data obtained with 3-D trajectography and High Frequency PIV. Fluidization law predicted by the numerical simulations is in very good agreement with the experimental curve and the main features of trajectories and Lagrangian velocity signal of the particles are well reproduced by the simulations. The evolution of particle and flow velocity variances as a function of bed solid volume fraction is also well captured by the simulations. In particular, the numerical simulations predict the right level of anisotropy of the dispersed phase fluctuations and its independence of bed solid volume fraction. They also confirm the high value of the ratio between the fluid and the particle phase fluctuating kinetic energy. A quick analysis suggests that the fluid velocity fluctuations are mainly driven by fluid–particle wake interactions (pseudo-turbulence) whereas the particle velocity fluctuations derive essentially from the large scale flow motion (recirculation). Lagrangian autocorrelation function of particle fluctuating velocity exhibits large-scale oscillations, which are not observed in the corresponding experimental curves, a difference probably due to a statistical averaging effect. Evolution as a function of the bed solid volume fraction and the collision frequency based upon transverse component of particle kinetic energy correctly matches the experimental trend and is well fitted by a theoretical expression derived from Kinetic Theory of Granular Flows.     

Wave propagation in alternating solid and fluid layers

The exact solution for harmonic plane wave propagation through periodically stratified media composed of alternating elastic solid and ideal fluid layers is derived for all frequencies and wavenumbers. The modal structure and the structure of the pass and stop bands is shown in the plane of slowness parallel to the stratification, s₁, and frequency, ω. In the low frequency limit, the slowness component perpendicular to the stratification, s₃, is found as a function of s₁, defining slowness surfaces in the s₁ and s₃ plane. At any angle of phase propagation through the medium (except perpendicular) there are two modes, a fast wave for which the motion in the solid and fluid layers is in phase and a slow wave for which the motion is 180° out of phase. An equivalent acoustic medium is found which has an anisotropic density tensor and a nonlocal constitutive relation between pressure and dilatation whose terms are fourth derivatives in time and in the space coordinate parallel to the stratification. This equivalent homogeneous anisotropic medium exhibits the indentical wave behavior as the layered solid fluid medium in the long wavelength limit as it has the same slowness surfaces.     

The sliding interface crack with friction between elastic and rigid bodies

The paper analyzes the frictional sliding crack at the interface between a semi-infinite elastic body and a rigid one. It gives solutions in complex form for non-homogeneous loading at infinity and explicit solutions for polynomial loading at the interface. It is found that the singularities at the crack tips are different and that they are related to distinct kinematics at the crack tips. Firstly, we postulate that the geometry of the equilibrium crack with crack-tip positions b and a is determined by the conditions of square integrable stresses and continuous displacement at both crack tips. The crack geometry solution is not unique and is defined by any compatible pair (b,a) belonging to a quasi-elliptical curve. Then we prove that, for an equilibrium crack under given applied load, the “energy release rate” G_tip, defined at each crack tip by the J_ε-integral along a semi-circular path, centered at the crack tip, with vanishing radius ε, vanishes. For arbitrarily shaped paths embracing the whole crack, with end points on the unbroken zone, the J-integral is path-independent and has the significance of the rate, with respect to the crack length, of energy dissipated by friction on the crack.     

Particle image velocimetry in a variable density flow: application to a dynamically evolving microburst

A fondamental difficulty in the experimental study of gravity-driven flows using particle image velocimetry (PIV) and other optical diagnostic techniques is the problem associated with variations in thé refractive index within the fluid. This paper discusses a method by which the refractive indices of two fluids are matched while maintaining density differences of up to 4%. Aqueous solutions of glycerol and potassium phosphate are used to achieve precise index matching in the presence of mixed and unmixed constituents. The effectiveness of the method is verified in a PIV study of a laboratory-scale model of an atmospheric microburst where planes of two-dimensional velocity vectors are obtained in thé evolving flow field.This work was sponsored by the National Science Foundation under grant CTS-9209948. We also thank TSI, Inc. for the use of its facility     

Application of CFD-DEM to the study of solid exchange in a dual-leg fluidized bed

The CFD-DEM model was developed to simulate solid exchange behavior between two half beds in a bench-scale two-dimensional dual-leg fluidized bed （DL-FB）. Power spectrum density （PSD） analysis was applied to obtain the dominant frequency （F） of the simulated differential particle number （APLR） between the two half beds. Effects of fluidization velocity （u） and bed material inventory （H） on the solid exchange behavior were studied using the CFD-DEM model. Not only snapshots of the simulated particle flow patterns using the OpenGL code but also the dominant frequency of APLR was similar to the experimental results. The simulation results show that higher fluidization velocity assists the exchange of more particles between the two half beds, but the dispersion of clusters on the bed surface into single particles decreases the cluster exchange frequency. A greater bed material inventory results in more intense cluster exchange. The cluster exchange frequency decreases with an increase of the bed material inventory.     

Particle image velocimetry study of the shock-induced single mode Richtmyer–Meshkov instability

Full field particle image velocimetry (PIV) measurements are obtained for the first time in Richtmyer–Meshkov instability shock tube experiments. The experiments are carried out in a vertical shock tube in which the light gas (air) and the heavy gas (SF₆) flow from opposite ends of the shock tube driven section and exit through narrow slots at the interface location. A sinusoidal perturbation is given to the interface by oscillating the shock tube in the horizontal direction. Richtmyer–Meshkov instability is then produced by the interaction with a weak shock wave (M _s  = 1.21). PIV measurements are obtained by seeding the flow with 0.30 μm polystyrene Latex spheres which are illuminated using a double-pulsed Nd:YAG laser. PIV measurements indicate the vorticity to be distributed in a sheet-like distribution on the interface immediately after shock interaction and that this distribution quickly rolls up into compact vortices. The integration of the vorticity distribution over one half wave length shows the circulation to increase with time in qualitative agreement with the numerical study of Peng et al. (Phys. Fluids, 15, 3730–3744, 2003).     

Particle convective heat transfer near the wall in a supercritical water fluidized bed by single particle model coupled with CFD-DEM

Supercritical water fluidized bed (SCWFB) is a promising reactor to gasify biomass or coal. Its optimization design is closely related to wall-to-bed heat transfer, where particle convective heat transfer plays an important role. This paper evaluates the particle convective heat transfer coefficient (h_pc) at the wall in SCWFB using the single particle model. The critical parameters in the single particle model which is difficult to get experimentally are obtained by the computational fluid dynamics-discrete element method (CFD-DEM). The contact statistics related to particle-to-wall heat transfer, such as contact number and contact distance, are also presented. The results show that particle residence time (τ), as the key parameter to evaluate h_pc, is found to decrease with rising velocity, while increase with larger thermal boundary layer thickness. τ follows a gamma function initially adopted in the gas–solid fluidized bed, making it possible to evaluate h_pc in SCWFB by a simplified single particle model. The theoretical predicted h_pc tends to increase with rising thermal gradient thickness at a lower velocity (1.5 U_mf), while first decreases and then increases at higher velocity (1.75 and 2 U_mf). h_pc occupies 30%–57% of the overall wall-to-bed heat transfer coefficient for a particle diameter of 0.25 mm. The results are helpful to predict the overall wall-to-bed heat transfer coefficient in SCWFB combined with a reasonable fluid convective heat transfer model from a theoretical perspective.     

A note on a boundary condition for spouted beds

The boundary condition, zero solids pressure at the top of a particle bed of maximum spoutable height, Hm, is shown to eliminate any resort to empiricism in the derivation of the fluid velocity in the annulus of a spouted bed for which both viscous and inertial effects are taken into account. The same boundary condition fails when applied to a spouted bed for which the bed height H 〈 Hm, especially when H 〈 0.8Hm.     

The size effect on void growth in ductile materials

We have extended the Rice-Tracey model (J. Mech. Phys. Solids 17 (1969) 201) of void growth to account for the void size effect based on the Taylor dislocation model, and have found that small voids tend to grow slower than large voids. For a perfectly plastic solid, the void size effect comes into play through the ratio εl/R₀, where l is the intrinsic material length on the order of microns, ε the remote effective strain, and R₀ the void size. For micron-sized voids and small remote effective strain such that εl/R₀?0.02, the void size influences the void growth rate only at high stress triaxialities. However, for sub-micron-sized voids and relatively large effective strain such that εl/R₀>0.2, the void size has a significant effect on the void growth rate at all levels of stress triaxiality. We have also obtained the asymptotic solutions of void growth rate at high stress triaxialities accounting for the void size effect. For εl/R₀>0.2, the void growth rate scales with the square of mean stress, rather than the exponential function in the Rice-Tracey model (1969). The void size effect in a power-law hardening solid has also been studied.     

Resonant non-linear vibrations of continuous systems—II Damped and transient behavior

The method presented in Part I is extended to cover the damped and transient behavior of nonlinear systems described by equations of the form u_ss?u_yy?εf(u,u_y, u_yys,…,s) = 0. The method is presented by application to the equation u_ss?u_yy?εu_yys?εu²_yu_yy= 0. Similar to the undamped case it is again shown that the PDE requires an infinite number of periodicity conditions to correctly characterize the resonant region. However, damping eliminates some of the branches of the amplitude-frequency spectrum of the undamped case. In fact, for ε = 0.25 all but the outermost branch disappear. A method of multiple time scales is presented for the study of the transient behavior and the stability of the branches for steady vibrations. The stability analysis yields an interior stable point in the amplitude-frequency spectrum which has no analog in the Duffing equation. Finally via the multiple scale procedure in the spirit of the early work of Zabusky and Kruskal one obtains forced Burgers and Korteweg-de Vries equations on a finite interval.     

An investigation of matched index of refraction technique and its application in optical measurements of fluid flow

Optical distortions caused by non-uniformities of the refractive index within the measurement volume is a major impediment for all laser diagnostic imaging techniques applied in experimental fluid dynamic studies. Matching the refractive indices of the working fluid and the test section walls and interfaces provides an effective solution to this problem. The experimental set-ups designed to be used along with laser imaging techniques are typically constructed of transparent solid materials. In this investigation, different types of aqueous salt solutions and various organic fluids are studied for refractive index matching with acrylic and fused quartz, which are commonly used in construction of the test sections. One aqueous CaCl₂·2H₂O solution (63 % by weight) and two organic fluids, Dibutyl Phthalate and P-Cymene, are suggested for refractive index matching with fused quartz and acrylic, respectively. Moreover, the temperature dependence of the refractive indices of these fluids is investigated, and the Thermooptic Constant is calculated for each fluid. Finally, the fluid viscosity for different shear rates is measured as a function of temperature and is applied to characterize the physical behavior of the proposed fluids.     

Resonant non-linear vibrations in continuous systems—I. Undamped case

A method is presented for obtaining periodic solutions to forced oscillations of non-linear systems governed by equations of the form u_ss?u_yy?εf(u,u,_y,u_yy…,s) = 0. The method is presented by application to the equation u_ss?u_yy?εu²_yu_yy= 0 which governs the vibrations of a soil layer that is free on the top surface and is forced harmonically at the bedrock. It is shown that unlike the ODE case (Duffing equation), the PDE requires an infinite number of periodicity conditions to correctly characterize the resonant region and these conditions lead to an infinite number of branches in the dispersion spectrum. Calculations indicate that these branches tend to an envelope curve. The uniform approach presented by Millmann and Keller is discussed in order to determine in what sense it can be viewed as an effective approximation for the fundamental mode.     

Asymptotic crack-tip fields for dynamic crack growth in compressive power-law hardening materials

The effect of material compressibility on the stress and strain fields for a mode-I crack propagating steadily in a power-law hardening material is investigated under plane strain conditions. The plastic deformation of materials is characterized by the J₂ flow theory within the framework of isotropic hardening and infinitesimal displacement gradient. The asymptotic solutions developed by the present authors [Zhu, X.K., Hwang K.C., 2002. Dynamic crack-tip field for tensile cracks propagating in power-law hardening materials. International Journal of Fracture 115, 323–342] for incompressible hardening materials are extended in this work to the compressible hardening materials. The results show that all stresses, strains, and particle velocities in the asymptotic fields are fully continuous and bounded without elastic unloading near the dynamic crack tip. The stress field contains two free parameters σ_eq0 and s₃₃₀ that cannot be determined in the asymptotic analysis, and can be determined from the full-field solutions. For the given values of σ_eq0 and s₃₃₀, all field quantities around the crack tip are determined through numerical integration, and then the effects of the hardening exponent n, the Poisson ratio ν, and the crack growth speed M on the asymptotic fields are studied. The approximate behaviors of the proposed solutions are discussed in the limit of ν → 0.5 or n → ∞.     

Hydrodynamic characteristics of a rectangular gas-driven inverse liquid-solid fluidized bed

The hydrodynamic characteristics of a rectangular gas-driven inverse liquid-solid fluidized bed (GDFB) using particles of different diameters and densities were investigated in detail. Rising gas bubbles cause a liquid upflow in the riser portion, enabling a liquid downflow that causes an inverse fluidization in the downer portion. Four flow regimes (fixed bed regime, initial fluidization regime, complete fluidization regime, and circulating fluidization regime) and three transition gas velocities (initial fluidization gas velocity, minimum fluidization gas velocity, and circulating fluidization gas velocity) were identified via visual observation and by monitoring the variations in the pressure drop. The axial local bed voidage (ε) of the downer first decreases and then increases with the increase of the gas velocity. Both the liquid circulation velocity and the average particle velocity inside the downer increase with the increase of the gas velocity in the riser, but decrease with the particle loading. An empirical formula was proposed to successfully predict the Richardson-Zaki index “n”, and the predicted ε obtained from this formula has a ±5% relative error when compared with the experimental ε.     

Molecular modeling of cracks at interfaces in nanoceramic composites

Toughness in Ceramic Matrix Composites (CMCs) is achieved if crack deflection can occur at the fiber/matrix interface, preventing crack penetration into the fiber and enabling energy-dissipating fiber pullout. To investigate toughening in nanoscale CMCs, direct atomistic models are used to study how matrix cracks behave as a function of the degree of interfacial bonding/sliding, as controlled by the density of C interstitial atoms, at the interface between carbon nanotubes (CNTs) and a diamond matrix. Under all interface conditions studied, incident matrix cracks do not penetrate into the nanotube. Under increased loading, weaker interfaces fail in shear while stronger interfaces do not fail and, instead, the CNT fails once the stress on the CNT reaches its tensile strength. An analytic shear lag model captures all of the micromechanical details as a function of loading and material parameters. Interface deflection versus fiber penetration is found to depend on the relative bond strengths of the interface and the CNT, with CNT failure occurring well below the prediction of the toughness-based continuum He–Hutchinson model. The shear lag model, in contrast, predicts the CNT failure point and shows that the nanoscale embrittlement transition occurs at an interface shear strength scaling as _τs~_εf,CNTσCNT${\tau }_s~{\epsilon }_{f,CNT}{\sigma }_{CNT}$ rather than _τs~_σCNT${\tau }_s~{\sigma }_{CNT}$ typically prevailing for micron scale composites, where _εf,CNT${\epsilon }_{f,CNT}$ and _σCNT${\sigma }_{CNT}$ are the CNT failure strain and stress, respectively. Interface bonding also lowers the effective fracture strength in SWCNTs, due to formation of defects, but does not play a role in DWCNTs having interwall coupling, which are weaker than SWCNTs but less prone to damage in the outerwall.