Comparison of particle dispersion obtained with simple stochastic models and les in pipe flow

A hybrid Eulerian-Lagrangian approach is employed to simulate heavy particle dispersion in turbulent pipe flow. The mean flow is provided by the Eulerian simulations developed by mean of JetCode, whereas the fluid fluctuations seen by particles are prescribed by a stochastic differential equation based on normalized Langevin. The statistics of particle velocity are compared to LES data which contain detailed statistics of velocity for particles with diameter equal to 20.4 µm. The model is in good agreement with the LES data for axial mean velocity whereas rms of axial and radial velocities should be adjusted.

Advances in global instability computations: from incompressible to hypersonic flow

Esta tesis constituye un gran avance en el conocimiento del estudio y análisis de inestabilidades hidrodinámicas desde un punto de vista físico y teórico, como consecuencia de haber desarrollado innovadoras técnicas para la resolución computacional eficiente y precisa de la parte principal del espectro correspondiente a los problemas de autovalores (EVP) multidimensionales que gobiernan la inestabilidad de flujos con dos o tres direcciones espaciales inhomogéneas, denominados problemas de estabilidad global lineal. En el contexto del trabajo de desarrollo de herramientas computacionales presentado en la tesis, la discretización mediante métodos de diferencias finitas estables de alto orden de los EVP bidimensionales y tridimensionales que se derivan de las ecuaciones de Navier-Stokes linealizadas sobre flujos con dos o tres direcciones espaciales inhomogéneas, ha permitido una aceleración de cuatro órdenes de magnitud en su resolución. Esta mejora de eficiencia numérica se ha conseguido gracias al hecho de que usando estos esquemas de diferencias finitas, técnicas eficientes de resolución de problemas lineales son utilizables, explotando el alto nivel de dispersión o alto número de elementos nulos en las matrices involucradas en los problemas tratados. Como más notable consecuencia cabe destacar que la resolución de EVPs multidimensionales de inestabilidad global, que hasta la fecha necesitaban de superordenadores, se ha podido realizar en ordenadores de sobremesa. Además de la solución de problemas de estabilidad global lineal, el mencionado desarrollo numérico facilitó la extensión de las ecuaciones de estabilidad parabolizadas (PSE) lineales y no lineales para analizar la inestabilidad de flujos que dependen fuertemente en dos direcciones espaciales y suavemente en la tercera con las ecuaciones de estabilidad parabolizadas tridimensionales (PSE-3D). Precisamente la capacidad de extensión del novedoso algoritmo PSE-3D para el estudio de interacciones no lineales de los modos de estabilidad, desarrollado íntegramente en esta tesis, permite la predicción de transición en flujos complejos de gran interés industrial y por lo tanto extiende el concepto clásico de PSE, el cuál ha sido empleado exitosamente durante las pasadas tres décadas en el mismo contexto para problemas de capa límite bidimensional. Típicos ejemplos de flujos incompresibles se han analizado en este trabajo sin la necesidad de recurrir a restrictivas presuposiciones usadas en el pasado. Se han estudiado problemas vorticales como es el caso de un vórtice aislado o sistemas de vórtices simulando la estela de alas, en los que la homogeneidad axial no se impone y así se puede considerar la difusión viscosa del flujo. Además, se ha estudiado el chorro giratorio turbulento, cuya inestabilidad se utiliza para mejorar las características de funcionamiento de combustores. En la tesis se abarcan adicionalmente problemas de flujos compresibles. Se presenta el estudio de inestabilidad de flujos de borde de ataque a diferentes velocidades de vuelo. También se analiza la estela formada por un elemento rugoso aislado en capa límite supersónica e hipersónica, mostrando excelentes comparaciones con resultados obtenidos mediante simulación numérica directa. Finalmente, nuevas inestabilidades se han identificado en el flujo hipersónico a Mach 7 alrededor de un cono elíptico que modela el vehículo de pruebas en vuelo HIFiRE-5. Los resultados comparan favorablemente con experimentos en vuelo, lo que subraya aún más el potencial de las metodologías de análisis de estabilidad desarrolladas en esta tesis. ABSTRACT The present thesis constitutes a step forward in advancing the frontiers of knowledge of fluid flow instability from a physical point of view, as a consequence of having been successful in developing groundbreaking methodologies for the efficient and accurate computation of the leading part of the spectrum pertinent to multi-dimensional eigenvalue problems (EVP) governing instability of flows with two or three inhomogeneous spatial directions. In the context of the numerical work presented in this thesis, the discretization of the spatial operator resulting from linearization of the Navier-Stokes equations around flows with two or three inhomogeneous spatial directions by variable-high-order stable finite-difference methods has permitted a speedup of four orders of magnitude in the solution of the corresponding two- and three-dimensional EVPs. This improvement of numerical performance has been achieved thanks to the high-sparsity level offered by the high-order finite-difference schemes employed for the discretization of the operators. This permitted use of efficient sparse linear algebra techniques without sacrificing accuracy and, consequently, solutions being obtained on typical workstations, as opposed to the previously employed supercomputers. Besides solution of the two- and three-dimensional EVPs of global linear instability, this development paved the way for the extension of the (linear and nonlinear) Parabolized Stability Equations (PSE) to analyze instability of flows which depend in a strongly-coupled inhomogeneous manner on two spatial directions and weakly on the third. Precisely the extensibility of the novel PSE-3D algorithm developed in the framework of the present thesis to study nonlinear flow instability permits transition prediction in flows of industrial interest, thus extending the classic PSE concept which has been successfully employed in the same context to boundary-layer type of flows over the last three decades. Typical examples of incompressible flows, the instability of which was analyzed in the present thesis without the need to resort to the restrictive assumptions used in the past, range from isolated vortices, and systems thereof, in which axial homogeneity is relaxed to consider viscous diffusion, as well as turbulent swirling jets, the instability of which is exploited in order to improve flame-holding properties of combustors. The instability of compressible subsonic and supersonic leading edge flows has been solved, and the wake of an isolated roughness element in a supersonic and hypersonic boundary-layer has also been analyzed with respect to its instability: excellent agreement with direct numerical simulation results has been obtained in all cases. Finally, instability analysis of Mach number 7 ow around an elliptic cone modeling the HIFiRE-5 flight test vehicle has unraveled flow instabilities near the minor-axis centerline, results comparing favorably with flight test predictions.

