Transport and deposition of neutral particles in magnetohydrodynamic turbulent channel flows at low magnetic Reynolds numbers

The effect of Lorentz force on particle transport and deposition is studied by using direct numerical simulation of turbulent channel flow of electrically conducting fluids combined with discrete particle simulation of the trajectories of uncharged, spherical particles. The magnetohydrodynamic equations for fluid flows at low magnetic Reynolds numbers are adopted. The particle motion is determined by the drag, added mass, and pressure gradient forces. Results are obtained for flows with particle ensembles of various densities and diameters in the presence of streamwise, wall-normal or spanwise magnetic fields. It is found that the particle dispersion in the wall-normal and spanwise directions is decreased due to the changes of the underlying fluid turbulence by the Lorentz force, while it is increased in the streamwise direction. The particle accumulation in the near-wall region is diminished in the magnetohydrodynamic flows. In addition, the tendency of small inertia particles to concentrate preferentially in the low-speed streaks near the walls is strengthened with increasing Hartmann number. The particle transport by turbophoretic drift and turbulent diffusion is damped by the magnetic field and, consequently, particle deposition is reduced.     

Drag models for simulating gas–solid flow in the turbulent fluidization of FCC particles

This paper examines the suitability of various drag models for predicting the hydrodynamics of the turbulent fluidization of FCC particles on the Fluent V6.2 platform. The drag models included those of Syamlal–O’Brien, Gidaspow, modified Syamlal–O’Brien, and McKeen. Comparison between experimental data and simulated results showed that the Syamlal–O’Brien, Gidaspow, and modified Syamlal–O’Brien drag models highly overestimated gas–solid momentum exchange and could not predict the formation of dense phase in the fluidized bed, while the McKeen drag model could not capture the dilute characteristics due to underestimation of drag force. The standard Gidaspow drag model was then modified by adopting the effective particle cluster diameter to account for particle clusters, which was, however, proved inapplicable for FCC particle turbulent fluidization. A four-zone drag model (dense phase, sub-dense phase, sub-dilute phase and dilute phase) was finally proposed to calculate the gas–solid exchange coefficient in the turbulent fluidization of FCC particles, and was validated by satisfactory agreement between prediction and experiment.     

Migration of settling particles in a horizontal viscous flow through a vertical slot with porous walls

Particle migration in a horizontal flow of dilute suspension through a vertical slot with porous walls is studied using the two-continua approach. The lateral migration is induced by two opposite effects: an inertial lift force due to particle settling and directed toward the slot centre-line, and a drag due to leak-off entraining particles toward the walls. An expression for the inertial lift on a settling particle in a horizontal channel flow found recently is generalized to the case of a low leak-off velocity. The evolution of an initial uniform particle concentration profile is studied within the full Lagrangian approach. Four migration regimes are found differing by the direction of particle migration and numbers of equilibrium positions. Conditions of the regime change and a critical value of dimensionless leak-off velocity for particle deposition on the walls are obtained analytically. Suspension flows with zones where the particle concentration is zero or increases infinitely, are studied numerically.     

On a Paradox in the Motion of a Rigid Particle along a Wall in a Fluid

The results of an experimental investigation of spherical particles with different surface roughnesses rolling under their own weight down an inclined pipe wall in a Newtonian fluid at low Reynolds numbers, both with (friction should be taken into account) and without contact with the wall, are presented. It is shown that a fixed particle moves differently in different fluids with similar viscosities and densities. This fact, as well as the possibility of particle motion without contact with the wall, cannot be explained within the framework of the usual hydrodynamic theories. An example is the dependence of the particle motion on the static pressure.     

Investigation on development of free-surface structures in turbulent liquid lithium flow

The liquid lithium film thickness facing the Deuterium beam of the International Fusion Material Irradiation Facility (IFMIF) determines the neutron flux to be generated. Hence, apart from its thickness also its spatio-temporal behaviour plays a decisive role in the performance of the target. Two aspects contributing to the free surface shape are the evolution of the viscous wall boundary layer in the nozzle and the development of turbulence downstream the nozzle exit, which are analysed here numerically by means of a fluid dynamic Large Eddy Simulation (LES). The numerical method is validated by experiments conducted at Osaka University with respect to mean and turbulent flow quantities in a broad spectrum of mean flow velocities. Thereby, both a qualitative and a quantitative agreement have been attained identifying different flow regimes and, moreover, allowing for a more refined, realistic IFMIF target prediction performance.     

On the friction drag reduction mechanism of streamwise wall fluctuations

