One-dimensional interfacial area transport of vertical upward bubbly flow in narrow rectangular channel

The design and safety analysis for miniature heat exchangers, the cooling system of high performance microelectronics, research nuclear reactors, fusion reactors and the cooling system of the spallation neutron source targets requires the knowledge of the gas–liquid two-phase flow in a narrow rectangular channel. In this study, flow measurements of vertical upward air–water flows in a narrow rectangular channel with the gap of 0.993 mm and the width of 40.0 mm were performed at seven axial locations by using the imaging processing technique. The local frictional pressure loss gradients were also measured by a differential pressure cell. In the experiment, the superficial liquid velocity and the void fraction ranged from 0.214 m/s to 2.08 m/s and from 3.92% to 42.6%, respectively. The developing two-phase flow was characterized by the significant axial changes of the local flow parameters due to the bubble coalescence and breakup in the tested flow conditions. The existing two-phase frictional multiplier correlations such as Chisholm, 1967, Mishima et al., 1993 and Lee and Lee (2001) were verified to give a good prediction for the measured two-phase frictional multiplier. The predictions of the drift-flux model with the rectangular channel distribution parameter correlation of Ishii (1977) and several existing drift velocity correlations of Ishii, 1977, Hibiki and Ishii, 2003 and Jones and Zuber (1979) agreed well with the measured void fractions and gas velocities. The interfacial area concentration (IAC) model of Hibiki and Ishii (2002) was modified by taking the channel width as the system length scale and the modified IAC model could predict the IAC and Sauter mean diameter acceptably.     

Wake structures and vortex-induced vibrations of a long flexible cylinder—Part 1: Dynamic response

Results showing the dynamic response of a vertical long flexible cylinder vibrating at low mode numbers are presented in this paper. The model had an external diameter of 16 mm and a total length of 1.5 m giving an aspect ratio of about 94, with Reynolds numbers between 1200 and 12 000. Only the lower 40% of its length was exposed to the water current in the flume and applied top tensions varied from 15 to 110 N giving fundamental natural frequencies in the range from 3.0 to 7.1 Hz. Reduced velocities based on the fundamental natural frequency up to 16 were reached. The mass ratio was 1.8 and the combined mass–damping parameter about 0.05. Cross-flow and in-line amplitudes, x–y trajectories and phase synchronisation, dominant frequencies and modal amplitudes are reported. Cross-flow amplitudes up to 0.7 diameters and in-line amplitudes over 0.2 were observed with dominant frequencies given by a Strouhal number of 0.16.     

Measurement method of complex viscoelastic material properties

A measurement technique of viscoelastic properties of polymers is proposed to investigate complex Poisson’s ratio as a function of frequency. The forced vibration responses for the samples under normal and shear deformation are measured with varying load masses. To obtain modulus of elasticity and shear modulus, the present method requires only knowledge of the load mass, geometrical characteristics of a sample, as well as both the amplitude ratio and phase lag of the forcing and response oscillations. The measured data were used to obtain the viscoelastic properties of the material based on a 2D numerical deformation model of the sample. The 2D model enabled us to exclude data correction by the empirical form factor used in 1D model. Standard composition (90% PDMS polymer + 10% catalyst) of silicone RTV rubber (Silastic^® S2) were used for preparing three samples for axial stress deformation and three samples for shear deformation. Comprehensive measurements of modulus of elasticity, shear modulus, loss factor, and both real and imaginary parts of Poisson’s ratio were determined for frequencies from 50 to 320 Hz in the linear deformation regime (at relative deformations 10^?6 to 10^?4) at temperature 25 °C. In order to improve measurement accuracy, an extrapolation of the obtained results to zero load mass was suggested. For this purpose measurements with several masses need to be done. An empirical requirement for the sample height-to-radius ratio to be more than 4 was found for stress measurements. Different combinations of the samples with different sizes for the shear and stress measurements exhibited similar results. The proposed method allows one to measure imaginary part of the Poisson’s ratio, which appeared to be about 0.04–0.06 for the material of the present study.     

Application of DEM modified with enlarged particle model to simulation of bead motion in a bead mill

