Flow structures around an equilateral triangle arrangement of three spheres

This paper represents the results of an experimental study on the flow structure around a single sphere and three spheres in an equilateral-triangular arrangement. Flow field measurements were performed using a Particle Image Velocimetry (PIV) technique and dye visualization in an open water channel for a Reynolds number of Re = 5 × 10³ based on the sphere diameter. The distributions and flow features at the critical locations of the contours of the velocity fluctuations, the patterns of sectional streamlines, the vorticity contours, the turbulent kinetic energy, the Reynolds stress correlations and shedding frequency are discussed. The gap ratios (G/D) of the three spheres were varied in the range of 1.0 ⩽ G/D ⩽ 2.5 where G was the distance between the sphere centers, and D was the sphere diameter which was taken as 30 mm. Due to the interference of the shedding shear layers and the wakes, more complex features of the flow patterns can be found in the wake region of the two downstream spheres behind the leading sphere. For G/D = 1.25, a jet-like flow around the leading sphere through the gap between the two downstream spheres occurred, which significantly enhanced the wake region. It was observed that a continuous flow development involving shearing phenomena and the interactions of shedding vortices caused a high rate of fluctuations over the whole flow field although most of the time-averaged flow patterns were almost symmetric about the two downstream spheres.     

Analytical solution of two intersecting cylindrical shells subjected to transverse moment on nozzle

In this paper a theoretical solution for two normally intersecting cylindrical shells subjected to transverse moment on the branch pipe is presented, which based on thin shell theory. The accurate shell equations, boundary conditions and calculating methods are adopted so that the solution presented can be applicable up to d/D⩽0.8 and λ=d/(DT)^1/2⩽8. The presented results are in very good agreement with experimental and numerical results for ORNL-1 Model. They are also in agreement with the results obtained by WRC Bulletin 297 when d/D is small.     

Heat transfer due to impinging co-axial jets and the jets’ fluid flow characteristics

This investigation had multiple goals. One goal was to obtain definitive information about the heat transfer characteristics of co-axial impinging jets, and this was achieved by measurements of the stagnation-point, surface-distribution and average heat transfer coefficients. These results are parameterized by the Reynolds number Re which ranged from 5000 to 25,000, the dimensionless separation distance between the jet exit and the impingement plate H/D (4–12), and the ratio of the inner diameters of the inner and outer pipes d/D (0–0.55). The d/D = 0 case corresponds to a single circular jet. The other major goal of this work was to quantify the velocity field of co-axial free jets (impingement plate removed). The velocity-field study included both measurements of the mean velocity and the turbulence intensity.It was found that the variation of the stagnation-point heat transfer coefficient with d/D attained a maximum at d/D = 0.55. Furthermore, the variation of the local heat transfer coefficient across the impingement surface was more peaked for d/D = 0 and became flatter with decreasing d/D. This suggests that for cooling a broad expanse of surface, co-axial jets of high d/D are preferable. On the other hand, for localized cooling, the single jet (d/D = 0) performed the best. In general, for a given Reynolds number, a co-axial jet yields higher heat transfer coefficients than a single jet. Off-axis velocity peaks were encountered for the jets with d/D = 0.105. The measurements of turbulence intensity yielded values as high as 18%.     

Low Reynolds number flow over a square cylinder with a detached flat plate

A numerical study of the alteration of a square cylinder wake using a detached downstream thin flat plate is presented. The wake is generated by a uniform flow of Reynolds number 150 based on the side length of the cylinder, D. The sensitivity of the near wake structure to the downstream position of the plate is investigated by varying the gap distance (G) along the wake centerline in the range 0 ⩽ G ⩽ 7D for a constant plate length of L = D. A critical gap distance is observed to occur at G_c ≈ 2.3D that indicates the existence of two flow regimes. Regime I is characterised by vortex formation occurring downstream of the gap while for regime II, formation occurs within the gap. By varying the plate length and gap distance, a condition is found where significant unsteady total lift reduction can occur. The root mean square lift reduction is limited by an unsteady stall process on the plate.     

CFD based mini- vs. micro-system delineation in elongated bubble flow regime

Delineation of mini- and micro-scale channels with respect to two-phase flow has been the subject of many research papers. There is no consensus on when the small channel can be characterized as a mini-channel or micro-channel. The idea proposed by this paper is to use the normalized bubble nose radius, liquid film thickness top over bottom ratio, and bubble shape contour, which are found under normal gravity conditions in slug flow through a horizontal adiabatic channel, as the delineation criteria. The input parameters are bubble nose radius and bubble nose velocity as the characteristic length scale and characteristic velocity scale respectively. 3D numerical simulation with ANSYS FLUENT was used to obtain the necessary data. Following CFD practice, a mesh independence study and a numerical model validation against published experimental data were both conducted. Analysis of the numerical simulation results showed that channels with D ⩽ 100 μm can be characterized as a micro-system, while channels with D ⩾ 400 μm belong to mini-systems. The region 200 μm ⩽ D ⩽ 300 μm represents a transition from the micro-scale to mini-scale.     