Understanding how to decrease the friction drag exerted by a fluid on a solid surface is becoming increasingly important to address key societal challenges, such as decreasing the carbon footprint of transport. Well-established techniques are not yet available for friction drag reduction. Direct numerical simulation results obtained by Józsa et al. (2019) previously indicated that a passive compliant wall can decrease friction drag by sustaining the drag reduction mechanism of an active control strategy. The proposed compliant wall is driven by wall shear stress fluctuations and responds with streamwise wall velocity fluctuations. The present study aims to clarify the underlying physical mechanism enabling the drag reduction of these active and passive control techniques. Analysis of turbulence statistics and flow fields reveals that both compliant wall and active control amplify streamwise velocity streaks in the viscous sublayer. By doing so, these control methods counteract dominant spanwise vorticity fluctuations in the near-wall region. The lowered vorticity fluctuations lead to an overall weakening of vortical structures which then mitigates momentum transfer and results in lower friction drag. These results might underpin the further development and practical implementation of these control strategies.     

双柱状颗粒在线性剪切流场中的动力学分析

为研究柱状颗粒在线性剪切流场中的运动状态和受力情况,本文以颗粒长径比为2,颗粒之间的初始距离ΔS_Py=4D为例,基于直接力浸入边界法数值模拟了双柱状颗粒在三维线性剪切流场中的运动过程。根据模拟结果分析了柱状颗粒周围流场参数分布,在考虑壁面对颗粒的影响和颗粒之间相互影响的条件下,研究了颗粒的受力和运动的变化,探索了流体曳力导致柱状颗粒迁移和转动的规律。研究结果表明,双柱状颗粒在线性剪切流场中易向速度大的流体区域运动;前后两颗粒运动状态和轨迹不同,颗粒之间距离较近时,曳力会产生较大的波动;只有当颗粒在壁面附近时,滞后颗粒才能追上领先颗粒,两颗粒发生牵引、翻滚和分离过程。     

On the origin of drag increase in varying-phase opposition control

This study investigates the physics underlying the drag increase in a low Reynolds number turbulent channel flow due to varying-phase opposition control by means of direct numerical simulation and modal analysis. The drag increase occurs for an extended region of the parameter space and we consider a controller with a positive phase shift in Fourier domain between sensor measurement and actuator response as a representative example for this regime. Analyses of instantaneous flow fields as well as spatial power spectra show that the structure of drag-increased flows is remarkably different from that of drag-reduced and canonical flows. In particular, the near-wall region is dominated by structures of short streamwise and large spanwise extent. Isolation of a representative control scale shows that these energetic structures can be characterized as spanwise rollers, which induce strong ejection and sweep motions and lead to drag increase. The presence of rollers, and therefore drag increase, in the full nonlinear system correlates well with the presence of an amplified eigenvalue in the eigenspectrum of the linearized Navier–Stokes operator. It is further shown that the scales responsible for drag increase at positive phase shifts are inactive at negative phase shifts and do not contribute to drag reduction. These scales can therefore be excluded from a controller aimed at drag reduction, which relaxes the spatial resolution requirements on the control hardware. The eigenspectrum may be used as a computationally cheap tool to identify such detrimental scales during an early design stage.     

Effect of coarse particles on the heat transfer in a particle-laden turbulent boundary layer

The temperature distribution in particle-laden turbulent flow, in a flume, was investigated both by DNS and experimentally. Simulations were performed at Re=171 and Pr=5.4 in order to study the interaction between the particle motion and flow turbulence. Two-way coupling was used to obtain various turbulence statistics, the grid resolution was sufficiently fine to resolve all essential turbulent scales. The effect of particle diameter on momentum, heat transfer and particle deposition was considered. The details of particle-turbulence interaction depend on the particle Stokes number and the particle Reynolds number.

The spatial structures of instantaneous flow and temperature fields were visualized. Low frequency small oscillations of deposited particles were observed. It was found that these small deviations from the initial position, caused strong changes in the instantaneous temperature field near the particle.

The experiments provided details of the temperature field on the heated wall close to the particle. In the front of the particle, a sharp increase in heat transfer coefficient was observed. The experimental results agree well with the computational predictions.     

On the liquid turbulence energy spectra in two-phase bubbly flow in a large diameter vertical pipe

