鼓泡流化床中流动特性的多尺度数值模拟

鼓泡流化床因其较高的传热特性以及较好的相间接触已经被广泛应用于工业生产中,而对鼓泡流态化气固流动特性的充分认知是鼓泡流化床设计的关键.在鼓泡流化床中,气泡相和乳化相的同时存在使得床中呈现非均匀流动结构,而这种非均匀结构给鼓泡流化床的数值模拟造成了很大的误差.基于此,以气泡作为介尺度结构,建立了多尺度曳力消耗能量最小的稳定性条件,构建了适用于鼓泡流化床的多尺度气固相间曳力模型.结合双流体模型,对A类和B类颗粒的鼓泡流化床中气固流动特性进行了模拟研究,分析了气泡速度、气泡直径等参数的变化规律.研究表明,与传统的曳力模型相比,考虑气泡影响的多尺度气固相间曳力模型给出的曳力系数与颗粒浓度的关系是一条分布带,建立了控制体内曳力系数与局部结构参数之间的关系.通过模拟得到的颗粒浓度和速度与实验的比较可以发现,考虑气泡影响的多尺度曳力模型可以更好地再现实验结果.通过A类和B类颗粒的鼓泡床模拟研究发现,A类颗粒的鼓泡床模拟受多尺度曳力模型的影响更为显著.      

Local heat transfer coefficients on the circumference of a tube in a gas fluidized bed

A heated horizontal heat transfer tube was installed 14.8 cm above the distributor plate in a square fluid bed measuring 30.5 × 30.5 cm. Four different Geldart B sized particle beds were used (sand of two different distributions, an abrasive and glass beads) and the bed was fluidized with cold air. The tube was instrumented with surface thermocouples around half of the tube circumference and with differential pressure ports that can be used to infer bubble presence. Numerical execution of the transient conduction equation for the tube allowed the local time-varying heat transfer coefficient to be extracted. Data confirm the presence of the stagnant zone on top of the tube associated with low superficial velocities. Auto-correlation of thermocouple data revealed bubble frequencies and the cross-correlation of thermal and pressure events confirmed the relationship between the bubbles and the heat transfer events. In keeping with the notion of a “Packet renewal” heat transfer model, the average heat transfer coefficient was found to vary in sympathy with the root-mean square amplitude of the transient heat transfer coefficient.     

DEM simulation of polydisperse systems of particles in a fluidized bed

Numerical simulations based on three-dimensional discrete element model (DEM) are conducted for mono-disperse, binary and ternary systems of particles in a fluidized bed. Fluid drag force acting on each particle depending on its size and relative velocity is assigned. The drag coefficient corresponding to Ergun’s correlation is applied to the system of fluidized bed with particle size ratios of 1:1 for the mono-disperse system, 1:1.2, 1:1.4 and 1:2 for the binary system and 1:1.33:2 for the ternary system b...     

Effect of temperature gradient within a solid particle on the rotation and oscillation modes in solid-dispersed two-phase flows

Liquid–solid two-phase flow with heat transfer is simulated, and the effect of temperature gradient within a solid particle on the particle behaviour and heat transfer is studied. The interaction between fluid and particles is considered with our original immersed solid approach on a rectangular grid system. The local heat flux at the fluid–solid interface is described with an anisotropic heat conductivity matrix, and the governing equation of temperature is time-updated with an implicit treatment for the diffusion term. The method is applied to a 2-D natural convection flow of a relatively low Rayleigh number including multiple particles. Heat transfer and particle behaviours are studied for different solid heat conductivities (ratio to the fluid conductivity ranging between 10⁻³ and 10³) and solid volume fractions. Under a condition of relatively low heat conductivity ratio, the particles show a simple circulating flow. By increasing the heat conductivity ratio, a transition of the particulate flow is observed to oscillation mode around the domain centre due to the buoyancy force as a restitution force. The oscillation period is found to vary with the heat conductivity ratio, and it is related to the time scales for the heat transfer via fluid and solid.     

Gas–liquid mass transfer in the gas–liquid–solid mini fluidized beds