We applied the discrete element method (DEM) of simulation modified by an enlarged particle model to simulate bead motion in a large bead mill. The stainless-steel bead mill has inner diameter of 102 mm and mill length of 198 mm. The bead diameter and filling ratio were fixed respectively at 0.5 mm and 85%. The agitator rotational speed was changed from 1863 to 3261 rpm. The bead motion was monitored experimentally using a high-speed video camera through a transparent mill body. For the simulation, enlarged particle sizes were set as 3–6 mm in diameter. With the DEM modified by the enlarged particle model, the motion of enlarged particles in a mill was simulated. The velocity data of the simulated enlarged particles were compared with those obtained in the experiment. The simulated velocity of the enlarged particles depends on the virtual frictional coefficient in the DEM model. The optimized value of the virtual frictional coefficient can be determined by considering the accumulated mean value. Results show that the velocity of the enlarged particles simulated increases with an increase in the optimum virtual frictional coefficient, but the simulated velocity agrees well with that determined experimentally by optimizing the virtual frictional coefficient in the simulation. The computing time in the simulation decreases with increased particle size.     

Multiple-zone sliding contact with friction on an anisotropic thermoelastic half-space

The paper studies a class of multiple-zone sliding contact problems. This class is general enough to include frictional and thermal effects, and anisotropic response of the indented material. In particular, a rigid die (indenter) slides with Coulomb friction and at constant speed over the surface of a deformable and conducting body in the form of a 2D half-space. The body is assumed to behave as a thermoelastic transversely isotropic material. Thermoelasticity of the Green–Lindsay type is assumed to govern. The solution method is based on integral transforms and singular integral equations. First, an exact transform solution for the auxiliary problem of multiple-zone (integer n > 1) surface tractions is obtained. Then, an asymptotic form for this auxiliary problem is extracted. This form can be inverted analytically, and the result applied to sliding contacts with multiple zones. For illustration, detailed calculations are provided for the case of two (n = 2) contact zones. The solution yields the contact zone width and location in terms of sliding speed, friction, die profile, and also the force exerted. Calculations for the hexagonal material zinc illustrate effects of speed, friction and line of action of the die force on relative contact zone size, location of maximal values for the temperature and the compressive stress, and the maximum temperature for a given maximum stress. Finally, from our general results, a single contact zone solution follows as a simple limit.     

Spatially resolved wall temperature measurements during flow boiling in microchannels

Spatial and temporal variations of channel wall temperature during flow boiling microchannel flows using infrared thermography are presented and analyzed. In particular, the top channel wall temperature in a branching microchannel silicon heat sink is measured non-intrusively. Using this technique, time-averaged temperature measurements, with a spatial resolution of 10 μm, are presented over an 18 mm × 18 mm area of the heat sink. Also presented, within a specific sub-region of the heat sink, are intensity maps that are recorded at a rate of 120 frames per second. Time series data at selected point locations in this sub-region are analyzed for their frequency content, and dominant temperature fluctuations are extracted using proper orthogonal decomposition.Results at low-vapor-quality boiling condition indicate that temperatures can be determined from recorded radiation intensities with a temperature uncertainty varying from 0.9 °C at 25 °C to 1.0 °C at 125 °C. The time series data indicate periodic wall temperature fluctuations of approximately 2 °C that are attributed to the passage of vapor slugs. A dominant band of frequencies around 2–4 Hz is suggested by the frequency analysis. Proper orthogonal decomposition results indicate that first six orthogonal modes account for approximately 90% of the variance in temperature. The first mode reconstruction accounts for temporal variations in the dataset in the sub-region analyzed; however the magnitude of fluctuations and spatial variations in temperature are not accurately captured. A reconstruction using the first 25 modes is considered sufficient to capture both the temporal and spatial variations in the data.     

Experimental analysis of the quasi-static and dynamic torsional behaviour of shape memory alloys

