Computational fluid dynamics simulation of hydrodynamics in the riser of an external loop airlift reactor

Local hydrodynamics in the riser of an external loop airlift reactor (EL-ALR) are identified and the performances of three drag models are evaluated in computational fluid dynamics simulation. The simulation results show that the Schiller–Naumann drag model underestimated the local gas holdup at lower superficial gas velocity whereas the Tomiyama drag model overestimated that at higher superficial gas velocity. By contrast, the dual-bubble-size (DBS)-local drag model gave more reasonable radial and axial distributions of gas holdup in all cases. The reason is that the DBS-local drag model gave correct values of the lumped parameter, i.e., the ratio of the drag coefficient to bubble diameter, for varying operating conditions and radial positions. This ratio is reasonably expected to decrease with increasing superficial gas velocity and be smaller in the center and larger near the wall. Only the DBS-local drag model correctly reproduced these trends. The radial profiles of the axial velocity of the liquid and gas predicted by the DBS-local model also agreed well with experimental data.     

Effect of a pulsating flow on hydrodynamics and fluid catalytic cracking chemical reaction in a circulating fluidized bed riser

The three-dimensional computational fluid dynamics (3D-CFD) of a pulsating flow applied to the fluid catalytic cracking (FCC) reaction was investigated in the riser of a circulating fluidized bed reactor. The kinetic parameters of the FCC and coke burning reactions for predicting the reactant conversion and product yield percentages were applied. To increase the reactant conversion level and product yield, the effect of the pulsating flow operating parameters was considered using a 2^k statistical experimental design with four factors (amplitude, frequency, types of the waveform, and amplitude ratio). The 3D-CFD simulation was successfully validated from the experimental literature data. The frequency and type of the waveform were found to be the significant operating parameters. The expression of the fitted regression model and response surface contour were derived and revealed that the pulsating flow provides a higher reactant conversion level and product yield percentages compared to a non-pulsating or steady flow.     

Computational investigation of hydrodynamics,coal combustion and NOx emissions in a tangentially fired pulverized coal boiler at various loads

This work presents a computational investigation of hydrodynamics, coal combustion and NOx emissions in a tangentially fired pulverized coal boiler at different loads (630, 440 and 300 MW; relative loads of 100%, 70% and 48%) to clarify the effect of load change on the furnace processes. A computational fluids dynamics model was established; the flow field, temperature profile, species concentration and NOx emissions were predicted numerically; and the influence of burner tilt angles was evaluated. Simulation results indicate that a decrease in boiler load decreases the gas velocity, attenuates the airflow rotations, and increases the tangent circle size. The high-temperature zone and flame moved toward the side walls. Such behaviors impair air–fuel mixing, heat transfer and steady combustion in the furnace. In terms of species concentrations, a decrease in boiler load increased the O₂ content, decreased the CO content, and decreased the char burnout rates only slightly. A change in boiler load from 630 to 440 and 300 MW increased the NOx emissions from 202 to 234 and 247 mg/m³, respectively. Burner tilt angles are important in coal combustion and NOx emissions. A burner angle of –15° favors heat transfer and low NOx emissions (<185 mg/m³) for the current tangentially fired boiler.     

Computational and experimental study of annular photo-reactor hydrodynamics

The performance of ultraviolet (UV) reactors used for water treatment is greatly influenced by the reactor hydrodynamics due to the non-homogeneity of the radiation field. Reliable modeling of the reactor flow structures is therefore crucial for the design process. In this study, the turbulent flow through two characteristic annular UV-reactor configurations, with inlets concentric (L-shape) and normal (U-shape) to the reactor axis, was investigated through computational fluid dynamic (CFD). The modeling results were evaluated with the velocity profiles from particle image velocimetry (PIV) experiments. The influence of mesh structure and density, as well as three turbulence models: Standard κ–, Realizable κ–, and Reynolds stress model (RSM), on the simulation results were evaluated. Mesh-independent solutions were achieved at mean cell volumes of 5 × 10⁻⁹ m³. The Realizable κ– displayed the best overall match to the experimental PIV measurements. In general, the CFD models showed a close agreement with the experimental data for the majority of the reactor domain and captured the influences of reactor configuration and internal reactor structures on the flow distribution. The validated CFD hydrodynamic models could be integrated with kinetic and radiation distribution models for UV-reactor performance simulation.     

Effect of baffles on performance of fluid catalytic cracking riser