Spatial linear global instability analysis of the HIFiRE-5 elliptic cone model flow

The linear instability of the three-dimensional boundary-layer over the HIFiRE-5 flight test geometry, i.e. a rounded-tip 2:1 elliptic cone, at Mach 7, has been analyzed through spatial BiGlobal analysis, in a effort to understand transition and accurately predict local heat loads on next-generation ight vehicles. The results at an intermediate axial section of the cone, Re x = 8x10 5, show three different families of spatially amplied linear global modes, the attachment-line and cross- ow modes known from earlier analyses, and a new global mode, peaking in the vicinity of the minor axis of the cone, termed \center-line mode". We discover that a sequence of symmetric and anti-symmetric centerline modes exist and, for the basic ow at hand, are maximally amplied around F* = 130kHz. The wavenumbers and spatial distribution of amplitude functions of the centerline modes are documented

Wind tunnel analysis of the flow topology over complex topographies

El principal objetivo de este trabajo es aportar conocimiento para contestar la pregunta: ¿hasta que punto los ensayos en túnel aerodinámico pueden contribuir a determinar las características que afectan la respuesta dinámica de los aerogeneradores operando en terreno complejo?. Esta pregunta no es nueva, de hecho, el debate en la comunidad científica comenzó en el primer tercio del siglo pasado y aún está intensamente vivo. El método generalmente aceptado para enfrentar el mencionado problema consiste en analizar un caso de estudio determinado en el cual se aplican tanto ensayos a escala real como análisis computacionales y ensayos en túnel aerodinámico. Esto no es ni fácil ni barato. Esta es la razón por la cual desde el experimento de Askervein en 1988, los modelizadores del flujo atmosférico tuvieron que esperar hasta 2007 a que el experimento de Bolund fuese puesto en marcha con un despliegue de medios técnicos equivalentes (teniendo en cuenta la evolución de las tecnologías de sensores y computación). El problema contempla tantos aspectos que ambas experiencias fueron restringidas a condiciones de atmósfera neutra con efectos de Coriolis despreciables con objeto de reducir la complejidad. Este es el contexto en el que se ha desarrollado la presente tesis doctoral. La topología del flujo sobre la isla de Bolund ha sido estudiada mediante la reproducción del experimento de Bolund en los túneles aerodinámicos A9 y ACLA16 del IDR. Dos modelos de la isla de Bolund fueron fabricados a dos escalas, 1:230 y 1:115. El flujo de entrada en el túnel aerodinámico simulando la capa límite sin perturbar correspondía a régimen de transición (transitionally rough regime) y fue usado como situación de referencia. El modelo a escala 1:230 fue ensayado en el túnel A9 para determinar la presión sobre su superficie. La distribución del coeficiente de presión sobre la isla proporcionó una visualización y estimación de una región de desprendimiento sobre el pequeño acantilado situado al frente de la misma. Las medidas de presión instantánea con suficiente grado de resolución temporal pusieron de manifiesto la no estacionariedad en la región de desprendimiento. El modelo a escala 1:115 fue ensayado utilizando hilo caliente de tres componentes y un sistema de velocimetría por imágenes de partículas de dos componentes. El flujo fue caracterizado por el ratio de aceleración, el incremento normalizado de energía cinética turbulenta y los ángulos de inclinación y desviación horizontal. Los resultados a lo largo de la dirección 270°y alturas de 2 m y 5 m presentaron una gran similitud con los resultados a escala real del experimento de Bolund. Los perfiles verticales en las localizaciones de las torres meteorológicas mostraron un acuerdo significativo con los resultados a escala real. El análisis de los esfuerzos de Reynolds y el análisis espectral en las localizaciones de los mástiles meteorológicos presentaron niveles de acuerdo variados en ciertas posiciones, mientras que en otras presentaron claras diferencias. El mapeo horizontal del flujo, para una dirección de viento de 270°, permitió caracterizar el comportamiento de la burbuja intermitente de recirculación sobre el pequeño acantilado existente al frente de la isla así como de la región de relajación y de la capa de cortadura en la región corriente abajo de Bolund. Se realizaron medidas de velocidad con alta resolución espacial en planos perpendiculares a la dirección del flujo sin perturbar. Estas medidas permitieron detectar y caracterizar una estructura de flujo similar a un torbellino longitudinal con regiones con altos gradientes de velocidad y alta intensidad de turbulencia. Esta estructura de flujo es, sin duda, un reto para los modelos computacionales y puede considerarse un factor de riesgo para la operación de los aerogeneradores. Se obtuvieron y analizaron distribuciones espaciales de los esfuerzos de Reynolds mediante 3CHW y PIV. Este tipo de parámetros no constituyen parte de los resultados habituales en los ensayos en túnel sobre topografías y son muy útiles para los modelizadores que utilizan simulación de grades torbellinos (LES). Se proporciona una interpretación de los resultados obtenidos en el túnel aerodinámico en términos de utilidad para los diseñadores de parques eólicos. La evolución y variación de los parámetros del flujo a lo largo de líneas, planos y superficies han permitido identificar como estas propiedades del flujo podrían afectar la localización de los aerogeneradores y a la clasificación de emplazamientos. Los resultados presentados sugieren, bajo ciertas condiciones, la robustez de los ensayos en túnel para estudiar la topología sobre terreno complejo y su comparabilidad con otras técnicas de simulación, especialmente considerando el nivel de acuerdo del conjunto de resultados presentados con los resultados a escala real. De forma adicional, algunos de los parámetros del flujo obtenidos de las medidas en túnel son difícilmente determinables en ensayos a escala real o por medios computacionales, considerado el estado del arte. Este trabajo fue realizado como parte de las actividades subvencionadas por la Comisión Europea como dentro del proyecto FP7-PEOPLE-ITN-2008WAUDIT (Wind Resource Assessment Audit and Standardization) dentro de la FP7 Marie-Curie Initial Training Network y por el Ministerio Español de Economía y Competitividad dentro del proyecto ENE2012-36473, TURCO (Determinación en túnel aerodinámico de la distribución espacial de parámetros estadísticos de la turbulencia atmosférica sobre topografías complejas) del Plan Nacional de Investigación (Subprograma de investigación fundamental no orientada 2012). El informe se ha organizado en siete capítulos y un conjunto de anexos. En el primer capítulo se introduce el problema. En el capítulo dos se describen los medios experimentales utilizados. Seguidamente, en el capítulo tres, se analizan en detalle las condiciones de referencia del principal túnel aerodinámico utilizado en esta investigación. En el capítulo tres se presentan resultados de ensayos de presión superficial sobre un modelo de la isla. Los principales resultados del experimento de Bolund se reproducen en el capítulo cinco. En el capítulo seis se identifican diferentes estructuras del flujo sobre la isla y, finalmente, en el capitulo siete, se recogen las conclusiones y una propuesta de lineas de trabajo futuras. ABSTRACT The main objective of this work is to contribute to answer the question: to which extend can the wind tunnel testing contribute to determine the flow characteristics that affect the dynamic response of wind turbines operating in highly complex terrains?. This question is not new, indeed, the debate in the scientific community was opened in the first third of the past century and it is still intensely alive. The accepted approach to face this problem consists in analysing a given case study where full-scale tests, computational modelling and wind tunnel testing are applied to the same topography. This is neither easy nor cheap. This is is the reason why since the Askervein experience in 1988, the atmospheric flow modellers community had to wait till 2007 when the Bolund experiment was setup with a deployment of technical means equivalent (considering the evolution of the sensor and computing techniques). The problem is so manifold that both experiences were restricted to neutral conditions without Coriolis effects in order to reduce the complexity. This is the framework in which this PhD has been carried out. The flow topology over the Bolund Island has been studied by replicating the Bolund experiment in the IDR A9 and ACLA16 wind tunnels. Two mock-ups of the Bolund island were manufactured at two scales of 1:230 and 1:115. The in-flow in the empty wind tunnel simulating the incoming atmospheric boundary layer was in the transitionally rough regime and used as a reference case. The 1:230 model was tested in the A9 wind tunnel to measure surface pressure. The mapping of the pressure coefficient across the island gave a visualisation and estimation of a detachment region on the top of the escarpment in front of the island. Time resolved instantaneous pressure measurements illustrated the non-steadiness in the detachment region. The 1:115 model was tested using 3C hot-wires(HW) and 2C Particle Image Velocimetry(PIV). Measurements at met masts M3, M6, M7 and M8 and along Line 270°were taken to replicate the result of the Bolund experiment. The flow was characterised by the speed-up ratio, normalised increment of the turbulent kinetic energy, inclination angle and turning angle. Results along line 270°at heights of 2 m and 5 m compared very well with the full-scale results of the Bolund experiment. Vertical profiles at the met masts showed a significant agreement with the full-scale results. The analysis of the Reynolds stresses and the spectral analysis at the met mast locations gave a varied level of agreement at some locations while clear mismatch at others. The horizontal mapping of the flow field, for a 270°wind direction, allowed to characterise the behaviour of the intermittent recirculation bubble on top of the front escarpment followed by a relaxation region and the presence of a shear layer in the lee side of the island. Further detailed velocity measurements were taken at cross-flow planes over the island to study the flow structures on the island. A longitudinal vortex-like structure with high mean velocity gradients and high turbulent kinetic energy was characterised on the escarpment and evolving downstream. This flow structure is a challenge to the numerical models while posing a threat to wind farm designers when siting wind turbines. Spatial distribution of Reynold stresses were presented from 3C HW and PIV measurements. These values are not common results from usual wind tunnel measurements and very useful for modellers using large eddy simulation (LES). An interpretation of the wind tunnel results in terms of usefulness to wind farm designers is given. Evolution and variation of the flow parameters along measurement lines, planes and surfaces indicated how the flow field could affect wind turbine siting. Different flow properties were presented so compare the level of agreement to full-scale results and how this affected when characterising the site wind classes. The results presented suggest, under certain conditions, the robustness of the wind tunnel testing for studying flow topology over complex terrain and its capability to compare to other modelling techniques especially from the level of agreement between the different data sets presented. Additionally, some flow parameters obtained from wind tunnel measurements would have been quite difficult to be measured at full-scale or by computational means considering the state of the art. This work was carried out as a part of the activities supported by the EC as part of the FP7- PEOPLE-ITN-2008 WAUDIT project (Wind Resource Assessment Audit and Standardization) within the FP7 Marie-Curie Initial Training Network and by the Spanish Ministerio de Economía y Competitividad, within the framework of the ENE2012-36473, TURCO project (Determination of the Spatial Distribution of Statistic Parameters of Flow Turbulence over Complex Topographies in Wind Tunnel) belonging to the Spanish National Program of Research (Subprograma de investigación fundamental no orientada 2012). The report is organised in seven chapters and a collection of annexes. In chapter one, the problem is introduced. In chapter two the experimental setup is described. Following, in chapter three, the inflow conditions of the main wind tunnel used in this piece of research are analysed in detail. In chapter three, preliminary pressure tests results on a model of the island are presented. The main results from the Bolund experiment are replicated in chapter five. In chapter six, an identification of specific flow strutures over the island is presented and, finally, in chapter seven, conclusions and lines for future works related to the presented one are included.

