Tunable hydrodynamic chromatography of microparticles localized in short microchannels.

This paper describes a new way to perform hydrodynamic chromatography (HDC) for the size separation of particles based on a unique recirculating flow pattern. Pressure-driven (PF) and electro-osmotic flows (EOF) are opposed in narrow glass microchannels that expand at both ends. The resulting bidirectional flow turns into recirculating flow because of nonuniform microchannel dimensions. This hydrodynamic effect, combined with the electrokinetic migration of the particles themselves, results in a trapping phenomenon, which we have termed flow-induced electrokinetic trapping (FIET). In this paper, we exploit recirculating flow and FIET to perform a size-based separation of samples of microparticles trapped in a short separation channel using a HDC approach. Because these particles have the same charge (same zeta potential), they exhibit the same electrophoretic mobility, but they can be separated according to size in the recirculating flow. While trapped, particles have a net drift velocity toward the low-pressure end of the channel. When, because of a change in the externally applied PF or electric field, the sign of the net drift velocity reverses, particles can escape the separation channel in the direction of EOF. Larger particles exhibit a larger net drift velocity opposing EOF, so that the smaller particles escape the separation channel first. In the example presented here, a sample plug containing 2.33 and 2.82 microm polymer particles was introduced from the inlet into a 3-mm-long separation channel and trapped. Through tuning of the electric field with respect to the applied PF, the particles could be separated, with the advantage that larger particles remained trapped. The separation of particles with less than 500 nm differences in diameter was performed with an analytical resolution comparable to that of baseline separation in chromatography. When the sample was not trapped in the separation channel but located further downstream, separations could be carried out continuously rather than in batch. Smaller particles could successfully pass through the separation channel, and particles were separated by size. One of the main advantages of exploiting FIET for HDC is that this method can be applied in quite short (a few millimeters) channel geometries. This is in great contrast to examples published to date for the separation of nanoparticles in much longer micro- and nanochannels.

Turbomachinery computational fluid dynamics: asymptotes and paradigm shifts.

This paper reviews the development of computational fluid dynamics (CFD) specifically for turbomachinery simulations and with a particular focus on application to problems with complex geometry. The review is structured by considering this development as a series of paradigm shifts, followed by asymptotes. The original S1-S2 blade-blade-throughflow model is briefly described, followed by the development of two-dimensional then three-dimensional blade-blade analysis. This in turn evolved from inviscid to viscous analysis and then from steady to unsteady flow simulations. This development trajectory led over a surprisingly small number of years to an accepted approach-a 'CFD orthodoxy'. A very important current area of intense interest and activity in turbomachinery simulation is in accounting for real geometry effects, not just in the secondary air and turbine cooling systems but also associated with the primary path. The requirements here are threefold: capturing and representing these geometries in a computer model; making rapid design changes to these complex geometries; and managing the very large associated computational models on PC clusters. Accordingly, the challenges in the application of the current CFD orthodoxy to complex geometries are described in some detail. The main aim of this paper is to argue that the current CFD orthodoxy is on a new asymptote and is not in fact suited for application to complex geometries and that a paradigm shift must be sought. In particular, the new paradigm must be geometry centric and inherently parallel without serial bottlenecks. The main contribution of this paper is to describe such a potential paradigm shift, inspired by the animation industry, based on a fundamental shift in perspective from explicit to implicit geometry and then illustrate this with a number of applications to turbomachinery.

Strong plasmonic enhancement of photovoltage in graphene.

From the wide spectrum of potential applications of graphene, ranging from transistors and chemical sensors to nanoelectromechanical devices and composites, the field of photonics and optoelectronics is believed to be one of the most promising. Indeed, graphene's suitability for high-speed photodetection was demonstrated in an optical communication link operating at 10 Gbit s(-1). However, the low responsivity of graphene-based photodetectors compared with traditional III-V-based ones is a potential drawback. Here we show that, by combining graphene with plasmonic nanostructures, the efficiency of graphene-based photodetectors can be increased by up to 20 times, because of efficient field concentration in the area of a p-n junction. Additionally, wavelength and polarization selectivity can be achieved by employing nanostructures of different geometries.