Increasing demand of automobile fuel and a need to process heavier crude oil makes it imperative to find improvements to the design of existing fluid catalytic cracking (FCC) units. Several modifications to the design of the riser section of FCC units have been suggested in previous studies including: improved feed nozzle designs, multiple nozzle configurations, internal baffles, and novel two-stage-riser systems. In this study, we investigate the effects of baffles on the performance of FCC risers using computational fluid dynamics simulations. In this study, predictions from a basis model (without baffles) are compared with those from four different configurations including: (i) 5-cm baffles at 5-m spacing, (ii) 7.5-cm baffles at 5-m spacing, (iii) 10-cm baffles with 5-m spacing, (iv) 10-cm baffles at 2.5-m spacing, and (v) 10-cm baffles at 1-m spacing. The baffles force the catalyst away from walls toward the center of the riser, enhancing the radial dispersion of the catalyst and the heat transfer inside the riser. The use of longer baffles and smaller spacings further increases the dispersion, yielding more homogeneous radial profiles. The changes in the radial dispersion result in variations in the conversion, yields, and pressure drops. The baffles increase conversion of vacuum gas oil (VGO) and the yield of gasoline. However, the simulations showed that longer baffles and a larger number of baffles did not always give a higher yield or higher conversion. Among the simulated configurations, the 5-cm baffles at 5-m spacing gave the highest conversion of VGO, whereas the 10-cm baffles at 1-m spacing resulted in the highest yield of the gasoline. Thus, rational optimization of baffle configurations is required to achieve optimal performance.     

Investigation on gas-solids hydrodynamics and matching enhancement in feedstock injection zone of double-level nozzle riser

With the development of current energy economy, it is necessary to improve the product distribution of fluid catalytic cracking process, which is achieved by a riser reactor with double-level of nozzles. The new riser is constructed by adding a level of secondary nozzle 0.5 m below the main nozzle of traditional riser. This paper investigates the gas-solids flow and oil-catalyst matching feature based on the optical fiber and tracer technologies. According to the distribution of solids holdup, particle velocity and dimensionless jet concentration, the feedstock injection zone can be divided into the upstream flow control region, the main flow control region, and the secondary flow control region in the radial direction. The size of the regions is changed by the jet gas velocity and axial height. There is a poor match of secondary nozzle jet to particles below the main nozzle. The jet gas from secondary nozzles can improve the matching effect of oil-catalyst near the wall and reduce the probability of coking above the main nozzle.     

Experimental investigation on vortex-induced vibrations of a hang-off evacuated drilling riser

Nonlinear Dynamics - A model experiment has been conducted for a hang-off evacuation riser considering the lower marine riser package (LMRP) based on the strain gauge test technology. The finite...     

Numerical simulation and experimental validation of gas-solid flow in the riser of a dense fluidized bed reactor

-Ladd models. As a result, the proposed drag force model can be used as an efficient approach for the dense gas-solid two-phase flow.     

CFD modeling of hydrodynamics and mass transfer of Rhodamine B in annular reactor

The hydrodynamics and mass transfer are the two crucial issues in annular reactors. An accurate prediction of these issues is required for design, optimization and scale-up applications. The present study deals with the modeling and simulation of flow through an annular reactor at different hydrodynamic conditions using computational fluid dynamics (CFD) to investigate the hydrodynamics and mass transfer. CFD modeling was utilized to predict velocity distribution, average velocity and average mass transfer coefficient in the annular geometry. The results of CFD simulations were compared with the mathematically derived equations and already developed correlations for validation purposes. CFD modeling was found suitable for predicting hydrodynamics and mass transfer for annular geometry under laminar flow conditions. It was observed that CFD also provides local values of the parameters of interest in addition to the average values for the simulated geometry.     

Experimental investigation of interparticle collision in the upper dilute zone of a cold CFB riser

Interparticle collision plays an important role in the mechanics of gas–solid two-phase flows. The paper presents direct measurements of collision rate as well as collision properties of spherical glass beads with sizes of 500 ± 50 μm in the upper dilute zone of a cold pilot-scale CFB riser, by using a high-speed imaging system. The recording rate of the high-speed digital camera is as high as 5000 fps with a resolution of 640 × 480. A large number of particle movement images at a height of 3.54 m above the distributor plates were taken. Manual inspection and automatic methods based on digital image processing algorithms were carried out to analyze particle image sequences. The experimental results show that the measured particle collision rate is proportional both to the particles’ average relative translational velocity and the square of the particle number density, which coincides with the collision theory derived according to the analogy of kinetic theory of gases. But the theoretical model is found to overestimate the real collision rates, and a coefficient a of 0.33 may be used to correct this discrepancy. The possible reasons for this discrepancy are also discussed. The measurement results of collision properties based on more than 50 particle collision events agrees well with Walton’s hard-sphere collision model. The three collision parameters, i.e., the average coefficient of friction μ, the normal and tangential coefficients of restitution e and β₀, for the glass beads used are measured to be 0.175 ± 0.005, 0.96 ± 0.02, and 0.43 ± 0.09, respectively.     