Quasi-cylindrical approximation to the swirling flow in an atomizer chamber

A quasi-cylindrical approximation is used to analyse the axisymmetric swirling flow of a liquid with a hollow air core in the chamber of a pressure swirl atomizer. The liquid is injected into the chamber with an azimuthal velocity component through a number of slots at the periphery of one end of the chamber, and flows out as an anular sheet through a central orifice at the other end, following a conical convergence of the chamber wall. An effective inlet condition is used to model the effects of the slots and the boundary layer that develops at the nearby endwall of the chamber. An analysis is presented of the structure of the liquid sheet at the end of the exit orifice, where the flow becomes critical in the sense that upstream propagation of long-wave perturbations ceases to be possible. This nalysis leads to a boundary condition at the end of the orifice that is an extension of the condition of maximum flux used with irrotational models of the flow. As is well known, the radial pressure gradient induced by the swirling flow in the bulk of the chamber causes the overpressure that drives the liquid towards the exit orifice, and also leads to Ekman pumping in the boundary layers of reduced azimuthal velocity at the convergent wall of the chamber and at the wall opposite to the exit orifice. The numerical results confirm the important role played by the boundary layers. They make the thickness of the liquid sheet at the end of the orifice larger than predicted by rrotational models, and at the same time tend to decrease the overpressure required to pass a given flow rate through the chamber, because the large axial velocity in the boundary layers takes care of part of the flow rate. The thickness of the boundary layers increases when the atomizer constant (the inverse of a swirl number, proportional to the flow rate scaled with the radius of the exit orifice and the circulation around the air core) decreases. A minimum value of this parameter is found below which the layer of reduced azimuthal velocity around the air core prevents the pressure from increasing and steadily driving the flow through the exit orifice. The effects of other parameters not accounted for by irrotational models are also analysed in terms of their influence on the boundary layers.