This paper investigates experimentally the quasi-static and dynamic torsional behaviour of shape memory alloys wires under cyclic loading. A specifically designed torsional pendulum made of a Ni–Ti wire is described. Results on the quasi-static behaviour of the wire obtained using this setup are presented, giving an overall view of the damping capacity of the material as function of the amplitude of the loading (imposed torsional angle), the frequency and the temperature. The dynamical behaviour is then presented through measured frequency response function between forcing angle at the top of the pendulum and the difference between top and bottom rotation angles in the vicinity of the first eigenfrequency of the wire, i.e. in the range [0.3 Hz, 1 Hz]. The softening-type non-linearity and its subsequent jump phenomenon, predicted theorically by the decrease of the effective stiffness when martensite transformation starts is clearly evidenced and analysed.     

Large-eddy simulations of chevron jet flows with noise predictions

Hybrid large-eddy type simulations for chevron nozzle jet flows are performed at Mach 0.9 and Re = 1.03 × 10⁶. Without using any subgrid scale model (SGS), the numerical approach applied in the present study is essentially implicit large-eddy simulation (ILES). However, a Reynolds-averaged Navier–Stokes (RANS) solution is patched into the near wall region. This makes the overall solution strategy hybrid RANS–ILES. The disparate turbulence length scales, implied by these different modeling approaches, are matched using a Hamilton–Jacobi equation. The complex geometry features of the chevron nozzles are fully meshed. With numerical fidelity in mind, high quality, hexahedral multi-block meshes of 12.5 × 10⁶ cells are used. Despite the modest meshes, the novel RANS–ILES approach shows encouraging performance. Computed mean and second-order fluctuating quantities of the turbulent near field compare favorably with measurements. The radiated far-field sound is predicted using the Ffowcs Williams and Hawkings (FW–H) surface integral method. Encouraging agreement of the predicted far-field sound directivity and spectra with measurements is obtained.     

Void fraction and pressure drop during external upward two-phase crossflow in tube bundles – part I: Experimental investigation

The present paper is the part I of a broad study concerning void fraction and pressure drop for air-water upward external flow across tube bundles. Experimental results were obtained for liquid and gas superficial velocities ranging from 0.02 to 1.50 m/s and 0.20 to 10.00 m/s, respectively. Void fraction measurements were performed for bubbly flow using a capacitive probe. The test section consisted of a triangular tube bundle counting with 19 mm OD tube and transverse pitch of 24 mm. Initially, the paper describes the test facility and the data regression and experimental procedures. Then, the pressure drop and void fraction measurements are validated based on tests for single-phase flow and quiescent liquid conditions, respectively. Finally, the experimental data are presented and analyzed. In the second part of this study (Part II), a literature review on predictive methods for void fraction and pressure drop is presented. Additionally, these methods are compared with the database presented in Part I and new predictive methods for void fraction and frictional pressure drop are proposed.     

Effect of magnetic field on 3D flow and heat transfer during solidification from a melt

We present the effect of a magnetic field on three-dimensional fluid flow and heat transfer during solidification from a melt in a cubic enclosure. The walls of the enclosure are considered perfectly electrically conducting and the magnetic field is applied separately in three directions. The finite-volume method with enthalpy formulation is used to solve the mathematical model in the solid and liquid phases. The results obtained by our computer code are compared with the numerical and experimental data found in the literature. For Gr = 5 × 10⁵ and Ha = 0, 25, 50, 75, and 100 (where Gr and Ha are the Grashof and Hartmann numbers, respectively), the effects of magnetic field on flow and thermal fields, and on solid/liquid interface shape are presented and discussed. The interface is localized with and without magnetic field. The results show a strong dependence between the interface shape and the intensity and orientation of magnetic field. When the magnetic field is applied along the X-direction, the magnetic stability diagrams (V_max–Ha) and (Nu_avg–Ha) show the strongest stabilization of the flow field and heat transfer.     

The effects of channel diameter on flow pattern,void fraction and pressure drop of two-phase air–water flow in circular micro-channels

Two-phase air–water flow characteristics are experimentally investigated in horizontal circular micro-channels. Test sections are made of fused silica. The experiments are conducted based on three different inner diameters of 0.53, 0.22 and 0.15 mm with the corresponding lengths of 320, 120 and 104 mm, respectively. The test runs are done at superficial velocities of gas and liquid ranging between 0.37–42.36 and 0.005–3.04 m/s, respectively. The flow visualisation is facilitated by systems mainly including stereozoom microscope and high-speed camera. The flow regime maps developed from the observed flow patterns are presented. The void fractions are determined based on image analysis. New correlation for two-phase frictional multiplier is also proposed for practical applications.     