NUMERICAL STUDY OF THE 3D FLOWFIELD IN A TRANSONIC COMPRESSOR ROTOR WITH A MODELLING OF THE TIP CLEARANCE FLOW.

This paper describes a computer code aimed at solving the equations of three dimensional viscous compressible flow in turbomachine geometries. The code is applied to the study of the flowfield in a transonic axial compressor rotor at design speed at both maximum flow and towards stall. The predicted flowfield is compared with the laser measurements and the performance of the code discussed. In addition the discussion highlights the change in the predicted endwall and tip clearance flows as the rotor operating point is moved towards stall.

Effect of blade tip geometry on the tip leakage flow in axial turbine cascades

The phenomenon of tip leakage has been studied in two linear cascades of turbine blades. The investigation includes an examination of the performance of the cascades with a variety of tip geometries. The effects of using plain tips, suction side squealers, and pressure side squealers are reported. Traverses of the exit flow field were made in order to determine the overall performance. A method of calculating the tip discharge coefficients for squealer geometries is put forward. In linking the tip discharge coefficient and cascade losses, a procedure for predicting the relative performance of tip geometries is developed. The model is used to examine the results obtained using the different tip treatments and to highlight the important aspects of the loss generation process.

The effect of surface hardening on the elastic shakedown of elliptical contacts

The effect of varying both the aspect ratio and the coefficient of friction of contacts with elliptical geometry on their elastic shakedown performance has been examined theoretically for surfaces with two types of subsurface hardness or strength profiles. In stepwise hardening the hard layer is of uniform strength while in linear hardening its strength reduces from a maximum at the surface to that of the core at the base of the hardened layer. The shakedown load is expressed as the ratio of the maximum Hertzian pressure to the strength of the core material. As the depth of hardening, expressed as a multiple of the elliptical semi-axis, is increased so the potential shakedown load increases from a level that is appropriate to a uniform half-space of unhardened material to a value reflecting the hardness of the surface and near-surface material. In a step-hardened material, the shakedown limit for a surface 'pummelled' by the passage of a sequence of such loads reaches a cut-off or plateau value, which cannot be exceeded by further increases in hardening depth irrespective of the value of the friction coefficient. For a linear-hardened material the corresponding plateau is approached asymptotically. The work confirms earlier results on the upper bounds on shakedown of both point and line contacts and provides numerical values of shakedown loads for intermediate geometries. In general, the case depth required to achieve a given shakedown limit reduces in moving from a transversely moving nominal line load to an axisymmetric point load.

Application of a laminar flamelet model to confined explosion hazards

The majority of computational studies of confined explosion hazards apply simple and inaccurate combustion models, requiring ad hoc corrections to obtain realistic flame shapes and often predicting an order of magnitude error in the overpressures. This work describes the application of a laminar flamelet model to a series of two-dimensional test cases. The model is computationally efficient applying an algebraic expression to calculate the flame surface area, an empirical correlation for the laminar flame speed and a novel unstructured, solution adaptive numerical grid system which allows important features of the solution to be resolved close to the flame. Accurate flame shapes are predicted, the correct burning rate is predicted near the walls, and an improvement in the predicted overpressures is obtained. However, in these fully turbulent calculations the overpressures are still too high and the flame arrival times too low, indicating the need for a model for the early laminar burning phase. Due to the computational expense, it is unrealistic to model a laminar flame in the complex geometries involved and therefore a pragmatic approach is employed which constrains the flame to propagate at the laminar flame speed. Transition to turbulent burning occurs at a specified turbulent Reynolds number. With the laminar phase model included, the predicted flame arrival times increase significantly, but are still too low. However, this has no significant effect on the overpressures, which are predicted accurately for a baffled channel test case where rapid transition occurs once the flame reaches the first pair of baffles. In a channel with obstacles on the centreline, transition is more gradual and the accuracy of the predicted overpressures is reduced. However, although the accuracy is still less than desirable in some cases, it is much better than the order of magnitude error previously expected.

