Turbulent transport of particles in a straight square duct

Turbulent flow through a duct of square cross-section gives rise to off-axis secondary flows, which are known to transfer momentum between fluid layers thereby flattening the velocity profile. The aim of this study is to investigate the role of the secondary flows in the transport and dispersion of particles suspended in a turbulent square duct flow. We have numerically simulated a flow through a square duct having a Reynolds number of Re_τ = 300 through discretization of the Navier–Stokes equations, and followed the trajectories of a large number of passive tracers and finite-inertia particles under a one-way coupling assumption. Snapshots of particle locations and statistics of single-particle and particle pair dispersion were analyzed. It was found that lateral mixing is enhanced for passive tracers and low-inertia particles due to the lateral advective transport that is absent in straight pipe and channels flows. Higher inertia particles accumulate close to the wall, and thus tend to mix more efficiently in the streamwise direction since a large number of the particles spend more time in a region where the mean fluid velocity is small compared to the bulk. Passive tracers tend to remain within the secondary swirling flows, circulating between the core and boundary of the duct.     

Initial section in a square duct rotating about the transverse axis

Laminar flow in a duct uniformly rotating about the transverse axis has been the subject of a number of studies. These are reviewed, for example, in [1]. However, in all these studies only the case of stablized flow was subjected to theoretical investigation. In this paper the development of laminar flow in a prismatic duct uniformly rotating about an axis perpendicular to one of the walls is analyzed on the basis of the results of a numerical integration of the parabolized system of Navier-Stokes equations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 41–46, September–October, 1985.     

Temperature statistics in a radiatively heated particle-laden turbulent square duct flow

Radiation absorption by preferentially concentrated particles in a turbulent square duct flow is studied experimentally. The particle-laden flow is exposed to near-infrared radiation, and the gas phase temperature statistics are measured along the wall bisector of the duct. It is found that the instantaneous temperature fluctuations are comparable to the overall mean temperature rise. The temperature statistics at the duct centerline and near the wall are qualitatively different. The former reflects preferential concentration in isotropic flows while the latter displays evidence of particle clustering into streamwise elongated streaks. Comparison of the experimental data to a simplified heat transfer model suggests that the Lagrangian evolution of particle clusters and voids, and turbulent mixing in the vicinity of particle clusters, are important. This work was motivated by particle solar receiver technology, but the findings are also relevant to systems where there is localized heat release or mass transfer from disperse particles or droplets. It shows that obtaining Lagrangian histories of particle trajectories is an important next step towards understanding thermal transport phenomena in particle-laden turbulent flows.     

Turbulent flow measurement in a square duct with hybrid holographic PIV

 A hybrid holographic system has been developed for three-dimensional particle image velocimetry. With unique high pass filters, the system combines advantages of both in-line and off-axis holography without having their draw-backs. It improves the signal to noise ratio of the reconstructed image, allows use of 3–15 μm particles in water at high population and achieves large dynamic ranges in both velocity and space. With an automated image acquisition and processing system it has been used for measuring the velocity distributions in a square duct at Re=1.23×10⁵. The data consists of 97×97×87 vectors (with 50% overlapping of adjacent interrogation windows). The quality of the results is evaluated using the continuity equation. The deviation from the equation decreases rapidly with increasing control volume and reaches a level of less than 10%. Mean velocities, r.m.s. velocity fluctuations and turbulence spectra are estimated using the data. Received: 16 December 1996/Accepted: 6 March 1997     

The structure of water-air bubble grid turbulence in a square duct

Local measurements of void fraction and continuous phase velocity field in water-air bubble, grid turbulence were conducted in a channel of vertical, square test section. The measured statistics indicate that, due mainly to the interaction of mean shear with the dispersed phase, the turbulence structure of the flow is modified. The observed change is characterized by a strong spatial dependence of void fraction and liquid flow properties, and the emergence of two spatial regions controlled by different physical processes. Intensity measurements indicate significant departure from isotropy in the flow. Two distinct regimes corresponding to low and high values of void fraction have been also identified. The autocorrelation and spectra measurements indicate that for low void fraction the scales of turbulence decrease while for higher values of void fraction increase again and inverse cascade is observed.     

Measurements of turbulent developing flow in a moderately curved square duct

Laser-Doppler measurements in the turbulent flow in a right-angled bend of square cross-section, radius/duct-width ratio 7.0, are presented and show the development of secondary circulation in cross-stream planes. Distribution of the streamwise and radial components of the mean velocity and turbulence intensity, and the corresponding Reynolds shear stress, are presented as contour plots and are intended for use in the further development of numerical flow prediction methods.     

Flow of second-order fluid in a curved duct with square cross-section