Prediction of symmetry during intermittent and annular horizontal two-phase flows

An optical method is developed for a 2.95 mm inner diameter circular mini-channel to estimate liquid film thicknesses. The greyscale pictures obtained with a high-speed camera are processed to determine the liquid-vapour interface positions for annular and intermittent flows. The experiments are thus performed for a large range of flow conditions. The saturation temperatures tested ranged from 20°C to 100°C in steps of 10°C and the mass velocities are 50, 100, 200, 300 and 400 kg m⁻² s⁻¹. A new parameter ranging from 0 to 1, the symmetry, is defined to account for the level of non uniformity of liquid distribution around the tube perimeter. New experimental data are presented (270 data points), that cover a range of symmetry parameter from 0.35 to 1.00. These data are merged with those available in the literature (406 data points with symmetry parameter from 0.71 to 1.00). A sensitivity analysis of the dimensionless numbers of major influence is run and a new correlation is proposed, that enables to predict over 90% of the data points in an error bandwidth of 10%. This correlation is proposed as criterion for asymmetry.     

Use of zero-net-mass-flow for separation control in diffusing S-duct

Flow control using zero-net-mass-flow jets in an S-shaped diffusing duct was investigated. Experiments were conducted in a channel flow facility at a Reynolds number, Re = 4.1 × 10⁴ with particle image velocimetry measurements in the symmetry plane of the duct. In the natural configuration, separation of the boundary layer occurs in a region of the duct with an high degree of curvature. A stability analysis of the wall normal base flow at the location of the applied control is presented and estimates the most effective frequency of the actuator. Time-averaged velocity fields show total reattachment of the boundary layer using active flow control.     

Positive frictional pressure gradient in vertical gas-high viscosity oil slug flow

The present study describes the wall shear stress and the falling liquid film behavior in upward vertical slug flow of air and high viscosity oil. The frictional pressure gradient is directly related to the wall shear stress, and it is usually negative (opposite to the overall flow direction). However, in vertical slug flow, the average total wall shear stress of a slug unit may be negative (in the same direction of the overall flow), resulting in a positive frictional pressure gradient. However, this does not mean, by any way, generation of additional energy or violation of the second law of thermodynamics.The positive frictional pressure gradient phenomenon, reasons and required conditions were explained in this paper. A simplified model was developed and validated against recent experimental data of air-high viscosity oil slug flow in a 50.8 mm ID vertical pipe. The oil viscosity was in the range of 127 mPa s to 580 mPa s. Positive frictional pressure gradient appears when the liquid film wall shear stress supersede the wall shear stress in the slug body. The rate of increase of both wall shear stresses (with respect to the mixture Reynolds number) depend, not only, on the mixture Reynolds number but also, highly, on the liquid viscosity.     

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.     

Technique for quantitative mapping of three-dimensional liquid–gas phase boundaries in microchannel flows

A diagnostic technique capable of characterizing interfaces between transparent, immiscible fluids is developed and demonstrated by investigating the morphology of liquid–gas interfaces in an adiabatic two-phase flow through a microchannel of 500 μm × 500 μm square cross section. Water seeded with 0.5 μm-diameter fluorescent polystyrene particles is pumped through the channel, and the desired adiabatic two-phase flow regime is achieved through controlled air injection. The diagnostic technique relies on obtaining particle position data through epifluorescent imaging of the flow at excitation and emission wavelengths of 532 nm and 620 nm, respectively. The particle position data are then used to resolve interface locations to within ±1 μm in the focal plane. By mapping the interface within individual focal planes at various depths within the channel, it is possible to determine the complete liquid–gas interface geometry across the channel cross section in a dynamic flow environment. Utilizing this approach, the liquid–gas phase boundaries of annular flows within a microchannel have been successfully characterized.     

Experiments with a Wire-Mesh Sensor for stratified and dispersed oil-brine pipe flow