Simulation of a valveless pump with an elastic tube

A valveless pump consisting of a pumping chamber with an elastic tube was simulated using an immersed boundary (IB) method. The interaction between the motion of the elastic tube and the pumping chamber generated a net flow toward the outlet throughout a full cycle of the pump. The net flow rate of the valveless pump was examined by varying the stretching coefficient (ϕ), bending coefficient (γ), the aspect ratio (l/d) of the elastic tube, and the frequency (f) of the pumping chamber. As the stretching and bending coefficients of the elastic tube increased, the net flow through the valveless pump decreased. Elastic tubes with aspect ratios in the range of 2 ⩽ l/d ⩽ 3 generated a higher flow rate than that generated for tubes with aspect rations of l/d = 1 or 4. As the frequency of the pumping chamber increased, the net flow rate of the pump for l/d = 2 increased. However, the net flow rate for l/d = 3 was nonlinearly related to the pumping frequency due to the complexity of the wave motions. Snapshots of the fluid velocity vectors and the wave motions of the elastic tube were examined over one cycle of the pump to gain a better understanding of the mechanism underlying the valveless pump. The relationship between the average gap in the elastic tube and the average flow rate of the pump was analyzed. A smaller gap in the elastic tube during the expansion mode and a wider gap in the elastic tube during the contraction mode played a dominant role in generating a high average flow rate in the pump, regardless of the stretching coefficient (ϕ), the aspect ratio (l/d) of the elastic tube, or the pumping frequency of the pumping chamber (f).     

The flow field around a micropillar confined in a microchannel

The flow field over a low aspect ratio (AR) circular pillar (L/D = 1.5) in a microchannel was studied experimentally. Microparticle image velocimetry (μPIV) was employed to quantify flow parameters such as flow field, spanwise vorticity, and turbulent kinetic energy (TKE) in the microchannel. Flow regimes of cylinder-diameter-based Reynolds number at 100 ⩽ Re_D ⩽ 700 (i.e., steady, transition from quasi-steady to unsteady, and unsteady flow) were elucidated at the microscale. In addition, active flow control (AFC), via a steady control jet (issued from the pillar itself in the downstream direction), was implemented to induce favorable disturbances to the flow in order to alter the flow field, promote turbulence, and increase mixing. Together with passive flow control (i.e., a circular pillar), turbulent kinetic energy was significantly increased in a controllable manner throughout the flow field.     

The influence of pipe length on thermal statistics computed from DNS of turbulent heat transfer

We present results from direct numerical simulation of turbulent heat transfer in pipe flow at a bulk flow Reynolds number of 5000 and Prandtl numbers ranging from 0.025 to 2.0 in order to examine the effect of streamwise pipe length (πδ ≡ πD/2 ⩽ L ⩽ 12πδ) on the convergence of thermal turbulence statistics. Various lower and higher order thermal statistics such as mean temperature, rms of fluctuating temperature, turbulent heat fluxes, two-point auto and cross-correlations, skewness and flatness were computed and it is found that the value of L required for convergence of the statistics depends on the Prandtl number: larger Prandtl numbers requires comparatively shorter pipe length for convergence of most of the thermal statistics.     

Experimental study of single bubble motion in a liquid metal column exposed to a DC magnetic field

The motion of single Argon bubbles rising in the eutectic alloy GaInSn under the influence of a DC longitudinal magnetic field (parallel to the direction of bubble motion) was examined. The magnetic field strength was varied up to 0.3 T corresponding to a magnetic interaction parameter N (which measures the ratio of electromagnetic forces to inertial forces) slightly greater than 1. The liquid metal was at rest in a cylindrical container. Bubble and liquid velocities were measured using ultrasound Doppler velocimetry (UDV). The measured bubble terminal velocity showed oscillations indicating a zigzag movement of ellipsoidal bubbles. For small bubbles (d_e ⩽ 4.6 mm) an increase of the drag coefficient with increasing magnetic interaction parameter N was observed, whereas for larger bubbles (d_e ⩾ 5.4 mm) the application of the magnetic field reduces the drag coefficient. The measurements revealed a distinct electromagnetic damping of the bubble induced liquid velocity leading to more rectilinear bubble trajectories when the magnetic field is applied. Moreover, significant modifications of the bubble wake structure were observed. Raising of the magnetic field strength caused an enlargement of the eddies in the wake. The Strouhal number decreases with increasing magnetic interaction parameter N.     