This paper investigates the inertial and creeping flow of a second-order fluid in a curved duct with a square cross-section. Numerical modeling is employed to analyze fluid flow, and the governing equations are discretized using the finite difference method on a staggered mesh. The marker-and-cell method is employed to allocate the parameters on the staggered mesh, and static pressure is calculated using the artificial compressibility approach. The effect of centrifugal force due to the curvature of the duct and the opposing effects of the first and second normal stress difference on the flow field are investigated. In addition, the order-of-magnitude technique is used to derive the force balance relations for the core region of flow. Based on these relations, the performance mechanism of centrifugal force and normal stress differences on the generation of secondary flows is considered. We also present an analytical relation for the axial velocity profile and flow resistance ratio of creeping flow. For this kind of flow, previous studies have investigated the effect of the first normal stress difference on the transition from one pair to two pairs of vortices while we show that the negative second normal stress difference has the opposite effect on this transition and can stabilize the flow.     

The near-corner mass transfer associated with turbulent flow in a square duct

Mass transfer to a nearly fully developed turbulent flow in a square duct is studied. Measurements are taken either on one or two adjacent walls in the region of a developing mass transfer boundary layer. Effects of the secondary flow on mass transfer are observed near the corners. The effects are initially small, but grow to a significant degree further downstream. The location of the maximum and its lateral spread is affected by the Reynolds number, as well as the adjacent wall mass transfer.

---

Stoffübertragung in der Nähe einer Ecke bei turbulenter Strömung in einem quadratischen Rohr

Zusammenfassung Es wird über den Massentransport in einer nahezu vollentwickelten turbulenten Strömung in einem Kanal mit quadratischem Querschnitt berichtet. Messungen werden an einer oder an zwei benachbarten Wänden im Bereich der sich entwickelnden Massenübergangsgrenzschicht durchgeführt. Im Bereich der Kanalecken werden Effekte der Sekundärströmung beobachtet. Diese Effekte sind anfänglich gering, wachsen jedoch weiter stromabwärts beträchtlich an. Die Lage des Maximums und seine laterale Ausbreitung werden sowohl durch die Reynoldszahl als auch durch den benachbarten Massenübergang zur Wand beeinflußt.

---

Nomenclature D Duct's span - Dh Hydraulic diameter - Re Reynolds number (U _bulk Dh/v) - Sh Sherwood number - U Axial velocity - U _bulk Bulk velocity - x, y, z Coordinate system, Fig. 1 - x Axial distance from the duct's inlet - Kinematic viscosity Dedicated to Prof. Dr.-Ing. U. Grigull's 80th birthday     

Numerical calculation of Reynolds stresses in a square duct with secondary flow

A new model for the Reynolds stress equations is presented. This model is used to obtain a theoretical solution for the problem of fully developed turbulent flow in a square duct. Nine governing equations for the axial velocity, lateral vorticity, lateral stream function and six components of the Reynolds stresses are simultaneously solved, by a finite-difference technique. To ensure numerical stability of the solution a special linearised implicit representation of the source terms is proposed, and simultaneous solution of the equations at each.mesh point is obtained. Near the wall a special procedure is used, by which the Reynolds stress equations are assumed to be in local equilibrium, and the velocity profile is assumed to be logarithmic. However, due to the secondary motion the logarithmic velocity profile is inclined to the axial direction. The results bear reasonable agreement with experimental data. Computer time requirements are moderate.     

Numerical simulation of magnetohydrodynamic flow in a toroidal duct of square cross-section

We present numerical simulation results of the quasi-static magnetohydrodynamic (MHD) flow in a toroidal duct of square cross-section with insulating Hartmann walls and conducting side walls. Both laminar and turbulent flows are considered. In the case of steady flows, we present a comprehensive analysis of the secondary flow. It consists of two counter-rotating vortex cells, with additional side wall vortices emerging at sufficiently high Hartmann number. Our results agree well with existing asymptotic analysis. In the turbulent regime, we make a comparison between hydrodynamic and MHD flows. We find that the curvature induces an asymmetry between the inner and outer side of the duct, with higher turbulence intensities occurring at the outer side wall. The magnetic field is seen to stabilize the flow so that only the outer side layer remains unstable. These features are illustrated both by a study of statistically averaged quantities and by a visualization of (instantaneous) coherent vortices.     

Analytic solution for secondary flow with heat transfer in a square duct

Analytic solutions are obtained by solving the conservation equations describing fully developed parabolic laminar flow with heat transfer to predict the axial velocity profile, the temperature distribution and the secondary flow behaviour in a square duct. The biharmonic equation arising in this problem has been solved by using the Navier method for simply supported rectangular plates. The analytic solutions substantially agree with the available experimental results.     

Experimental investigation of local heat transfer in a square duct with various-shaped ribs