The liquid turbulence structure of air–water bubbly flow in a 200 mm diameter vertical pipe was experimentally investigated. A dual optical probe was used to measure the bubble characteristics, while the liquid turbulence was measured using hot-film anemometry. Experiments were performed at two liquid superficial velocities of 0.2 and 0.68 m/s for gas superficial velocities in the range of 0–0.18 m/s, corresponding to an area averaged void fraction up to 13.6%. In general, there is an increase in the liquid turbulence energy when the bubbles are introduced into the liquid flow. The increase in the energy mainly occurs over a range of length scales that are on the order of the bubble diameter. A suppression of the turbulence was observed close to the wall at very low void fraction flows. Initially, the suppression occurs in the low wave number range and then extends to higher wave numbers as the suppression is increased.     

Measurement of the liquid film thickness in micro tube slug flow

Slug flow is one of the representative flow regimes of two-phase flow in micro tubes. It is well known that the thin liquid film formed between the tube wall and the vapor bubble plays an important role in micro tube heat transfer. In the present study, experiments are carried out to clarify the effects of parameters that affect the formation of the thin liquid film in micro tube two-phase flow. Laser focus displacement meter is used to measure the thickness of the thin liquid film. Air, ethanol, water and FC-40 are used as working fluids. Circular tubes with five different diameters, D = 0.3, 0.5, 0.7, 1.0 and 1.3 mm, are used. It is confirmed that the liquid film thickness is determined only by capillary number and the effect of inertia force is negligible at small capillary numbers. However, the effect of inertia force cannot be neglected as capillary number increases. At relatively high capillary numbers, liquid film thickness takes a minimum value against Reynolds number. The effects of bubble length, liquid slug length and gravity on the liquid film thickness are also investigated. Experimental correlation for the initial liquid film thickness based on capillary number, Reynolds number and Weber number is proposed.     

Drag models for simulating gas-solid flow in the turbulent fluidization of FCC particles

in the turbulent fiuidization of FCC particles, and was validated by satisfactory agreement between prediction and experiment.     

Research on the effect of cylinder particles on the turbulent properties in particulate flows

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｔｗｏ＿ｐｈａｓｅｆｌｏｗｗｉｔｈｃｙｌｉｎｄｅｒｐａｒｔｉｃｌｅｓｉｓａｄｉｆｆｉｃｕｌｔｍａｔｔｅｒｉｎｔｈｅｏｒｅｔｉｃａｌｒｅｓｅａｒｃｈｏｆｍｕｌｔｉｐｈａｓｅｆｌｏｗａｎｄｎｏｎ＿Ｎｅｗｔｏｎｉａｎｆｌｏｗ ,ｓｏｉｔｉｓａｃａｄｅｍｉｃａｌｌｙｖａｌｕａｂｌｅｔｏｓｔｕｄｙｉｔ.Ｍｅａｎｗｈｉｌｅ ,ｓｕｃｈｆｌｏｗｓｈａｖｅａｂｒｏａｄｂａｃｋｇｒｏｕｎｄｉｎｉｎｄｕｓｔｒｙ .Ｉｎｔｈｅｍａｔｅｒｉａｌｓｃｉｅｎｃｅ,ｔｈｅｍｏｌｄｉｎｇａｎｄｐｒｏｃｅｓ…     

TRPIV investigation of space-time correlation in turbulent flows over flat and wavy walls

The present experimental work is devoted to in- vestigate a new space-time correlation model for the turbulent boundary layer over a flat and a wavy walls. A turbulent boundary layer flow at Reo = 2460 is measured by tomographic time-resolved particle image velocimetry （Tomo-TRPIV）. The space-time correlations of instantaneous streamwise fluctuation velocity are calculated at 3 different wall-normal locations in logarithmic layer. It is found that the scales of coherent structure increase with moving far away from the wall. The growth of scales is a manifestation of the growth of prevalent coherent structures in the turbulent boundary layer like hairpin vortex or hairpin packets when they lift up. The resulting contours of the space-time correlation exhibit elliptic-like shapes rather than straight lines. It is suggested that, instead of Taylor hypothesis, the elliptic model of the space-time correlation is valid for the wallbounded turbulent flow over either a flat wall or a wavy wall. The elliptic iso-correlation curves have a uniform preferred orientation whose slope is determined by the convection velocity. The convection velocity derived from the space-time correlation represents the velocity at which the large-scale eddies carry small-scale eddies. The sweep velocity rep- resents the distortions of the small-scale eddies and is intimately associated with the fluctuation velocity in the logarithmic layer of turbulent boundary layers. The nondimensionalized correlation curves confirm that the elliptic model is more proper for approximating the space-time correlation than Taylor hypothesis, because the latter can not embody the small-scale motions which have non-negligible distortions. A second flow over a wavy wall is also recorded using TRPIV. Due to the combined effect of shear layers and the adverse pressure gradient, the space-time correlation does not show an elliptic-like shape at some specific heights over the wavy wall, but in the outer region of the wavy wallbounded flow, the elliptic model remai     