Study on premixed combustion in cylindrical micro combustors: Transient flame behavior and wall heat flux

The micro combustor is a key component of the micro thermophotovoltaic (TPV) system. Improving the wall temperature of the micro combustor is an effective way to elevate the system efficiency. An experimental study on the wall temperature and radiation heat flux of a series of cylindrical micro combustors (with a backward-facing step) was carried out. For the micro combustors with d = 2 mm, the regime of successful ignition (under the cold wall condition) was identified for different combustor lengths. Acoustic emission was detected for some cases and the emitted sound was recorded and analyzed. Under the steady-state condition, the effects of the combustor diameter (d), combustor length (L), flow velocity (u₀) and fuel–air equivalence ratio (Ф) on the wall temperature distribution were investigated by measuring the detailed wall temperature profiles. In the case that the micro combustor is working as an emitter, the optimum efficiency was found at Ф ≈ 0.8, independent of the combustor dimensions (d and L) and the flow velocity. Under the experimental conditions employed in the present study, the positions of the peak wall temperature were found to be about 8–11 mm and 4–6 mm from the step for the d = 3 mm and d = 2 mm micro combustors, respectively, which are 8–11 and 8–12 times of their respective step heights. This result suggests that the backward-facing step employed in the combustor design is effective in stabilizing the flame position.     

Surface responses induced by point load or uniform traction moving steadily on an anisotropic half-plane

Surface responses induced by point load or uniform traction moving steadily with subsonic speed on an anisotropic half-plane boundary are investigated. It is found that the effects of the material constant on surface displacements are through matrices L^?1(v) and S(v)L^?1(v), while those on surface stress components are through matrices Ω(v) and Γ(v). Explicit expressions for the elements of these four matrices are expressed in terms of elastic stiffness for general anisotropic materials. The special cases of monoclinic materials with symmetry plane at x₁ = 0, x₂ = 0 and x₃ = 0, and the case for orthotropic materials are all deduced. Results for isotropic material may be recovered from present results. For monoclinic materials with a plane of symmetry at x₃ = 0, two of the elements of matrix Ω(v) are found to be independent of subsonic speed.     

Computational and experimental characterization of a liquid jet plunging into a quiescent pool at shallow inclination

A circular water jet (Re = 1.6 × 10⁵; We = 8.8 × 10³) plunging at shallow angles (θ ≈ 12.5°) into a quiescent pool is investigated computationally and experimentally. A surprising finding from the work is that cavities, of the order of jet diameter, are formed periodically in the impact location, even though the impinging flow is smooth and completely devoid of such a periodicity. Computational prediction of these frequencies was compared with experimental findings, yielding excellent agreement. The region in the vicinity of the impact is characterized by strong churning due to splashing and formation of air cavities. Measured velocity profiles indicate a concentration of momentum beneath the free surface slightly beyond the impact location (X/D_j ≈ 14), with a subsequent shift towards the free surface further downstream of this point (X/D_j ≈ 30). This shift is due primarily to the action of buoyancy on the cavity/bubble population. Comparisons of the mean velocity profile between simulations and experiments are performed, yielding good agreement, with the exception of the relatively small churning flow region. Further downstream (X/D_j ≳ 40), the flow develops mostly due to diffusion and the location of peak velocity coincides with the free surface. In this region, the free surface acts as an adiabatic boundary and restricts momentum diffusion, causing the peak velocity to occur at the free surface.     

Hydrodynamic forces on a rotating sphere

The wake dynamics of a rotating sphere with prescribed rotation axis angles are quantitatively analysed by carrying out numerical simulations at Reynolds numbers of Re = 100, 250 and 300, non-dimensional rotational rates Ω¹ = 0–1 and rotation axis angles α = 0, π/6, π/3 and π/2 measured from the free stream axis. These parameters are the same as those in an earlier study (Poon et al., 2010, Int. J. Heat Fluid Flow) where the instantaneous flow structures were discussed qualitatively. This study extends the findings of the earlier study by employing phase diagrams (C_Lx, C_Ly) and (C_D, C_L) to provide a quantitative analysis of the time-dependent behaviour of the flow structures. At Re = 300 and Ω¹ = 0.05, the phase diagrams (C_Lx, C_Ly) show ‘saw tooth’ patterns for both α = 0 and π/6. The ‘saw tooth’ pattern indicates that the flow structures comprise a higher frequency oscillation component at a Reynolds number of 300 which is not observed until Re ≈ 800 for a stationary sphere. This ‘saw tooth’ pattern disappears as Ω¹ increases. The employment of the phase diagrams also reveals that different flow structures induce different oscillation amplitudes on both lateral force coefficients. With the exception of the vortices formed from a shear layer instability, all other flow regimes show larger fluctuations in C_L than C_D.     