Wake, shock and potential field interactions in a 1.5 stage turbine: Part I: Vane-rotor and rotor-vane interaction

The composition of the time-resolved surface pressure field around a high-pressure rotor blade caused by the presence of neighboring blade rows was studied, with the individual effects of wake, shock and potential field interaction being determined. Two test geometries were considered: first, a high-pressure turbine stage coupled with a swan-necked diffuser exit duct; secondly, the same high-pressure stage but with a vane located in the downstream duct. Both tests were carried out at engine-representative Mach and Reynolds numbers. By comparing the results to time-resolved computational predictions of the flowfield, the accuracy with which the computation predicts blade interaction was determined. It was found that in addition to upstream vane-rotor and rotor-downstream vane interactions, a new interaction mechanism was found resulting from the interaction between the downstream vane's potential field and the upstream vane's trailing edge potential field and shock.

Wake, shock and potential field interactions in a 1.5 stage turbine: Part II: Vane-vane interaction and discussion of results

The composition of the time-resolved surface pressure field around a high-pressure rotor blade caused by the presence of neighboring blade rows was studied, with the individual effects of wake, shock and potential field interaction being determined. Two test geometries were considered: first, a high-pressure turbine stage coupled with a swan-necked diffuser exit duct; secondly, the same high-pressure stage but with a vane located in the downstream duct. Both tests were carried out at engine-representative Mach and Reynolds numbers. By comparing the results to time-resolved computational predictions of the flowfield, the accuracy with which the computation predicts blade interaction was determined. Evidence was obtained that for a large downstream vane, the flow conditions in the rotor passage, at any instant in time, are close to being steady state.

Heat transfer and aerodynamics of turbine blade tips in a linear cascade

Local measurements of the heat transfer coefficient and pressure coefficient were conducted on the tip and near tip region of a generic turbine blade in a five-blade linear cascade. Two tip clearance gaps were used: 1.6% and 2.8% chord. Data was obtained at a Reynolds number of 2.3 × 10 5 based on exit velocity and chord. Three different tip geometries were investigated: a flat (plain) tip, a suction-side squealer, and a cavity squealer. The experiments reveal that the flow through the plain gap is dominated by flow separation at the pressure-side edge and that the highest levels of heat transfer are located where the flow reattaches on the tip surface. High heat transfer is also measured at locations where the tip-leakage vortex has impinged onto the suction surface of the aerofoil. The experiments are supported by flow visualisation computed using the CFX CFD code which has provided insight into the fluid dynamics within the gap. The suction-side and cavity squealers are shown to reduce the heat transfer in the gap but high levels of heat transfer are associated with locations of impingement, identified using the flow visualisation and aerodynamic data. Film cooling is introduced on the plain tip at locations near the pressure-side edge within the separated region and a net heat flux reduction analysis is used to quantify the performance of the successful cooling design. copyright © 2005 by ASME.

Viscous layer meshes from level sets on cartesian meshes

The long term goal of our work is to enable rapid prototyping design optimization to take place on geometries of arbitrary size in a spirit of a real time computer game. In recent papers we have reported the integration of a Level Set based geometry kernel with an octree-based cut-Cartesian mesh generator, RANS flow solver and post-processing all within a single piece of software - and all implemented in parallel with commodity PC clusters as the target. This work has shown that it is possible to eliminate all serial bottlenecks from the CED Process. This paper reports further progress towards our goal; in particular we report on the generation of viscous layer meshes to bridge the body to the flow across the cut-cells. The Level Set formulation, which underpins the geometry representation, is used as a natural mechanism to allow rapid construction of conformal layer meshes. The guiding principle is to construct the mesh which most closely approximates the body but remains solvable. This apparently novel approach is described and examples given.

