Rack Level Modeling of Air Flow Through Perforated Tile in a Data Center

Effective air flow distribution through perforated tiles is required to efficiently cool servers in a raised floor data center. We present detailed computational fluid dynamics (CFD) modeling of air flow through a perforated tile and its entrance to the adjacent server rack. The realistic geometrical details of the perforated tile, as well as of the rack are included in the model. Generally, models for air flow through perforated tiles specify a step pressure loss across the tile surface, or porous jump model based on the tile porosity. An improvement to this includes a momentum source specification above the tile to simulate the acceleration of the air flow through the pores, or body force model. In both of these models, geometrical details of tile such as pore locations and shapes are not included. More details increase the grid size as well as the computational time. However, the grid refinement can be controlled to achieve balance between the accuracy and computational time. We compared the results from CFD using geometrical resolution with the porous jump and body force model solution as well as with the measured flow field using particle image velocimetry (PIV) experiments. We observe that including tile geometrical details gives better results as compared to elimination of tile geometrical details and specifying physical models across and above the tile surface. A modification to the body force model is also suggested and improved results were achieved.

COMPARISON AND CROSS-VALIDATION OF OPTICAL TECHNIQUES IN DIFFERENT SWIRL SPRAY REGIMES

This paper deals with an experimental study of pressure-swirl hydraulic injector nozzles using non-intrusive optical techniques. Experiments were conducted to study atomization characteristics using two nozzles with different orifice diameters, 0.3 mm and 0.5 mm, and injection pressures, 0.3-3.5 Mpa, which correspond to Reynolds number (Re-p) = 7,000-45,000, depending on nozzle utilized. Three laser diagnostic techniques were utilized: Shadowgraph, PIV (Particle Image Velocimetry), and PDPA (Phase Doppler Particle Anemometry). Measurements made in the spray in both axial and radial directions indicate that velocity, average droplet diameter profiles, and spray dynamics are highly dependent on the nozzle characteristics and injection pressure. Limitations of these techniques in the different flow regimes, related to the primary and secondary breakups as well as coalescence, are provided. Results indicate that all three techniques provide similar results throughout the different regimes. Shadowgraph and PDPA were possible in the secondary atomization and coalescence regimes while PIV measurements could be made only at the end of secondary atomization and coalescence.

Flexibility in flapping foil suppresses meandering of induced jet in absence of free stream