Influence of wall proximity on characteristics of wake behind a square cylinder: PIV measurements and POD analysis

Influence of wall proximity on characteristics of the wake behind a two-dimensional square cylinder was experimentally studied in the present work. A low-speed recirculation water channel was established for the experiment; the Reynolds number based on the free-stream velocity and cylinder width (D) was kept at Re_D = 2250. Four cases with different gap width, e.g., G/D = 0.1, 0.2, 0.4 and 0.8, were chosen for comparison. Two experimental techniques, e.g., the standard PIV with high image-density CCD camera and TR-PIV with a high-speed camera were employed in measuring the wake field, enabling a comprehensive view of the time-averaged wake pattern at high spatial resolution and the instantaneous flow field at high temporal resolution, respectively. For the four cases, the difference in spatial characteristics of the wake in the vicinity of the plane wall was analyzed in terms of the time-averaged quantities measured by the standard PIV, e.g., the streamline pattern, the vector field, the streamwise velocity fluctuation intensity and the reverse-flow intermittency. The proper orthogonal decomposition (POD) method was extensively used to decompose the TR-PIV measurements, giving a close-up view of the energetic POD modes buried in the wake. The low-order flow model of the wake at G/D = 0.8 and 0.4 was constructed by using the linear combination of the first two POD modes and the time-mean flow field, which reflected well the vortex shedding process in the sense of the phase-dependent patterns. The intermittent appearance of the weakly separated region near the wall was found at G/D = 0.4. On going from G/D = 0.8 to 0.4, the remarkable variation of the instantaneous wake in the longitudinal direction confirmed that the wall constraint stretches the vortices in the plane of the wall and transfers the energy to the longitudinal component at the expense of the lateral one.     

Turbulent boundary layers over sparsely-spaced rod-roughened walls

Direct numerical simulations (DNSs) of spatially developing turbulent boundary layers (TBLs) over sparsely-spaced two-dimensional (2D) rod-roughened walls were performed. The rod elements were periodically arranged along the streamwise direction with pitches of p_x/k = 8, 16, 32, 64 and 128, where p_x is the streamwise spacing of the rods, and k is the roughness height. The Reynolds number based on the momentum thickness was varied from Re_θ = 300–1400, and the height of the roughness element was k = 1.5θ_in, where θ_in is the momentum thickness at the inlet. The characteristics of the TBLs, such as the friction velocity, mean velocity, and Reynolds stresses over the rod-roughened walls, were examined by varying the spacing of the roughness features (8 ⩽ p_x/k ⩽ 128). The outer-layer similarity between the rough and smooth walls was established for the sparsely-distributed rough walls (p_x/k ⩾ 32) based on the profiles of the Reynolds stresses, whereas those are not for p_x/k = 8 and 16. Inspection of the interaction between outer-layer large-scale motions and near-wall small-scale motions using two-point amplitude modulation (AM) covariance showed that modulation effect of large-scale motions on near-wall small-scale motions was strongly disturbed over the rough wall for p_x/k = 8 and 16. For p_x/k ⩾ 32, the flow that passed through the upstream roughness element transitioned to a smooth wall flow between the consecutive rods. The strong influence of the surface roughness in the outer layer for p_x/k = 8 and 16 was attributed to large-scale erupting motions by the surface roughness, creating both upward shift of the near-wall turbulent energy and active energy production in the outer layer with little influence on the near-wall region.     

The effect of incidence angle on the mean wake of surface-mounted finite-height square prisms

The mean wake of a surface-mounted finite-height square prism was studied experimentally in a low-speed wind tunnel to explore the combined effects of incidence angle (α) and aspect ratio (AR). Measurements of the mean wake velocity field were made with a seven-hole pressure probe for finite square prisms of AR = 9, 7, 5 and 3, at a Reynolds number of Re = 3.7 × 10⁴, for incidence angles from α = 0° to 45°. The relative thickness of the boundary layer on the ground plane, compared to the prism width, was δ/D = 1.5. As the incidence angle increases from α = 0° to 15°, the mean recirculation zone shortens and the mean wake shifts in the direction opposite to that of the mean lift force. The downwash is also deflected to this side of the wake and the mean streamwise vortex structures in the upper part of the wake become strongly asymmetric. The shortest mean recirculation zone, and the greatest asymmetry in the mean wake, is found at the critical incidence angle of α_critical ≈ 15°. As the incidence angle increases from α = 15° to 45°, the mean recirculation zone lengthens and the mean streamwise vortex structures regain their symmetry. These vortices also elongate in the wall-normal direction and become contiguous with the horseshoe vortex trailing arms. The mean wake of the prism of AR = 3 has some differences, such as an absence of induced streamwise vorticity near the ground plane, which support its classification as lying below the critical aspect ratio for the present flow conditions.     