The influence of the flow of the reacting gas on the conditions for a thermal explosion

The classical problem of thermal explosion is modified so that the chemically active gas is not at rest but is flowing in a long cylindrical pipe. Up to a certain section the heat-conducting walls of the pipe are held at low temperature so that the reaction rate is small and there is no heat release; at that section the ambient temperature is increased and an exothermic reaction begins. The question is whether a slow reaction regime will be established or a thermal explosion will occur. The mathematical formulation of the problem is presented. It is shown that when the pipe radius is larger than a critical value, the solution of the new problem exists only up to a certain distance along the axis. The critical radius is determined by conditions in a problem with a uniform axial temperature. The loss of existence is interpreted as a thermal explosion; the critical distance is the safe reactor’s length. Both laminar and developed turbulent flow regimes are considered. In a computational experiment the loss of the existence appears as a divergence of a numerical procedure; numerical calculations reveal asymptotic scaling laws with simple powers for the critical distance.

Wind Shear, Gust, and Yaw-Induced Dynamic Stall on Wind-Turbine Blades

This study examined the effect of a spanwise angle of attack gradient on the growth and stability of a dynamic stall vortex in a rotating system. It was found that a spanwise angle of attack gradient induces a corresponding spanwise vorticity gradient, which, in combination with spanwise flow, results in a redistribution of circulation along the blade. Specifically, when modelling the angle of attack gradient experienced by a wind turbine at the 30% span position during a gust event, the spanwise vorticity gradient was aligned such that circulation was transported from areas of high circulation to areas of low circulation, increasing the local dynamic stall vortex growth rate, which corresponds to an increase in the lift coefficient, and a decrease in the local vortex stability at this point. Reversing the relative alignment of the spanwise vorticity gradient and spanwise flow results in circulation transport from areas of low circulation generation to areas of high circulation generation, acting to reduce local circulation and stabilise the vortex. This circulation redistribution behaviour describes a mechanism by which the fluctuating loads on a wind turbine are magnified, which is detrimental to turbine lifetime and performance. Therefore, an understanding of this phenomenon has the potential to facilitate optimised wind turbine design.