In the present study, experimental studies are carried out to investigate the heat transfer and friction characteristics in a square duct roughened by various-shaped ribs on one wall. The ribs are oriented transversely to the main stream in a periodic arrangement. Liquid crystal thermography is employed to measure the local and average heat transfer coefficient on the ribbed surface. The rib height-to-duct hydraulic diameter ratio is fixed at 0.1; the rib pitch-to-height ratio varies from 8 to 15 and the test Reynolds number spans from 8,000 to 20,000. The results show that the trapezoidal-shaped ribs with decreasing height in the flow direction (case C) provide the highest heat transfer enhancement factor and are likely to be used to suppress the local hot spot which usually occurs in the region just behind the ribs.     

Heat transfer enhancement through a square duct fitted with twisted tape inserts

Experimental investigations of friction factor and heat transfer characteristics of a square duct fitted with twisted tapes of different twist ratios have been reported at nearly uniform wall temperature conditions. The experimental results indicate that the friction factor and Nusselt number increases with decreasing twist ratio. The maximum heat transfer enhancement was observed for a minimum twist ratio. The thermohydraulic performance analysis is made to identify potential benefits of using a twisted tape.     

Application of photobleaching molecular tagging velocimetry to turbulent bubbly flow in a square duct

To evaluate turbulence energy budget in bubbly flows, an image processing method in a photobleaching molecular tagging velocimetry is improved for accurate evaluation of velocity gradients. Turbulence properties in single-phase and two-phase dilute-bubbly flows in a square duct are measured using the improved method. As a result, the following conclusions are obtained: (1) The axial velocity and axial turbulent intensity measured by the present method agree well with those measured by laser Doppler velocimetry not only for the single-phase flow but also for the dilute-bubbly flow. (2) The present method can measure velocity components and velocity gradients in the vicinity of the wall, and therefore the present method is of great use in understanding the mechanism of turbulence generation and dissipation near the wall. (3) The present method can provide detailed information on turbulence structure such as turbulence kinetic energy budget. (4) Bubbles tend to increase not only the turbulence production but also the turbulence dissipation.     

Numerical investigation of fully developed turbulent fluid flow and heat transfer in a square duct

Three-dimensional fully developed turbulent fluid flow and heat transfer in a square duct are numerically investigated with the author's anisotropic low-Reynolds-number k-ε turbulence model. Special attenton has been given to the regions close to the wall and the corner, which are known to influence the characteristics of secondary flow a great deal. Hence, instead of the common wall function approach, the no-slip boundary condition at the wall is directly used. Velocity and temperature profiles are predicted for fully developed turbulent flows with constant wall temperature. The predicted variations of both local wall shear stress and local wall heat flux are shown to be in close agreement with available experimental data. The present paper also presents the budget of turbulent kinetic energy equation and the systematic evaluation for existing wall function forms. The commonly adopted wall function forms that are valid for two-dimensional flows are found to be inadequate for three-dimensional turbulent flows in a square duct.     

Impingement cooling in a rotating curved square annular duct with crossflow effect from Rib-roughened surfaces

An experimental investigation was carried out to determine the effects of jets in crossflow on impingement heat transfer from rib-roughened curved duct with rotational speeds of 120 and 240 rpm. The jet impinged on curved surface in crossflow having repeated square ribs (3 2 3mm). The curvature of the duct was fixed and has a value of 300 mm. The rib height (e) was fixed at 1.5 mm and pitch-to-height ratio (p/e) kept at 2. The study covered jet Reynolds number in the range 6500 to 26000, and duct stream crossflow Reynolds number from 3250 to 13000 based on jet nozzle diameter, respectively, which gave M = 0.12 to 2. Results were presented for rotational, crossflow, roughened surfaces and curvature effects on local Nusselt numbers. Significant heat transfer enhancement was found for the present physical geometry and within the ranges of operating parameters considered in the study.     

Experimental studies on mixed convection heat transfer in laminar flow through a plain square duct

This paper reports the findings of experimental studies on combined free and forced convection through a plain square duct in laminar region. The test fluid flows through an inner square duct, hot water at high flow rate circulated through a annular channel formed between square duct and circular tube, in counter current fashion to attain a nearly uniform wall temperature conditions. The importance of mixed convection is judged by the value of the Richardson number (Ri). It was observed that at low Reynolds number, heat transfer was mainly governed by mixed convection. However at higher values of Reynolds number, heat transfer was significantly dominated by forced convection. It was found that Reynolds number higher than 1050 for water and 480 for ethylene glycol resulted in laminar forced convention heat transfer. The empirical correlation developed for Nusselt number in terms of Grashoff number and Graez number, was found to fit with experimental Nusselt number within ±11 and ±12?% for water and ethylene glycol respectively.