An experimental study of circular sand–water wall jets

Laboratory experiments were carried out to study the effects of sand particles on circular sand–water wall jets. Mean and turbulence characteristics of sand particles in the sand–water wall jets were measured for different sand concentrations c_o ranging from 0.5% to 2.5%. Effects of sand particle size on the centerline sand velocity of the jets were evaluated for sand size ranging from 0.21 mm to 0.54 mm. Interesting results with the range of measurements are presented in this paper. It was found that the centerline sand velocity of the wall jets with larger particle size were 15% higher than the jets with smaller particle size. Concentration profiles in the vertical direction showed a peak value at x/d = 5 (where x is the longitudinal distance from the nozzle and d is the nozzle diameter) and the sand concentration decreased linearly for x/d > 5. Experimental results showed that the turbulence level enhanced from the nozzle to x/d = 10. For sand–water wall jets with a higher concentration (c_o = 1.5–2.5%), the turbulence intensity became smaller than the corresponding single-phase wall jets by 34% due to turbulent modulation. A modified logarithmic formulation was introduced to model the longitudinal turbulent intensity at the centerline and along the axis of the jet.     

Large-eddy simulation of thermal fatigue in a mixing tee

Temperature fluctuations occur due to thermal mixing of hot and cold streams in the T-junctions of the piping system in nuclear power plants, which may cause thermal fatigue of piping system. In this paper, three-dimensional, unsteady numerical simulations of coolant temperature fluctuations at a mixing T-junction of equal diameter pipes were performed using the large eddy simulation (LES) turbulent model. The experiments used in this paper to benchmark the simulations were performed by Hitachi Ltd. The calculated normalized mean temperatures and fluctuating temperatures are in good agreement with the measurements. The influence of the time-step ranging from 100 Hz to 1000 Hz on the numerical simulation results was explored. The simulation results indicate that all the results with different frequencies agree well with the experimental data. Finally, the attenuation of fluctuation of fluid temperature was also investigated. It is found that, drastic fluctuation occurs within the range of less than L/D = 4.0; the fluctuation of fluid temperature does not always attenuate from the pipe center to the wall due to the continuous generation of vortexes. At the top wall, the position of L/D = 1.5 has a minimum normalized mean temperature and a peak value of root-mean square temperature, whereas at the bottom wall, the position having the same characteristics is L/D = 2.0.     

Innovative iteration technique for nonlinear ordinary differential equations of large deflection problem of circular plates

The large deflection problem of a bending circular plate is reduced to a nonlinear ordinary differential equation, and a type of pseudo-linearization is used to obtain the final solution. An innovative iteration technique is suggested. The initial iteration values of some functions for a given loading are adopted from the previous solution for a lower loading. This will significantly extend the range of solution for the non-dimensional loading. Previously, the solution for the non-dimensional loading Q = 100 was obtained. However, in this paper the solution for the non-dimensional loading Q = 1000 is achieved.     

Heat transfer enhancement in a tube with equilateral triangle cross sectioned coiled wire inserts

The heat transfer and pressure drop were experimentally investigated in a coiled wire inserted tube in turbulent flow regime. The coiled wire has equilateral triangular cross section and was inserted separately from the tube wall. The experiments were carried out with three different pitch ratios (P/D = 1, 2 and 3) and two different ratio of equilateral triangle length side to tube diameter (a/D = 0.0714 and 0.0892) at a distance (s) of 1 mm from the tube wall in the range of Reynolds number from 3500 to 27,000. Uniform heat flux was applied to the external surface of the tube and air was selected as fluid. The experimental results obtained from a smooth tube were compared with those from the studies in literature for validation of experimental set-up. The use of coiled wire inserts leads to a considerable increase in heat transfer and pressure drop over the smooth tube. The Nusselt number rises with the increase of Reynolds number and wire thickness and the decrease of pitch ratio. The highest overall enhancement efficiency of 36.5% is achieved for the wire with a/D = 0.0892 and P/D = 1 at Reynolds number of 3858. Consequently, the experimental results reveal that the best operating regime of all coiled wire inserts is detected at low Reynolds number, leading to more compact heat exchanger.