The gas–liquid–solid mini fluidized bed (GLSMFB) combines the advantages of fluidized bed and micro-reactor, and meets the requirements for safety and efficiency of green development of process industry. However, there are few studies on its flow performance and no studies on its mass and heat transfer performance. In this paper, the characteristics of gas–liquid mass transfer in a GLSMFB were studied in order to provide basic guidance for the study of GLSMFB reaction performance and application. Using CO₂ absorption by NaOH as the model process, the gas–liquid mass transfer performance of GLSMFB was investigated. The results show that the liquid volumetric mass transfer coefficient and the gas–liquid interfacial area both increase with the increase of the superficial gas velocity within the experimental parameter range under the same given superficial liquid velocity. At the same ratio of superficial gas to liquid velocity, the liquid volumetric mass transfer coefficient increases with the increase of the superficial liquid velocity. Fluidized solid particles strengthen the liquid mass transfer process, and the liquid volumetric mass transfer coefficient is about 13% higher than that of gas–liquid mini bubble column.     

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.     

Fluidization and bubbling behavior of potash particles in a deep fluidized bed

Most existing models for predicting bubble size and bubble frequency have been developed for freely bubbling fluidized beds. Accurate prediction of bubbling behavior in deep fluidized beds, however, has been a challenge due to the higher degree of bubble coalescence and break up, high probability of the slugging regime, partial fluidization, and chaotic behavior in the bubbling regime. In this work, the bubbling and fluidization behavior of potash particles was investigated in a deep fluidized bed employing a twin-plane electrical capacitance tomography (ECT) system. Solid volume fraction, average bubble velocity, average bubble diameter, and bubble frequency in both bubbling and slugging regimes were measured at two different bed height ratios (H/D = 3.5 and H/D = 3.78). This work is the first to illustrate a sequential view of bubbles at different superficial gas velocities in a fluidized bed. The results show that both the bubble diameter and rising velocity increased with increasing the superficial gas velocity for the two bed heights, with larger values observed in the deeper bed compared to the shallower one. Predicted values for bubble diameter, bubble rise velocity and bubble frequency from different models are compared with the experimental data obtained from the ECT system in this work. Good agreement has been achieved between the values predicted by the previous models and the experimental data for the bubble diameter and bubble rise velocity with an average absolute deviation of 16% and 15% for the bed height of 49 cm and 13% and 8% for the bed height of 53 cm, respectively.     

Bubble behavior in a three-phase fluidized bed

Bubble size measurements by off-axis holography are reported for the first time in three-phase fluidized beds. The solid material was glass cylinders, with the liquid in one case having matching refractive index. Two different columns, 76 and 152 mm dia, were employed. An interpolative correlation was developed for the solid volume fraction of the bed as a function of liquid flux. A generalized dimensionless model was developed to predict the equilibrium bubble size in both fixed and fluidized beds.     

Experimental study and modeling of particle drying in a continuously-operated horizontal fluidized bed

Multistage fluidized beds are frequently used for product drying in industry. One advantage of these fluidized beds is that they can achieve a high throughput, when operated continuously. In this study, γ-Al₂O₃ particles were dried in a pilot-scale horizontal fluidized bed, without considering any comminution effects. For each experiment, the particle moisture content distribution and residence time distribution were determined. To take into account particle back mixing in our experiments, a one-dimensional population balance model that considers particle residence time was introduced into a fluidized bed-drying model. Experimental particle residence time distributions were reproduced using a tank-in-series model. Subsequently, the moisture content distribution was implemented, as a second dimension to the population balance in this model. These two-dimensional simulations were able to describe the experimental data, especially the spread in the residual particle moisture distribution, much more accurately than one-dimensional simulations. Using this novel two-dimensional model, the effects of different operating parameters (process gas temperature, solid feed rate, superficial air velocity) on the particle moisture content distribution were systematically studied.     

Transient behavior of heat removal from a cylindrical heat storage vessel packed with spherical porous particles