Thrust-generating flapping foils are known to produce jets inclined to the free stream at high Strouhal numbers St = fA/U-infinity, where f is the frequency and A is the amplitude of flapping and U-infinity is the free-stream velocity. Our experiments, in the limiting case of St —> infinity (zero free-stream speed), show that a purely oscillatory pitching motion of a chordwise flexible foil produces a coherent jet composed of a reverse Benard-Karman vortex street along the centreline, albeit over a specific range of effective flap stiffnesses. We obtain flexibility by attaching a thin flap to the trailing edge of a rigid NACA0015 foil; length of flap is 0.79 c where c is rigid foil chord length. It is the time-varying deflections of the flexible flap that suppress the meandering found in the jets produced by a pitching rigid foil for zero free-stream condition. Recent experiments (Marais et al., J. Fluid Mech., vol. 710, 2012, p. 659) have also shown that the flexibility increases the St at which non-deflected jets are obtained. Analysing the near-wake vortex dynamics from flow visualization and particle image velocimetry (PIV) measurements, we identify the mechanisms by which flexibility suppresses jet deflection and meandering. A convenient characterization of flap deformation, caused by fluid-flap interaction, is through a non-dimensional effective stiffness', EI* = 8 EI/(rho V-TEmax(2) s(f) c(f)(3)/2), representing the inverse of the flap deflection due to the fluid-dynamic loading; here, EI is the bending stiffness of flap, rho is fluid density, V-TEmax is the maximum velocity of rigid foil trailing edge, s(f) is span and c(f) is chord length of the flexible flap. By varying the amplitude and frequency of pitching, we obtain a variation in EI* over nearly two orders of magnitude and show that only moderate EI*. (0.1 less than or similar to EI * less than or similar to 1 generates a sustained, coherent, orderly jet. Relatively `stiff' flaps (EI* greater than or similar to 1), including the extreme case of no flap, produce meandering jets, whereas highly `flexible' flaps (EI* less than or similar to 0.1) produce spread-out jets. Obtained from the measured mean velocity fields, we present values of thrust coefficients for the cases for which orderly jets are observed.

On flame-turbulence interaction in constant-pressure expanding flames

In this paper we present one of the first high-speed particle image velocimetry measurements to quantify flame-turbulence interaction in centrally-ignited constant-pressure premixed flames expanding in nearisotropic turbulence. Measurements of mean flow velocity and rms of fluctuating flow velocity are provided over a range of conditions both in the presence and absence of the flame. The distributions of stretch rate contributions from different terms such as tangential straining, normal straining and curvature are also provided. It is found that the normal straining displays non-Gaussian pdf tails whereas the tangential straining shows near Gaussian behavior. We have further tracked the motion of the edge points that reside and co-move with the edge of the flame kernel during its evolution in time, and found that within the measurement conditions, on average the persistence time scales of stretch due to pure curvature exceed that due to tangential straining by at least a factor of two. (C) 2014 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Acoustic response of vortex breakdown modes in a coaxial isothermal unconfined swirling jet

The present experimental work is concerned with the study of amplitude dependent acoustic response of an isothermal coaxial swirling jet. The excitation amplitude is increased in five distinct steps at the burner's Helmholtz resonator mode (i.e., 100 Hz). Two flow states are compared, namely, sub-critical and super-critical vortex breakdown (VB) that occur before and after the critical conical sheet breakdown, respectively. The geometric swirl number is varied in the range 2.14-4.03. Under the influence of external pulsing, global response characteristics are studied based on the topological changes observed in time-averaged 2D flow field. These are obtained from high resolution 2D PIV (particle image velocimetry) in the longitudinal-mid plane. PIV results also illustrate the changes in the normalized vortex core coordinates (r(vcc)/(r(vcc))(0) (Hz), y(vcc)/(y(vcc))(0) (Hz)) of internal recirculation zone (IRZ). A strong forced response is observed at 100 Hz (excitation frequency) in the convectively unstable region which get amplified based on the magnitude of external forcing. The radial extent of this forced response region at a given excitation amplitude is represented by the acoustic response region (b). The topological placement of the responsive convectively unstable region is a function of both the intensity of imparted swirl (characterized by swirl number) and forcing amplitude. It is observed that for sub-critical VB mode, an increase in the excitation amplitude till a critical value shifts the vortex core centre (particularly, the vortex core moves downstream and radially outwards) leading to drastic fanning-out/widening of the IRZ. This is accompanied by similar to 30% reduction in the recirculation velocity of the IRZ. It is also observed that b < R (R: radial distance from central axis to outer shear layer-OSL). At super-critical amplitudes, the sub-critical IRZ topology transits back (the vortex core retracts upstream and radially inwards) and finally undergoes a transverse shrinkage ((r(vcc))/(r(vcc))(0 Hz) decreases by similar to 20%) when b >= R. In contrast, the vortex core of super-critical breakdown mode consistently spreads radially outwards and is displaced further downstream. Finally, the IRZ fans-out at the threshold excitation amplitude. However, the acoustic response region b is still less than R. This is explained based on the characteristic geometric swirl number (S-G) of the flow regimes. The super-critical flow mode with higher S-G (hence, higher radial pressure drop due to rotational effect which scales as Delta P similar to rho u theta(2) and acts inwards towards the center line) compared to sub-critical state imposes a greater resistance to the radial outward spread of b. As a result, the acoustic energy supplied to the super-critical flow mode increases the degree of acoustic response at the pulsing frequency and energizes its harmonics (evident from power spectra). As a disturbance amplifier, the stronger convective instability mode within the flow structure of super-critical VB causes the topology to widen/fan-out severely at threshold excitation amplitude. (C) 2015 AIP Publishing LLC.

Deformation field evolution in indentation of a porous brittle solid

An experimental study of plane strain wedge indentation of a model porous brittle solid has been made to understand the effect of indentation parameters on the evolution of the deformation field and the accompanying volume change. A series of high-speed, high-resolution images of the indentation region and simultaneous measurements of load response were captured through the progression of the indentation process. Particle image velocimetry analysis of the images facilitated in situ measurement of the evolution of the resulting plastic zone in terms of incremental material displacement (velocity), strain rate, strain and volume change (e.g., local pore collapse). These measurements revealed initiation and propagation of flow localizations and fractures, as well as enabled estimate of volume changes occurring in the deformation zone. The results were directly compared with theoretical estimates of indentation pressure and deformation zone geometry and were used to validate a modified cavity expansion solution that incorporates effects of volume changes in the plastic zone. (C) 2015 Elsevier Ltd. All rights reserved.

After transition in a soft-walled microchannel

In comparison to the flow in a rigid channel, there is a multifold reduction in the transition Reynolds number for the flow in a microchannel when one of the walls is made sufficiently soft, due to a dynamical instability induced by the fluid-wall coupling, as shown by Verma & Kumaran (J. Fluid Mech., vol. 727, 2013, pp. 407-455). The flow after transition is characterised using particle image velocimetry in the x-y plane, where x is the streamwise direction and y is the cross-stream coordinate along the small dimension of the channel of height 0.2-0.3 mm. The flow after transition is characterised by a mean velocity profile that is flatter at the centre and steeper at the walls in comparison to that for a laminar flow. The root mean square of the streamwise fluctuating velocity shows a characteristic sharp increase away from the wall and a maximum close to the wall, as observed in turbulent flows in rigid-walled channels. However, the profile is asymmetric, with a significantly higher maximum close to the soft wall in comparison to that close to the hard wall, and the Reynolds stress is found to be non-zero at the soft wall, indicating that there is a stress exerted by fluid velocity fluctuations on the wall. The maximum of the root mean square of the velocity fluctuations and the Reynolds stress (divided by the fluid density) in the soft-walled microchannel for Reynolds numbers in the range 250-400, when scaled by suitable powers of the maximum velocity, are comparable to those in a rigid channel at Reynolds numbers in the range 5000-20 000. The near-wall velocity profile shows no evidence of a viscous sublayer for (y upsilon(*)/nu) as low as two, but there is a logarithmic layer for (y upsilon(*)/nu) up to approximately 30, where the von Karman constants are very different from those for a rigid-walled channel. Here, upsilon(*) is the friction velocity, nu is the kinematic viscosity and y is the distance from the soft surface. The surface of the soft wall in contact with the fluid is marked with dye spots to monitor the deformation and motion along the fluid-wall interface. Low-frequency oscillations in the displacement of the surface are observed after transition in both the streamwise and spanwise directions, indicating that the velocity fluctuations are dynamically coupled to motion in the solid.

Quantitative In Situ Analysis of Deformation in Sliding Metals: Effect of Initial Strain State

Using in situ, high-speed imaging of a hard wedge sliding against pure aluminum, and image analysis by particle image velocimetry, the deformation field in sliding is mapped at high resolution. This model system is representative of asperity contacts on engineered surfaces and die-workpiece contacts in deformation and machining processes. It is shown that large, uniform plastic strains of 1-5 can be imposed at the Al surface, up to depths of 500 mu m, under suitable sliding conditions. The spatial strain and strain rate distributions are significantly influenced by the initial deformation state of the Al, e.g., extent of work hardening, and sliding incidence angle. Uniform straining occurs only under conditions of steady laminar flow in the metal. Large pre-strains and higher sliding angles promote breakdown in laminar flow due to surface fold formation or flow localization in the form of shear bands, thus imposing limits on uniform straining by sliding. Avoidance of unsteady sliding conditions, and selection of parameters like sliding angle, thus provides a way to control the deformation field. Key characteristics of the sliding deformation such as strain and strain rate, laminar flow, folding and prow formation are well predicted by finite element simulation. The deformation field provides a quantitative basis for interpreting wear particle formation. Implications for engineering functionally graded surfaces, sliding wear and ductile failure in metals are discussed.

Advancing Edge Speeds of Epithelial Monolayers Depend on Their Initial Confining Geometry

Collective cell migrations are essential in several physiological processes and are driven by both chemical and mechanical cues. The roles of substrate stiffness and confinement on collective migrations have been investigated in recent years, however few studies have addressed how geometric shapes influence collective cell migrations. Here, we address the hypothesis that the relative position of a cell within the confinement influences its motility. Monolayers of two types of epithelial cells-MCF7, a breast epithelial cancer cell line, and MDCK, a control epithelial cell line-were confined within circular, square, and cross-shaped stencils and their migration velocities were quantified upon release of the constraint using particle image velocimetry. The choice of stencil geometry allowed us to investigate individual cell motility within convex, straight and concave boundaries. Cells located in sharp, convex boundaries migrated at slower rates than those in concave or straight edges in both cell types. The overall cluster migration occurred in three phases: an initial linear increase with time, followed by a plateau region and a subsequent decrease in cluster speeds. An acto-myosin contractile ring, present in the MDCK but absent in MCF7 monolayer, was a prominent feature in the emergence of leader cells from the MDCK clusters which occurred every similar to 125 mu m from the vertex of the cross. Further, coordinated cell movements displayed vorticity patterns in MDCK which were absent in MCF7 clusters. We also used cytoskeletal inhibitors to show the importance of acto-myosin bounding cables in collective migrations through translation of local movements to create long range coordinated movements and the creation of leader cells within ensembles. To our knowledge, this is the first demonstration of how bounding shapes influence long-term migratory behaviours of epithelial cell monolayers. These results are important for tissue engineering and may also enhance our understanding of cell movements during developmental patterning and cancer metastasis.

Shear banding in a yield stress bearing Langmuir monolayer

We examine the shear-thinning behaviour of a two dimensional yield stress bearing monolayer of sorbitan tristearate at air/water interface. The flow curve consists of a linear region at low shear stresses/shear rates, followed by a stress plateau at higher values. The velocity profile obtained from particle imaging velocimetry indicates that shear banding occurs, showing coexistence of the fluidized region near the rotor and solid region with vanishing shear-rate away from the rotor. In the fluidized region, the velocity profile, which is linear at low shear rates, becomes exponential at the onset of shear-thinning, followed by a time varying velocity profile in the plateau region. At low values of constant applied shear rates, the viscosity of the film increases with time, thus showing aging behaviour like in soft glassy three-dimensional (3D) systems. Further, at the low values of the applied stress in the yield stress regime, the shear-rate fluctuations in time show both positive and negative values, similar to that observed in sheared 3D jammed systems. By carrying out a statistical analysis of these shear-rate fluctuations, we estimate the effective temperature of the soft glassy monolayer using the Galavatti-Cohen steady state fluctuation relation.

Velocity and rotation measurements in acoustically levitated droplets

The velocity scale inside an acoustically levitated droplet depends on the levitator and liquid properties. Using Particle Imaging Velocimetry (PIV), detailed velocity measurements have been made in a levitated droplet of different diameters and viscosity. The maximum velocity and rotation are normalized using frequency and amplitude of acoustic levitator, and droplet viscosity. The non-dimensional data are fitted for micrometer- and millimeter-sized droplets levitated in different levitators for different viscosity fluids. It is also shown that the rotational speed of nanosilica droplets at an advanced stage of vaporization compares well with that predicted by exponentially fitted parameters. (C) 2012 Elsevier B.V. All rights reserved.

Airblast spray in crossflow - Structure, trajectory and droplet sizing

This study reports results of an experimental investigation of airblast spray of water and ethanol in crossflow. Laser shadowgraphy and Particle/Droplet Imaging Analysis (PDIA) are used to derive spray trajectory and drop size information while Particle Tracking Velocimetry (PTV) is used to measure droplet velocities. A new phenomenon of spray bifurcation is observed for low Gas to Liquid Ratio (GLR) cases. The reasons for the spatial bifurcation can be attributed to a combination of reasons. These are (a) presence of large ligaments and droplets in the near-nozzle region for low GLRs (b) secondary breakup experienced by ligaments/droplets leading to formation of a large number of small droplets, and (c) the crossflow causing differential dispersion of the small and large droplets. A novel correlation for spray trajectory is proposed incorporating the momentum ratio and liquid surface tension. This correlation is shown to be effective in predicting the non-linear spray trajectory over a large range of conditions for not only water but ethanol and Jet-A also. It is observed that the larger droplets penetrate further into the crossflow, in the direction of injection. Thus, with increase in height of the measurement location from the injection plane, the droplet Sauter Mean Diameter (SMD) is found to increase. Moreover, as the droplets travel downstream in the crossflow direction, the droplet SMD is observed to decrease. The effect of drag is assessed by comparing velocity of different sizes of droplets at various locations. Smaller droplets are entrained into the crossflow at much lower elevations, whereas larger droplets tend to penetrate further into the crossflow. (C) 2015 Elsevier Ltd. All rights reserved.

Effect of H-2/CO Composition on Extinction Strain Rates of Counterflow Syngas Flames

This work presents a detailed experimental and numerical investigation of the effect of H-2/CO composition on extinction characteristics of premixed and nonpremixed syngas flames. Experimental measurements of local and global extinction strain rates in counterflow diffusion flames have been reported at atmospheric pressure for six different compositions of syngas fuel. The concentration of H-2 was varied from 5 to 20% with a 3% increment, and correspondingly, CO was decreased from 35 to 20% in steps of 3%. Particle imaging velocimetry has been used to determine the local extinction strain rates. Local extinction strain rates increased with an increase in the H-2/CO ratio in both nonpremixed and premixed flames. The predicted extinction strain rates for both nonpremixed and premixed counterflow flames using five different mechanisms available in the literature were compared with measurements. The Davis H-2/CO and Ranzi H-2/CO mechanisms predicted extinction strain rates within 10% of experimental values irrespective of the H-2/CO ratio. In the nonpremixed case, the Cl mechanism by Li et al., GRI 3.0, and the Ranzi H-2/CO mechanism predicted extinction strain rates well for low H-2/CO ratios (from 5:35 to 14:26) but deviated from experiments for higher H-2/CO values (17:23 and 20:20). In addition to kinetics, preferential diffusion effects were found to affect the reaction zone significantly and create distinct localized reaction zone structures in nonpremixed flames, which could contribute to discrepancies in extinction predictions.

Multiobjective discrete particle swarm optimization for multisensor image alignment

A new technique is proposed for multisensor image registration by matching the features using discrete particle swarm optimization (DPSO). The feature points are first extracted from the reference and sensed image using improved Harris corner detector available in the literature. From the extracted corner points, DPSO finds the three corresponding points in the sensed and reference images using multiobjective optimization of distance and angle conditions through objective switching technique. By this, the global best matched points are obtained which are used to evaluate the affine transformation for the sensed image. The performance of the image registration is evaluated and concluded that the proposed approach is efficient.

Particle field imaging using digital in-line holography

Digital holography is the direct recording of holograms using a CCD camera and is an alternative to the use of a film or a plate. In this communication in-line digital holographic microscopy has been explored for its application in particle imaging in 3D. Holograms of particles of about 10 mu m size have been digitally reconstructed. Digital focusing was done to image the particles in different planes along the depth of focus. Digital holographic particle imaging results were compared with conventional optical microscope imaging. A methodology for dynamic analysis of microparticles in 3D using in-line digital holography has been proposed.