Flow patterns in granulating systems

Particle flow patterns were investigated for wet granulation and dry powder mixing in ploughshare mixers using Positron Emission Particle Tracking (PEPT). In a 4-1 mixer, calcium carbonate with mean size 45 mum was granulated using a 50 wt.% solution of glycerol and water as binding fluid, and particle movement was followed using a 600-mum calcium hydroxy-phosphate tracer particle. In a 20-1 mixer, dry powder flow was studied using a 600-mum resin bead tracer particle to simulate the bulk polypropylene powder with mean size 600 mum. Important differences were seen between particle flow patterns for wet and dry systems. Particle speed relative to blade speed was lower in the wet system than in the dry system, with the ratios of average particle speed to blade tip speed for all experiments in the range 0.01-015. In the axial plane, the same particle motion was observed around each blade; this provides a significant advance for modelling flow in ploughshare mixers. For the future, a detailed understanding of the local velocity, acceleration and density variations around a plough blade will reveal the effects of flow patterns in granulating systems on the resultant distribution of granular product attributes such as size, density and strength. (C) 2002 Elsevier Science B.V All rights reserved.

Ocular blood flow measurements in healthy human myopic eyes

Background. To evaluate the haemodynamic features of young healthy myopes and emmetropes, in order to ascertain the perfusion profile of human myopia and its relationship with axial length prior to reaching a degenerative state. Methods The retrobulbar, microretinal and pulsatile ocular blood flow (POBF) of one eye of each of twenty-two high myopes (N=22, mean spherical equivalent (MSE) =-5.00D), low myopes (N=22, MSE-1.00 to-4.50D) and emmetropes (N=22, MSE±0.50D) was analyzed using color Doppler Imaging, Heidelberg retinal flowmetry and ocular blood flow analyser (OBF) respectively. Intraocular pressure, axial length (AL), systemic blood pressure, and body mass index were measured. Results. When compared to the emmetropes and low myopes, the AL was greater in high myopia (p<0.0001). High myopes showed higher central retinal artery resistance index (CRA RI) (p=0.004), higher peak systolic to end diastolic velocities ratio (CRA ratio) and lower end diastolic velocity (CRA EDv) compared to low myopes (p=0.014, p=0.037). Compared to emmetropes, high myopes showed lower OBFamplitude (OBFa) (p=0.016). The POBF correlated significantly with the systolic and diastolic blood velocities of the CRA (p=0.016, p=0.036). MSE and AL correlated negatively with OBFa (p=0.03, p=0.003), OBF volume (p=0.02, p<0.001), POBF (p=0.01, p<0.001) and positively with CRA RI (p=0.007, p=0.05). Conclusion. High myopes exhibited significantly reduced pulse amplitude and CRA blood velocity, the first of which may be due to an OBF measurement artefact or real decreased ocular blood flow pulsatility. Axial length and refractive error correlated moderately with the ocular pulse and with the resistance index of the CRA, which in turn correlated amongst themselves. It is hypothesized that the compromised pulsatile and CRA haemodynamics observed in young healthy myopes is an early feature of the decrease in ocular blood flow reported in pathological myopia. Such vascular features would increase the susceptibility for vascular and age-related eye diseases.