Numerical simulation and experimental validation of gas–solid flow in the riser of a dense fluidized bed reactor

Gas–solid flow in the riser of a dense fluidized bed using Geldart B particles (sand), at high gas velocity (7.6–15.5 m/s) and with comparatively high solid flux (140–333.8 kg/m² s), was investigated experimentally and simulated by computational fluid dynamics (CFD), both two- and three-dimensional and using the Gidaspow, O’Brien-Syamlal, Koch-Hill-Ladd and EMMS drag models. The results predicted by EMMS drag model showed the best agreement with experimental results. Calculated axial solids hold-up profiles, in particular, are well consistent with experimental data. The flow structure in the riser was well represented by the CFD results, which also indicated the cause of cluster formation. Complex hydrodynamical behaviors of particle cluster were observed. The relative motion between gas and solid phases and axial heterogeneity in the three subzones of the riser were also investigated, and were found to be consistent with predicted flow structure. The model could well depict the difference between the two exit configurations used, viz., semi-bend smooth exit and T-shaped abrupt exit. The numerical results indicate that the proposed EMMS method gives better agreement with the experimental results as compared with the Gidaspow, O’Brien-Syamlal, Koch-Hill-Ladd models. As a result, the proposed drag force model can be used as an efficient approach for the dense gas–solid two-phase flow.     

Computational and experimental investigation of flow over a transient pitching hydrofoil

The present study is developed within the framework of marine structure design operating in transient regimes. It deals with an experimental and numerical investigation of the time–space distribution of the wall-pressure field on a NACA66 hydrofoil undergoing a transient up-and-down pitching motion from 0° to 15° at four pitching velocities and a Reynolds number Re = 0.75 × 10⁶. The experimental investigation is performed using an array of wall-pressure transducers located on the suction side and by means of time–frequency analysis and Empirical Modal Decomposition method. The numerical study is conducted for the same flow conditions. It is based on a 2D RANS code including mesh reconstruction and an ALE formulation in order to take into account the foil rotation and the tunnel walls. Due to the moderate Reynolds number, a laminar to turbulent transition model was also activated. For the operating flow conditions of the study, experimental and numerical flow analysis revealed that the flow experiences complex boundary layer events as leading-edge laminar separation bubble, laminar to turbulent transition, trailing-edge separation and flow detachment at stall. Although the flow is relatively complex, the calculated wall pressure shows a quite good agreement with the experiment provided that the mesh resolution and the temporal discretization are carefully selected depending on the pitching velocity. It is particularly shown that the general trend of the wall pressure (low frequency) is rather well predicted for the four pitching velocities with for instance a net inflection of the wall pressure when transition occurs. The inflection zone is reduced as the pitching velocity increases and tends to disappear for the highest pitching velocity. Conversely, high frequency wall-pressure fluctuations observed experimentally are not captured by the RANS model. Based on the good agreement with experiment, the model is then used to investigate the effects of the pitching velocity on boundary layer events and on hydrodynamic loadings. It is shown that increasing the pitching velocity tends to delay the laminar-to-turbulence transition and even to suppress it for the highest pitching velocity during the pitch-up motion. It induces also an increase of the stall angle (compared to quasi-static one) and an increase of the hysteresis effect during pitch-down motion resulting to a significant increase of the hydrodynamic loading.     

Modeling of the hydrodynamics of a facility for removing mechanical impurities from oil

The separation of a viscous fluid flow by means of a device operating on the centrifugal principle was modeled. Parameters of the device with a specified separation size were obtained in a series of numerical experiments. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 4, pp. 108–112, July–August, 2008.     

Experimental investigation of oil transport in a compound vortex