The transient heat transfer behavior in the case of heat removal from a cylindrical heat storage vessel packed with spherical particles was investigated experimentally for various factors (flow rate, diameter of spherical particles packed, temperature difference between flowing cold air and spherical particles accumulating heat, and physical properties of spherical particles). The experiments were covered in ranges of Reynolds number based on the mean diameter of spherical particles packed Re_d = 10.3–2200, porosity?=0.310 to 0.475, ratio of spherical particle diameter to cylinder diameterd/D = 0.0075–0.177 and ratio of length of the cylinder to cylinder diameterL/D=2.5–10. It was found that especially the flow rate and the dimension of spherical particles played an important role in estimating the transient local heat transfer characteristics near the wall of the cylindrical vessel in the present heat storage system. As flow rate and diameter of spherical particles were increased under a given diameter of the cylinder heat storage vessel, the mean heat transfer coefficient between the flow cold air and the hot spherical particles increased and the time period to finish removing heat from the vessel reduced. In addition, the useful experimental correlation equations of mean heat transfer coefficient between both phases and the time period to finish removing heat from the vessel were derived with the functional relationship of Nusselt numberNu _d=f [modified Prandtl numberPr ^* (d/D), Re_d) and Fourier numberFo = f(d/D, L/D, Pr^*, Re_d).     

Influence of vertical internals on a bubbling fluidized bed characterized by X-ray tomography

An ultra-fast X-ray tomographic scanner is applied to study the hydrodynamics in a bubbling fluidized bed with and without vertical internals (e.g., heat exchanger tubes). The objective of this study is to understand the influence of vertical internals on hydrodynamic properties such as bubble volume, size and velocity and to provide measurement data for the design and scale-up of catalytic bubbling fluidized bed reactors with vertical internals. With these new measurements, correlations of bubble properties can be developed to reliably scale-up bubbling fluidized beds with vertical internals. For the investigated reactor with Geldart A/B particles, no relation between bubble size and velocity was observed for individual bubbles, i.e.; smaller bubbles tend to rise with higher velocities. A significant reduction in bubble size and sharpening of the bubble size distribution was generally obtained for a bed with vertical internals.     

The mathematical modelling of thermal radiation in high temperature fluidized bed

The aim of this study is composed of two parts. One of them is to calculate the radiation heat flux and the other is to determine the overall heat transfer coefficient for the gas-fluidized bed. The radiative heat transfer model is developed for predicting the total heat transfer coefficients between submerged surfaces and fluidized beds for several working temperatures. The role of radiation heat transfer in the overall heat transfer process at an immersed surface in a gas-fluidized bed at high temperatures is investigated. Analytical results are compared with the previously done experiments and a good agreement between the two, is obtained.

---

Bestimmung der Wärmeübertragungs-Koeffizienten in Gas-Wirbelschichten

Zusammenfassung Diese Untersuchung besteht aus folgenden zwei Teilen: 1. Kalkulation des Radiationswärmeübergangs in Gas-Wirbelschichten. 2. Bestimmung des Wärmeübergangs-Koeffizienten in Gas-Wirbelschichten. Dieses Radiationswärmeübergangsmodell wurde entwickelt, um die Wärmeübertragungs-Koeffizienten zwischen der eingetauchten Oberfläche und der Wirbelschicht bei verschiedener Wärme schätzungsweise zu bestimmen. Es wurde das Verhältnis der Radiationswärmeübertragung in Gas-Wirbelschichten zum totalen Wärmeübergang untersucht. Die Meßwerte wurden mit theoretischen Resultaten verglichen.

---