Variations of local heat transfer coefficient in piped flow of viscous liquids

Measurements were carried out to determine local coefficients of heat transfer in short lengths of horizontal pipe, and in the region of an discontinuity in pipe diameter. Laminar, transitional and turbulent flow regimes were investigated, and mixtures of propylene glycol and water were used in the experiments to give a range of viscous fluids. Theoretical and empirical analyses were implemented to find how the fundamental mechanism of forced convection was modified by the secondary effects of free convection, temperature dependent viscosity, and viscous dissipation. From experiments with the short tube it was possible to determine simple empirical relationships describing the axial distribution of the local 1usselt number and its dependence on the Reynolds and Prandtl numbers. Small corrections were made to account for the secondary effects mentioned above. Two different entrance configurations were investigated to demonstrate how conditions upstream could influence the heat transfer coefficients measured downstream In experiments with a sudden contraction in pipe diameter the distribution of local 1u3se1t number depended on the Prandtl number of the fluid in a complicated way. Graphical data is presented describing this dependence for a range of fluids indicating how the local Nusselt number varied with the diameter-ratio. Ratios up to 3.34:1 were considered. With a sudden divergence in pipe diameter, it was possible to derive the axial distribution of the local Nusse1t number for a range of Reynolds and Prandtl numbers in a similar way to the convergence experiments. Difficulty was encountered in explaining some of the measurements obtained at low Reynolds numbers, and flow visualization techniques wore used to determine the complex flow patterns which could lead to the anomalous results mentioned. Tests were carried out with divergences up to 1:3.34 to find the way in which the local Nusselt number varied with the diameter ratio, and a few experiments were carried out with very large ratios up .to 14.4. A limited amount of theoretical analysis of the 'divergence' system was carried out to substantiate certain explanations of the heat transfer mechanisms postulated.