Gas turbulence modulation in the pneumatic conveying of massive particles in vertical tubes

The present work is concerned with the interaction between large particles and gas phase turbulence. Gas turbulence modulation in these systems is considered to be dominated by a generation mechanism which arises due to the presence of wakes behind particles. Following a recent proposal, a closure for gas turbulence modulation accounting for the effect of wakes is employed within the context of a mathematical model for particle-laden, turbulent flows. The model accounts for particle particle and particle-wall interactions associated with larger particles based on concepts from gas kinetic theory. It is shown that due to the significant flattening of the mean gas velocity profile with the addition of particles, and the corresponding decrease in turbulent energy production, a generation mechanism must be present in order to produce gas velocity fluctuation predictions which are consistent with the experimental measurements, even in the case where the experimental results indicate a net suppression of gas phase turbulence in the presence of particles.     

Numerical analysis of turbulent flow in a rotating U-bend with roughened walls

A numerical analysis has been performed for a developing turbulent flow in a rotating U-bend of strong curvature with rib-roughened walls using an anisotropic turbulent model. In this calculation, an algebraic Reynolds stress model is used to precisely predict Reynolds stresses, and a boundary-fitted coordinate system is introduced as a method of coordinate transformation to set the exact boundary conditions along the complicated shape of U-bend with rib-roughened walls. Calculated results for mean velocity and Reynolds stresses are compared to the experimental data in order to validate the proposed numerical method and the algebraic Reynolds stress model. Although agreement is certainly not perfect in all details, the present method can predict characteristic velocity profiles and reproduce the separated flow generated near the outer wall, which is located just downstream of the curved duct. The Reynolds stresses predicted by the proposed turbulent model agree well with the experimental data, except in regions of flow separation.     

On the Direct and Radiated Components of the Collisional Particle Pressure in Liquid-Solid Flows

In a recent study the collisional particle pressure was measured for liquid fluidized beds and liquid-solid flows. The particle pressure was defined as the additional pressure generated by the presence of the particulate-solid phase in a liquid-solid mixture. The particle pressure generated by collisions of particles was found to be composed of two main contributions: one from pressure pulses generated by direct collisions of particles against the containing walls (direct component), and a second one from pressure pulses due to collisions between individual particles that are transmitted through the liquid (radiated component). This paper presents a summary of the technique to measure the particle pressure and the main results of that study.Additional experiments were performed to further study each one of the components of the particle pressure. The direct component was studied by impacting particles on the active face of the pressure transducer. The magnitude of the measured impulse was found to be related to the impact velocity, the mass and the size of the impacting particle. By comparing the measurements with the predictions from Hertzian theory, a quantification of the interstitial fluid effects could be obtained. The radiated component was investigated by generating binary collisions of particles in the vicinity of the transducer. The magnitude of the measured impulse was found to be a function of fluid density, particle size and impact velocity. Predictions based on impulse-pressure theory were obtained and compared with the experimental measurements. The model results showed good agreement with the experimental measurements.     

Microcavitation in the slow motion of a solid spherical particle along a wall in a fluid

In the slow motion of a spherical particle along a wall, a cavity in the form of a non-spherical bubble is detected visually in the narrow particle-wall clearance. The nondimensional parameters associated with cavity formation are formulated and the effect of these parameters on the particle motion is investigated.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, 2004, pp. 110–118. Original Russian Text Copyright © 2004 by Prokunin.     

Experimental studies on the dynamic stability of liquid in a spherical tank covered with diaphragm under vertical excitation

Experimental studies are conducted on the liquid sloshing characteristics in a spherical tank covered with a flexible diaphragm. A spherical acrylic tank with 145.2 mm radius is used as a test tank, and is filled with water. Silicon diaphragms, plane or hemispherical type, with 0.2 mm thickness are used as test diaphragms. The test tank is harmonically excited in the vertical direction by an electro-dynamic exciter. During the test, vibrations due to parametric instability occur when the excitation frequency is twice the natural frequency. Parametric instability regions for some natural modes are measured and are presented in the excitation frequency–excitation acceleration diagram for three cases: liquid surface is uncovered (i.e., free surface), covered with a plane diaphragm, and covered with a hemispherical diaphragm, with the volume of filling water being changed appropriately.