Towards jet flow LES of conceptual nozzles for acoustic predictions

Progress in simulating chevron nozzle jet flows using ILES/RANS-ILES approaches and using the Ffowcs Williams and Hawkings (FW-H) surface integral method to predict the radiated far field sound is presented in this paper. With the focus on the realistic chevron geometries, SMC001 and SMC006, coarse and fine meshes are generated in the range of 3∼13 million mesh cells. Throughout this work, to minimize numerical dissipation introduced by mesh quality issues, the hexahedral cell type is used. Numerical simulations are then carried out with cell-vertex and cell-centered codes. Despite the modest grids, mean velocities and turbulent statistics are found to be in reasonable accord with measurements. Also, far field sound levels predicted by the FW-H post processor are encouraging. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc.

Flame imaging of gas-turbine relight

High-altitude relight inside a lean-direct-injection gas-turbine combustor is investigated experimentally by highspeed imaging. Realistic operating conditions are simulated in a ground-based test facility, with two conditions being studied: one inside and one outside the combustor ignition loop. The motion of hot gases during the early stages of relight is recorded using a high-speed camera. An algorithm is developed to track the flame movement and breakup, revealing important characteristics of the flame development process, including stabilization timescales, spatial trajectories, and typical velocities of hot gas motion. Although the observed patterns of ignition failure are in broad agreement with results from laboratory-scale studies, other aspects of relight behavior are not reproduced in laboratory experiments employing simplified flow geometries and operating conditions. For example, when the spark discharge occurs, the air velocity below the igniter in a real combustor is much less strongly correlated to ignition outcome than laboratory studies would suggest. Nevertheless, later flame development and stabilization are largely controlled by the cold flowfield, implying that the location of the igniter may, in the first instance, be selected based on the combustor cold flow. Copyright © 2010.

Investigation of the effect of combustor cooling geometry on acoustic energy absorption

To control combustion instabilities occurring in LPP gas turbine combustors, several active and passive systems have been developed in recent years. The combustion chamber cooling geometry has the potential to influence instability feedback loops by absorbing acoustical energy inside the combustor. The design of the cooling liner and the geometry of the cooling plenum and the cooling air flow rate have a significant influence on the absorption characteristics of the system. This paper presents the results of a cold flow study which was carried out in the course of a comprehensive study on the influence of the cooling geometry on combustor thermoacoustics. Absorption characteristics of three different cooling liner geometries and non-perforated plates were determined over a frequency range from 50 Hz to 600 Hz for different cooling flow rates and different cooling plenum volumes. The experimental results compared well with results from a low order thermoacoustic network model. The acoustic energy absorption spectrum of a cooling liner with 90°-hole configuration was found to be strongly dependent on cooling flow rate and cooling plenum volume, whereas the absorption spectrum of cooling liners with 25°-holes were found to be strongly dependent on the cooling plenum volume, but less dependent on the cooling air flow rate. All cooling liner setups with perforations were capable of increased acoustic absorption over a broad band of frequencies compared to the case of non-perforated combustor walls. © 2010 by Johannes Schmidt.

Fluid-structure interaction simulation of an inflatable aerodynamic tension-cone decelerator

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. In this paper, we focus on an inflatable tension cone design that has potential advantages over other geometries. A computational fluid-structure interaction model of a tension cone is employed to investigate the behavior of the inflatable aeroshell at supersonic speeds for conditions matching recent experimental results. A parametric study is carried out to investigate the deflections of the tension cone as a function of inflation pressure of the torus at a Mach of 2.5. Comparison of the behavior of the structure, amplitude of deformations, and determined loads are reported. © 2010 by the American Institute of Aeronautics and Astronautics, Inc.