Spray drier simulation and air flow pattern studies

The literature pertaining to the key stages of spray drying has been reviewed in the context of the mathematical modelling of drier performance. A critical review is also presented of previous spray drying models. A new mathematical model has been developed for prediction of spray drier performance. This is applicable to slurries of rigid, porous crust-forming materials to predict trajectories and drying profiles for droplets with a distribution of sizes sprayed from a centrifugal pressure nozzle. The model has been validated by comparing model predictions to experimental data from a pilot-scale counter-current drier and from a full-scale co-current drier. For the latter, the computed product moisture content was within 2%, and the computed air exit temperature within 10oC of experimental data. Air flow patterns have been investigated in a 1.2m diameter transparent countercurrent spray tower by flow visualisation. Smoke was introduced into various zones within the tower to trace the direction, and gauge the intensity, of the air flow. By means of a set of variable-angle air inlet nozzles, a variety of air entry configurations was investigated. The existence of a core of high rotational and axial velocity channelling up the axis of the tower was confirmed. The stability of flow within the core was found to be strongly dependent upon the air entry arrangement. A probe was developed for the measurement of air temperature and humidity profiles. This was employed for studying evaporation of pure water drops in a 1.2m diameter pilot-scale counter-current drier. A rapid approach to the exit air properties was detected within a 1m distance from the air entry ports. Measured radial profiles were found to be virtually flat but, from the axial profiles, the existence of plug-flow, well-mixed-flow and some degree of air short-circuiting can be inferred. The model and conclusions should assist in the improved design and optimum operation of industrial spray driers.

Flow patterns, performance and scale-up of distillation trays

This work is concerned with the nature of liquid flow across industrial sieve trays operating in the spray, mixed, and the emulsified flow regimes. In order to overcome the practical difficulties of removing many samples from a commercial tray, the mass transfer process was investigated in an air water simulator column by heat transfer analogy. The temperature of the warm water was measured by many thermocouples as the water flowed across the single pass 1.2 m diameter sieve tray. The thermocouples were linked to a mini computer for the storage of the data. The temperature data were then transferred to a main frame computer to generate temperature profiles - analogous to concentration profiles. A comprehensive study of the existing tray efficiency models was carried out using computerised numerical solutions. The calculated results were compared with experimental results published by the Fractionation Research Incorporation (FRl) and the existing models did not show any agreement with the experimental results. Only the Porter and Lockett model showed a reasonable agreement with the experimental results for cenain tray efficiency values. A rectangular active section tray was constructed and tested to establish the channelling effect and the result of its effect on circular tray designs. The developed flow patterns showed predominantly flat profiles and some indication of significant liquid flow through the central region of the tray. This comfirms that the rectangular tray configuration might not be a satisfactory solution for liquid maldistribution on sieve trays. For a typical industrial tray the flow of liquid as it crosses the tray from the inlet to the outlet weir could be affected by the mixing of liquid by the eddy, momentum and the weir shape in the axial or the transverse direction or both. Conventional U-shape profiles were developed when the operating conditions were such that the froth dispersion was in the mixed regime, with good liquid temperature distribution while in the spray regime. For the 12.5 mm hole diameter tray the constant temperature profiles were found to be in the axial direction while in the spray regime and in the transverse direction for the 4.5 mm hole tray. It was observed that the extent of the liquid stagnant zones at the sides of the tray depended on the tray hole diameter and was larger for the 4.5 mm hole tray. The liquid hold-up results show a high liquid hold-up at the areas of the tray with low liquid temperatures, this supports the doubts about the assumptions of constant point efficiency across an operating tray. Liquid flow over the outlet weir showed more liquid flow at the centre of the tray at high liquid loading with low liquid flow at both ends of the weir. The calculated results of the point and tray efficiency model showed a general increase in the calculated point and tray efficiencies with an increase in the weir loading, as the flow regime changed from the spray to the mixed regime the point and the tray efficiencies increased from approximately 30 to 80%.Through the mixed flow regime the efficiencies were found to remain fairly constant, and as the operating conditions were changed to maintain an emulsified flow regime there was a decrease in the resulting efficiencies. The results of the estimated coefficient of mixing for the small and large hole diameter trays show that the extent of liquid mixing on an operating tray generally increased with increasing capacity factor, but decreased with increasing weir loads. This demonstrates that above certain weir loads, the effect of eddy diffusion mechanism on the process of liquid mixing on an operating tray to be negligible.