Nomenclature c (x) specific heat capacity of packet [J/kg K] - c _p specific heat capacity of particle [J/kg K] - c _pg specific heat capacity of gas [J/kg K] - d _p average diameter of the bed particles [m] - f ₀ the fraction of time that a unit surface exposed to the bubble phase - 1–f ₀ the fraction of time that a unit surface exposed to the packet phase - g acceleration due to gravity [m/s²] - h _b heat transfer coefficient for the surface in contact with bubble [W/m² K] - h _bc conduction heat transfer coefficient for the surface/bubble [W/m²K] - h _br radiation heat transfer coefficient for the surface/bubble [W/m²K] - h _p heat transfer coefficient for the surface in contact with packet [W/m²K] - h _pc conduction heat transfer coefficient for the surface/packet [W/m² K] - h _pr radiation heat transfer coefficient for the surface/packet [W/m² K] - h _T total heat transfer coefficient between bed and surface [W/m² K] - k ⁰ thermal conductivity of the emulsion phase for fixed bed [W/m K] - k(x) thermal conductivity of packet [W/m K] - k _e the logarithmic mean of conductivity for first layer in packet [W/m K] - k _g the logarithmic mean of conductivity for the first layer in packet [W/m K] - K extinction coefficient [1/m] - m mass [kg] - n number of layers - p air pressure [pa] - q _pc mean local conduction heat transfer for packet [kW/m²] - q _pr mean local radiation heat transfer for packet [kW/m²] - Q _p average heat flux during packet contact with surface [kW/m²] - Q _b average heat flux during bubble contact with surface [kW/m²] - R gas constant [287.04 J/kg K] - t time [s] - t _g residence time for gas bubble [s] - t _k residence time for packet [s] - T temperature [K] - T _b bed temperature [K] - T _W surface temperature [K] - V _mf minimum fluidization velocity [m/s] - v _t terminal velocity [m/s] - x distance [m] Greek symbols t time increment - x thickness of the layer - emissivity - thermal diffusivity [m²/s] - (x) voidage of fluidized bed - _mf void ratio of the bed at minimum fluidization - ₀ voidage of fixed bed - _g dynamic viscosity of gas [kg/m s] - _g kinematic viscosity of gas [m²/s] - (x) density of packet [kg/m³] - _p density of particles [kg/m³] - _g density of gas [kg/m³] - Stefan-Boltzmann constant [5.66·10^–8 W/m²K⁴] - geometric shape factor for particles Dimensionless numbers Ar Archimedes numberAr=g d _p ³ ( _p – _g ) _g / _g ² - Nu Nusselt numberNu=h·d/k - Re Reynolds numberRe=d _p ·V _mf / _g - Pr Prandtl numberPr=C _{pg g} /k _g      

Local heat transfer coefficients around a horizontal heated tube immersed in a gas fluidized bed

The heat transfer characteristics around a single horizontal heated tube immersed in air fluidized bed was investigated, to clarify the mechanism of heat transfer in a fluidized bed heat exchanger. The local heat transfer coefficient around the tube was measured at various fluidization velocities and five different solid particles. The experimental values of the local heat transfer coefficient at the minimum fluidization velocity condition were correlated with the particle size in two empirical equations. The predicted results were in good agreement with the experiment data.     

Analysis of temperature field in brake disc for fade assessment

Simulation of the temperature distribution in a brake disc of a passenger car during fade cycles (repeated brake application) is carried out using Finite Element Method. Verification of the predicted values is done using a sophisticated inertia dynamometer. Effect of variable convective heat transfer coefficient, effect of mass of the mounting flange on the temperature distribution are also studied. From the investigation it is concluded that the present model can be used for the simulation of temperature distribution in rigid brake dises during fade test.

---

Berechnung des Temperaturfeldes in Bremsscheiben zur Fading-Abschätzung

Zusammenfassung Mit einer Finite-Elemente-Methode wird die Temperaturverteilung in einer Bremsscheibe eines Personenkraftwagens während der Fading-Zyklen (wiederholtes Abbremsen) ermittelt. Mittels eines hochentwickelten Schwungmassen-Bremsprüfstandes werden die berechneten Werte bestätigt. Der Einfluß der veränderlichen konvektiven Wärmeübergangskoeffizienten, sowie der Einfluß der Masse des Befestigungsflansches auf die Temperaturverteilung werden betrachtet. Aus den Untersuchungen wird der Schluß gezogen, daß das vorliegende Modell für die Simulierung der Temperaturverteilung in starren Bremsscheiben bei Fading-Tests verwendet werden kann.

---

Nomenclature c specific heat - h convective heat transfer coefficient - k _r thermal conductivity inr direction - k thermal conductivity in direction - Mu coefficient of friction of the brake lining - T temperature - t time - density - Q Heat input - r radius - § Heat partition ratio - A Surface area - M Vehicle Mass - U _i Initial velocity - U _f Final velocity - BR Braking ratio - V Volume of heat generation - K _z Thermal conductivity inZ direction     

Parametric studies of wall heat transfer in a fluidized bed circulating through a packed bed