Experiments were performed to study oil transport in the compound vortex produced by a uniformly rotating disk in a cylindrical container. It was found that on the liquid surface in the neighborhood of the axis of rotation, light oil formed a spot with spiral branches or broke up into separate elongated drops, and the spiral branches stretched out in the direction opposite to the direction of rotation of the liquid. In the bulk of the liquid in the vicinity of the axis of its rotation, oil formed a compact body. The geometric parameters of the structural components of the flows were determined.     

Numerical comparison of two modes of gas-solid riser operation: Fluid catalytic cracking vs CFB combustor

Two modes of gas-solid riser operation, i.e., fluid catalytic cracking (FCC) and circulating fluidized bed combustor (CFBC), have been recognized in literature; particularly in the understanding of choking phenomena. This work compares these two modes of operation through computational fluid dynamics (CFD) simulation. In CFD simulations, the different operations are represented by fixing appropriate boundary conditions: solids flux or solids inventory. It is found that the FCC and CFBC modes generally have the same dependence of solids flux on the mean solids volume fraction or solids inventory. However, during the choking transition, the FCC mode of operation needs more time to reach a steady state; thus the FCC system may have insufficient time to respond to valve adjustments or flow state change, leading to the choking. The difference between FCC and CFBC systems is more pronounced for the systems with longer risers. A more detailed investigation of these two modes of riser operation may require a three-dimensional full loop simulation with dynamic valve adjustment.     

Effect of static bed height in the upper fluidized bed on flow behavior in the lower riser section of a coupled reactor

To study olefin reduction by using an auxiliary reactor for FCC naphtha upgrading, a large-scale cold model of a riser-bed coupled to an upper fluidized bed was established. The effect of static bed height in the upper fluidized bed on particle flow behavior in the lower riser was investigated experimentally. A restriction index of solids holdup was used to evaluate quantitatively the restrictive effect of the upper fluidized bed. Experimental results show that, under the restrictive effect of the upper fluidized bed, the riser could be divided into three regions in the longitudinal direction: accelerating, fully developed and restriction. The axial distribution of solids holdup in the riser is characterized by large solids holdup in the top and bottom sections and small solids holdup in the middle section. Overall solids holdup increased with increasing static bed height in the upper fluidized bed, while particle velocity decreased. Such restrictive effect of the upper fluidized bed could extend from the middle and top sections to the whole riser volume when riser outlet resistance is increased, which increases with increasing static bed height in the upper fluidized bed. The upper bed exerts the strongest restriction on the area close to the riser outlet.     

Effect of static bed height in the upper fluidized bed on flow behavior in the lower riser section of a coupled reactor

To study olefin reduction by using an auxiliary reactor for FCC naphtha upgrading, a large-scale cold model of a riser-bed coupled to an upper fluidized bed was established. The effect of static bed height in the upper fluidized bed on narticle flow behavior in the lower riser was investigated experimentally. A restriction index of solids holdup was used to evaluate quantitatively the restrictive effect of the upper fluidized bed. Experimental results show that, under the restrictive effect of the upper fluidized bed, the riser could be divided into three regions in the longitudinal direction： accelerating, fully developed and restriction. The axial distribution of solids holdup in the riser is characterized by large solids holdup in the top and bottom sections and small solids holdup in the middle section. Overall solids holdup increased with increasing static bed height in the upper fluidized bed, while particle velocity decreased. Such restrictive effect of the upper fluidized bed could extend from the middle and top sections to the whole riser volume when riser outlet resistance is increased, which increases with increasing static bed height in the upper fluidized bed. The upper bed exerts the strongest restriction on the area close to the riser outlet.     

Computational aeroelastic investigation of a transonic limit-cycle-oscillation experiment at a transport aircraft wing model

Aeroelastic measurements of a three-dimensional wing model, the so-called Aerostabil wing, were conducted in the Transonic Windtunnel Göttingen. This clean, backward-swept wing allowed the experimental investigation of limit cycle oscillations in a certain transonic parameter range. In this paper, a detailed insight into the observed physical phenomena, especially the measured limit cycle oscillations, is presented by means of CFD–CSM coupled simulations. These simulations on the basis of a detailed structural finite element model reveal the specific properties of the Aerostabil wing and furthermore allow investigating the unstable behavior of this windtunnel model for transonic flow settings. The aerodynamic characteristics include a two-shock system and large flow separation areas, further increasing the complexity of the aeroelastic problem. A structural single degree-of-freedom system is used for the prediction of the experimental stability range and the limit cycle oscillation investigations. Due to the good agreement of simulation and experiment the limit cycle oscillations can be explained by means of nonlinear aerodynamic effects.