The application of robotics to the turbine blade encapsulation process

A survey of the existing state-of-the-art of turbine blade manufacture highlights two operations that have not been automated namely that of loading of a turbine blade into an encapsulation die, and that of removing a machined blade from the encapsulation block. The automation of blade decapsulation has not been pursued. In order to develop a system to automate the loading of an encapsulation die a prototype mechanical handling robot has been designed together with a computer controlled encapsulation die. The robot has been designed as a mechanical handling robot of cylindrical geometry, suitable for use in a circular work cell. It is the prototype for a production model to be called `The Cybermate'. The prototype robot is mechanically complete but due to unforeseen circumstances the robot control system is not available (the development of the control system did not form a part of this project), hence it has not been possible to fully test and assess the robot mechanical design. Robot loading of the encapsulation die has thus been simulated. The research work with regard to the encapsulation die has focused on the development of computer controlled, hydraulically actuated, location pins. Such pins compensate for the inherent positional inaccuracy of the loading robot and reproduce the dexterity of the human operator. Each pin comprises a miniature hydraulic cylinder, controlled by a standard bidirectional flow control valve. The precision positional control is obtained through pulsing of the valves under software control, with positional feedback from an 8-bit transducer. A test-rig comprising one hydraulic location pin together with an opposing spring loaded pin has demonstrated that such a pin arrangement can be controlled with a repeatability of +/-.00045'. In addition this test-rig has demonstrated that such a pin arrangement can be used to gauge and compensate for the dimensional error of the component held between the pins, by offsetting the pin datum positions to allow for the component error. A gauging repeatability of +/- 0.00015' was demonstrated. This work has led to the design and manufacture of an encapsulation die comprising ten such pins and the associated computer software. All aspects of the control software except blade gauging and positional data storage have been demonstrated. Work is now required to achieve the accuracy of control demonstrated by the single pin test-rig, with each of the ten pins in the encapsulation die. This would allow trials of the complete loading cycle to take place.

Hydrodynamics and mass transfer studies in a scheibel extractor

A 10 cm diameter four-stage Scheibel column with dispersed phase wetted packing sections has been constructed to study the hydrodynamics and mass transfer using the system toluene-acetone-water. The literature pertaining to the above extractor has been examined and the important phenomena such as droplet break-up and coalescence, mass transfer and backmixing have been reviewed. A critical analysis of the backmixing or axial mixing models and the corresponding techniques for parameter estimation was applied and an optimization technique based on Marquardt's algorithm was implemented. A single phase sampling technique was developed to estimate the acetone concentration profile in both phases along the column. Column flooding characteristics were investigated under various operating conditions and it was found that, when the impellers were located at about DI/5cm from the upper surface of the pads, the limiting flow rates increased with impeller speed. This unusual behaviour was explained in terms of the pumping effect created by the turbine impellers. Correlations were developed to predict Sauter mean drop diameters. A five-cell with backflow model was used to estimate the column performance (stage efficiency) and phases non-ideality (backflow parameters). Overall mass transfer coefficients were computed using the above model and compared with those calculated using the correlations based on single drop mechanism.

Hydrodynamic phenomena in the downcomer during flow rate transients in the primary circuit of a PWR

Mixing phenomena observed when the flow rate in a single loop of the primary circuit is changed can influence the operation of pressurized water reactor (PWR) by inducing local gradients of boron concentration or coolant temperature. Analysis of one-dimensional Laser Doppler Anemometry (LDA) measurements during the start-up and shutdown of pump on a single loop of the ROCOM test facility has been performed. The effect of a step change and a ramped change in the flow rate on the axial and azimuthal velocities was examined. Numerical simulations were also performed for the step change in the flow rate that gave quantitative agreement with the axial velocities. Phenomenological agreement was made on the turbulent kinetic energy; however, observed values were a factor of 2.5 less than the turbulent kinetic energy derived from the measurements. © 2007.