In the analysis presented, the circulating fluidized bed is treated as equivalent to a gas-solid suspension flow. The heat transfer is modelled in terms of the effective increase in thermal conductivity of the gas. The results show that the increase in thermal conductivity of the continuum due to the molecular conductivity of the circulating particles has a much lower contribution than particle eddy conductivity. High mass ratio and low size of the circulating solids are beneficial. The results are a lower bound to the heat transfer increase possible. Additional contribution due to gas convection at high velocities can be simply modelled, by increasing the particle eddy conductivity suitably.

---

Parameter-Studien von Wärmeübertragung in einer Wirbelschicht die durch ein Festbett zirkuliert

Zusammenfassung In der vorliegenden Untersuchung wird die zirkulierende Wirbelschicht wie eine Gas-Festkörper-Suspensions-strömung behandelt. Die Wärmeübertragung wird in Abhängigkeit vom effektiven Anstieg der thermischen Leitfähigkeit des Gases dargestellt. Die Ergebnisse zeigen, daß der Anstieg der thermischen Leitfähigkeit des Kontinuums in Folge von molekularer Leitung der zirkulierenden Partikel weniger bedeutend ist als die Wirbelleitfähigkeit. Große Massenverhältnisse und kleine Größen der zirkulierenden Festkörper sind vorteilhaft. Die Ergebnisse zeigen die mögliche Verbesserung der Wärmeübertragungswerte. Diese können in Folge von Gaskonvektion bei hohen Geschwindigkeiten durch geeignetes Ansteigen der Wirbelleitfähigkeit verbessert werden.

---

Nomenclature B non-dimensional conductivity - C _p specific heat - d _p bed particle diameter - f solid volume fraction - Fo Fourier number - h heat transfer coefficient - k _B bulk conductivity - k _e equivalent conductivity - M mass ratio of the circulating bed - Nu Nusselt number - q heat flux - t _i initial temperature - t _o wall temperature - x distance from the wall Greek symbols _B bulk voidage of the circulating bed - _g density of gas - _s solid density - time     

Bubble wake solids content in three-phase fluidized beds

It is shown that existing equations for predicting the holdups of wakes behind bubbles in three-phase fluidized beds are not entirely satisfactory. A new model is then developed whereby the wake is treated as the sphere-completing volume of a spherical cap bubble, due allowance being made for hydrodynamic interactions between bubbles. The generalized wake equations of Bhatia & Epstein (1974) are applied to compute the ratio of solids holdup in the wakes to that in the remaining liquid of the bed. Using experimental data from the literature, a rational equation is then generated for predicting this ratio from measured variables, and a mechanism for wake solids entrainment is proposed which is consistent with this equation.     

A novel dual-material probe for in situ measurement of particle charge densities in gas–solid fluidized beds

Particle charge density is vitally important for monitoring electrostatic charges and understanding particle charging behavior in fluidized beds. In this paper, a dual-material probe was tested in a gas–solid fluidized bed for measuring the charge density of fluidized particles. The experiments were conducted in a two-dimensional fluidized bed with both single bubble injection and freely bubbling, at various particle charge densities and superficial gas velocities. Uniformly sized glass beads were used to eliminate complicating factors at this early stage of probe development. Peak currents, extracted from dynamic signals, were decoupled to determine charge densities of bed particles, which were found to be qualitatively and quantitatively consistent with charge densities directly measured by Faraday cup from the freely bubbling fluidized bed. The current signals were also decoupled to estimate bubble rise velocities, which were found to be in reasonable agreement with those obtained directly by analyzing video images.     

A novel dual-material probe for in situ measurement of particle charge densities in gas-solid fluidized beds

Particle charge density is vitally important for monitoring electrostatic charges and understanding particle charging behavior in fluidized beds.In this paper,a dual-material probe was tested in a gas-solid fluidized bed for measuring the charge density of fluidized particles.The experiments were conducted in a two-dimensional fluidized bed with both single bubble injection and freely bubbling,at various particle charge densities and superficial gas velocities.Uniformly sized glass beads were used to eliminate complicating factors at this early stage of probe development.Peak currents,extracted from dynamic signals,were decoupled to determine charge densities of bed particles,which were found to be qualitatively and quantitatively consistent with charge densities directly measured by Faraday cup from the freely bubbling fluidized bed.The current signals were also decoupled to estimate bubble rise velocities,which were found to be in reasonable agreement with those obtained directly by analyzing video images.