Two-phase oil–water flow was studied in a 15 m long horizontal steel pipe, with 8.28 cm internal diameter, using mineral oil (having 830 kg/m³ density and 7.5 mPa s viscosity) and brine (1073 kg/m³ density and of 0.8 mPa s viscosity). Measurements of the holdup and of the cross-sectional phase fraction distribution were obtained for stratified flow and for highly dispersed oil–water flows, applying a capacitive Wire-Mesh Sensor specially designed for that purpose. The applicability of this measurement technique, which uses a circuit for capacitive measurements that is adapted to conductive measurements, where one of the fluids is water with high salinity (mimicking sea water), was assessed. Values for the phase fraction values were derived from the raw data obtained by the Wire-Mesh Sensor using several mixture permittivity models. Two gamma-ray densitometers allowed the accurate measurement of the holdups, which was used to validate the data acquired with the capacitive Wire-Mesh Sensor. The measured time-averaged distribution of the phase fraction over the cross-sectional area was used to investigate the details of the observed two-phase flow patterns, including the interface shape and water height. The experiments were conducted in the multiphase-flow test facility of Shell Global International B.V. in the Netherlands.     

Influence of sudden contractions on in situ volume fractions for oil–water flows in horizontal pipes

Oil–water two-phase flow experiments were conducted in horizontal ducts made of Plexiglas® to determine the in situ oil fraction (holdup) by means of the closing valves technique, using mineral oil (viscosity: 0.838 Pa s at 20 °C; density: 890 kg m⁻³) and tap water. The ducts present sudden contractions from 50 mm to 40 mm i.d. and from 50 mm to 30 mm i.d., with contraction ratios of 0.64 and 0.36, respectively. About 200–320 tests were performed by varying the flow rates of the phases. Flow patterns were investigated for both the up- and downstream pipe in order to assess whether relevant variations of the flow patterns across the sudden contraction take place. Data were then compared with predictions of a specific correlation for oil–water flow and some correlations for gas–water flow. A drift-flux model was also applied to determine the distribution parameter.     

High efficiency fuel sleds for polar traverses

We describe here the evolution of lightweight, high-efficiency fuel sleds for Polar over-snow traverses. These sleds consist of flexible bladders strapped to sheets of high molecular weight polyethylene. They cost 1/6th, weigh 1/10th and triple the fuel delivered per towing tractor compared with steel sleds. An eight-tractor fleet has conducted three 3400-km roundtrips to South Pole with each traverse delivering ～320,000 kg of fuel while emitting <1% the pollutants, consuming 1/2 the fuel and saving ～$1.6 M compared with aircraft resupply. A two-tractor fleet in Greenland recently delivered ～83,000 kg of fuel in bladder sleds to Summit with similar benefits. Performance monitoring has revealed that bladder-sled towing resistance is largely governed by sliding friction, which can start high and drop in half over the first 30 min of travel. Frictional heating probably produces a thin water layer that lubricates the sled–snow interface. Consequently, towing resistance depends on the thermal budget of the sled. For example, black fuel bladders increase solar gain and thus decrease sled resistance; data suggest they could double again the fuel delivered per tractor. The outstanding efficiency and low cost of these sleds has transformed fuel delivery to Polar research stations.     

Hydrodynamics of vertical falling films in a large-scale pilot unit – a combined experimental and numerical study

The hydrodynamics of vertical falling films in a large-scale pilot unit are investigated experimentally and numerically. We study a broad range of operating conditions with Kapitza and Reynolds numbers ranging from Ka = 191–3394 and Re 24–251, respectively. We compare film thickness measurements, conducted by a laser triangulation scanner, with those obtained by directly solving the full Navier–Stokes equations in two dimensions and using the volume of fluid (VOF) numerical framework. We examine the evolution of the liquid film at multiple locations over a vertical distance of 4.5 m. In both our experiments and simulations we identify a natural wave frequency of the system of approximately 10 Hz. We investigate the formulation of the inlet boundary condition and its effects on wave formation. We show how potentially erroneous conclusions can be made if the simulated domain is shorter than 1000 film thicknesses, by mistaking the forced inlet frequency for the natural wave frequency. We recommend an inlet disturbance consisting of a multitude of frequencies to achieve the natural wave frequency over relatively short streamwise distances.