Strong Down-Valley Low-Level Jets over the Atacama Desert: Observational Characterization

The near-surface wind and temperature regime at three points in the Atacama Desert of northern Chile is described using two-year multi-level measurements from 80-m towers located in an altitude range between 2100 and 2700 m ASL. The data reveal the frequent development of strong nocturnal drainage flows at all sites. Down-valley nose-shaped wind speed profiles are observed with maximum values occurring at heights between 20 m and 60 m AGL. The flow intensity shows considerable inter-daily variability and a seasonal modulation of maximum speeds, which in the cold season can attain hourly average values larger than 20 m s−1. Turbulent mixing appears significant over the full tower layer, affecting the curvature of the nighttime temperature profile and possibly explaining the observed increase of surface temperatures in the down-valley direction. Nocturnal valley winds and temperatures are weakly controlled by upper-air conditions observed at the nearest aerological station. Estimates of terms in the momentum budget for the development and the quasi-stationary phases of the down-valley flows suggest that the pressure gradient force due to the near-surface cooling along the sloping valley axes plays an important role in these drainage flows. A scale for the jet nose height of equilibrium turbulent down-slope jets is proposed, based on surface friction velocity and surface inversion intensity. At one of the sites this scale explains about 70% of the case-to-case observed variance of jet nose heights. Further modeling and observational work is needed, however, in order to better define the dynamics, extent and turbulence structure of this flow system, which has significant wind-energy, climatic and environmental implications.

Effect of meal ingestion on liver stiffness in patients with cirrhosis and portal hypertension.

BACKGROUND AND AIMS Liver stiffness is increasingly used in the non-invasive evaluation of chronic liver diseases. Liver stiffness correlates with hepatic venous pressure gradient (HVPG) in patients with cirrhosis and holds prognostic value in this population. Hence, accuracy in its measurement is needed. Several factors independent of fibrosis influence liver stiffness, but there is insufficient information on whether meal ingestion modifies liver stiffness in cirrhosis. We investigated the changes in liver stiffness occurring after the ingestion of a liquid standard test meal in this population. METHODS In 19 patients with cirrhosis and esophageal varices (9 alcoholic, 9 HCV-related, 1 NASH; Child score 6.9±1.8), liver stiffness (transient elastography), portal blood flow (PBF) and hepatic artery blood flow (HABF) (Doppler-Ultrasound) were measured before and 30 minutes after receiving a standard mixed liquid meal. In 10 the HVPG changes were also measured. RESULTS Post-prandial hyperemia was accompanied by a marked increase in liver stiffness (+27±33%; p<0.0001). Changes in liver stiffness did not correlate with PBF changes, but directly correlated with HABF changes (r = 0.658; p = 0.002). After the meal, those patients showing a decrease in HABF (n = 13) had a less marked increase of liver stiffness as compared to patients in whom HABF increased (n = 6; +12±21% vs. +62±29%,p<0.0001). As expected, post-prandial hyperemia was associated with an increase in HVPG (n = 10; +26±13%, p = 0.003), but changes in liver stiffness did not correlate with HVPG changes. CONCLUSIONS Liver stiffness increases markedly after a liquid test meal in patients with cirrhosis, suggesting that its measurement should be performed in standardized fasting conditions. The hepatic artery buffer response appears an important factor modulating postprandial changes of liver stiffness. The post-prandial increase in HVPG cannot be predicted by changes in liver stiffness.

Unusual Southern Hemisphere tree growth patterns induced by changes in the Southern Annular Mode

Recent changes in the summer climate of the Southern Hemisphere extra-tropics are primarily related to the dominance of the positive phase of the Southern Annular Mode1, 2. This shift in the behaviour of the Southern Annular Mode—essentially a measure of the pressure gradient between Southern Hemisphere mid and high latitudes—has been predominantly induced by polar stratospheric ozone depletion2, 3, 4. The concomitant southward expansion of the dry subtropical belts5, 6 could have consequences for forest growth. Here, we use tree-ring records from over 3,000 trees in South America, Tasmania and New Zealand to identify dominant patterns of tree growth in recent centuries. We show that the foremost patterns of growth between 1950 and 2000 differed significantly from those in the previous 250 years. Specifically, growth was higher than the long-term average in the subalpine forests of Tasmania and New Zealand, but lower in the dry-mesic forests of Patagonia. We further demonstrate that variations in the Southern Annular Mode can explain 12–48% of the tree growth anomalies in the latter half of the twentieth century. Tree-ring-based reconstructions of summer Southern Annular Mode indices suggest that the high frequency of the positive phase since the 1950s is unprecedented in the past 600 years. We propose that changes in the Southern Annular Mode have significantly altered tree growth patterns in the Southern Hemisphere.

New Self-Expanding Transcatheter Valve for Off-Pump Transatrial Mitral Valve-In-Ring Implantation.

OBJECTIVES To validate a self-expanding transcatheter valve for off-pump transatrial mitral valve-in-ring (VIR) implantation via a left thoracotomy. METHODS Mitral valve annuloplasty was performed via sternotomy during cardiopulmonary bypass on 9 pigs. After successful weaning from extracorporal circulation, the custom-made, self-expanding transcatheter VIR device was deployed under fluoroscopic guidance within the annuloplasty ring via a left thoracotomy. Hemodynamic data before and after the implantation were recorded. Mitral annulus diameter and valve area were measured by echocardiography. Transvalvular and left-ventricular outflow-tract pressure gradient were measured invasively. RESULTS Eight successful implantations were performed. Implantation failed in 1 pig because of difficulty with technical delivery of the sheath. Mean transatrial procedure time was 12.6 ± 1.7 min. Hemodynamic status during transatrial implantation was stable, and differences were not statistically significant. Mean mitral annulus diameter and mean mitral orifice area were 2.32 ± 0.2 and 3.84 ± 0.55 cm2, respectively. Mild regurgitation was detected in 7 animals and moderate regurgitation in 1. Mean gradients were 6.1 ± 5.0 mm Hg across the device. Postmortem examination confirmed adequate positioning of devices within the annuloplasty ring. CONCLUSIONS This custom-made transcatheter device allows for safe and reproducible off-pump transatrial mitral VIR implantations. Transatrial access is a promising route to facilitate VIR implantations. Our custom-made stent-valve may be suitable for VIR procedures.

Exaggerated pulmonary hypertension and right ventricular dysfunction in high-altitude dwellers with patent foramen ovale.

BACKGROUND There is considerable interindividual variability in pulmonary artery pressure among high-altitude (HA) dwellers, but the underlying mechanism is not known. At low altitude, a patent foramen ovale (PFO) is present in about 25% of the general population. Its prevalence is increased in clinical conditions associated with pulmonary hypertension and arterial hypoxemia, and it is thought to aggravate these problems. METHODS We searched for a PFO (transesophageal echocardiography) in healthy HA dwellers (n = 22) and patients with chronic mountain sickness (n = 35) at 3,600 m above sea level and studied its effects (transthoracic echocardiography) on right ventricular (RV) function, pulmonary artery pressure, and vascular resistance at rest and during mild exercise (50 W), an intervention designed to further increase pulmonary artery pressure. RESULTS The prevalence of PFO (32%) was similar to that reported in low-altitude populations and was not different in participants with and without chronic mountain sickness. Its presence was associated with RV enlargement at rest and an exaggerated increase in right-ventricular-to-right-atrial pressure gradient (25 ± 7 mm Hg vs 15 ± 9 mm Hg, P < .001) and a blunted increase in fractional area change of the right ventricle (3% [-1%, 5%] vs 7% [3%, 16%], P = .008) during mild exercise. CONCLUSIONS These findings show, we believe for the first time, that although the prevalence of PFO is not increased in HA dwellers, its presence appears to facilitate pulmonary vasoconstriction and RV dysfunction during a mild physical effort frequently associated with daily activity. TRIAL REGISTRY ClinicalTrials.gov; No.: NCT01182792; URL: www.clinicaltrials.gov.

Serum tests, liver stiffness and artificial neural networks for diagnosing cirrhosis and portal hypertension

BACKGROUND The diagnostic performance of biochemical scores and artificial neural network models for portal hypertension and cirrhosis is not well established. AIMS To assess diagnostic accuracy of six serum scores, artificial neural networks and liver stiffness measured by transient elastography, for diagnosing cirrhosis, clinically significant portal hypertension and oesophageal varices. METHODS 202 consecutive compensated patients requiring liver biopsy and hepatic venous pressure gradient measurement were included. Several serum tests (alone and combined into scores) and liver stiffness were measured. Artificial neural networks containing or not liver stiffness as input variable were also created. RESULTS The best non-invasive method for diagnosing cirrhosis, portal hypertension and oesophageal varices was liver stiffness (C-statistics=0.93, 0.94, and 0.90, respectively). Among serum tests/scores the best for diagnosing cirrhosis and portal hypertension and oesophageal varices were, respectively, Fibrosis-4, and Lok score. Artificial neural networks including liver stiffness had high diagnostic performance for cirrhosis, portal hypertension and oesophageal varices (accuracy>80%), but were not statistically superior to liver stiffness alone. CONCLUSIONS Liver stiffness was the best non-invasive method to assess the presence of cirrhosis, portal hypertension and oesophageal varices. The use of artificial neural networks integrating different non-invasive tests did not increase the diagnostic accuracy of liver stiffness alone.

Prevention of Rebleeding From Esophageal Varices in Patients With Cirrhosis Receiving Small-Diameter Stents Versus Hemodynamically Controlled Medical Therapy

BACKGROUND & AIMS Patients with cirrhosis and variceal hemorrhage have a high risk of rebleeding. We performed a prospective randomized trial to compare the prevention of rebleeding in patients given a small-diameter covered stent vs those given hepatic venous pressure gradient (HVPG)-based medical therapy prophylaxis. METHODS We performed an open-label study of patients with cirrhosis (92% Child class A or B, 70% alcoholic) treated at 10 medical centers in Germany. Patients were assigned randomly more than 5 days after variceal hemorrhage to groups given a small covered transjugular intrahepatic portosystemic stent-shunt (TIPS) (8 mm; n = 90), or medical reduction of portal pressure (propranolol and isosorbide-5-mononitrate; n = 95). HVPG was determined at the time patients were assigned to groups (baseline) and 2 weeks later. In the medical group, patients with an adequate reduction in HVPG (responders) remained on the drugs whereas nonresponders underwent only variceal band ligation. The study was closed 10 months after the last patient was assigned to a group. The primary end point was variceal rebleeding. Survival, safety (adverse events), and quality of life (based on the Short Form-36 health survey) were secondary outcome measures. RESULTS A significantly smaller proportion of patients in the TIPS group had rebleeding within 2 years (7%) than in the medical group (26%) (P = .002). A slightly higher proportion of patients in the TIPS group experienced adverse events, including encephalopathy (18% vs 8% for medical treatment; P = .05). Rebleeding occurred in 6 of 23 patients (26%) receiving medical treatment before hemodynamic control was possible. Per-protocol analysis showed that rebleeding occurred in a smaller proportion of the 32 responders (18%) than in nonresponders who received variceal band ligation (31%) (P = .06). Fifteen patients from the medical group (16%) underwent TIPS placement during follow-up evaluation, mainly for refractory ascites. Survival time and quality of life did not differ between both randomized groups. CONCLUSIONS Placement of a small-diameter, covered TIPS was straightforward and prevented variceal rebleeding in patients with Child A or B cirrhosis more effectively than drugs, which often required step-by-step therapy. However, TIPS did not increase survival time or quality of life and produced slightly more adverse events. Clinical Trial no: ISRCTN 16334693.

Eolian records from sediment core RC27-61 in the Arabian Sea

The modern Indian Ocean summer monsoon is driven by differential heating between the Asian continent and the Indian Ocean to the south. This differential heating produces a strong pressure gradient which drives southwest monsoon winds during June, July, and August. Satellite and meteorological observations, aerosol measurements, sediment trap studies, and mineralogical studies indicate an atmospheric mode of transport for modern lithogenic sediments in the northwest Arabian Sea. Analyses of lithogenic grain size and mass accumulation rate (MAR) records from the Owen Ridge indicate that eolian transport has been the primary mode of transport for the past 370 kyr. Visual inspection shows that the MAR record is positively correlated with global ice volume as indicated by the marine delta18O record. In contrast, the grain-size record varies at a much higher frequency, showing little correlation to either the MAR or the delta18O records. Spectral analyses confirm these relationships, indicating that the lithogenic grain-size and MAR records are coherent only over the precession band whereby the grain size leads the MAR by 124° (~8 kyr). We conclude that an eolian transport mechanism is the only mechanism that allows for this phase difference and at the same time is supported by comparison of the grain size and MAR with independent eolian records. We use lithogenic grain size as a paleoclimatic indicator of summer monsoon wind strength and lithogenic MAR as a paleoclimatic indicator of source-area aridity. These interpretations are supported by comparison of the lithogenic records to independent indicators of wind strength (Globigerina bulloides upwelling record) and aridity (a loess record from central China). Such comparisons indicate high coherence and zero phase relationships. Our work supports the findings of previous studies which have documented the link between monsoon strength and the Earth's axial precession cycles. Both the lithogenic MAR and the grain-size records have high coherency with precessional insolation. Maximum lithogenic MAR (source-area aridity) is in phase with delta18O (global ice volume) and leads maximum precessional insolation by 88° (~6 kyr). We attribute this lead to the influence of glacial conditions on the aridity, and therefore the deflation potential, of the source areas. Maximum lithogenic grain size (summer monsoon wind strength) lags maximum precession by 148° (~9 kyr). We attribute this lag both to the influence of global and/or local ice volume and to the availability of latent heat from the southern hemisphere Indian Ocean, the two of which combine to determine the strength of the Indian Ocean monsoon.

Methylation index of Branched Tetraethers (MBT) and Cyclization ratio of Branched Tetraethers (CBT), annual mean air temperature and soil pH of sediment core GeoB6518-1 in the Congo River basin

We analyzed the distribution of branched tetraether membrane lipids derived from soil bacteria in a marine sediment record that was recovered close to the Congo River outflow, and the results enabled us to reconstruct large-scale continental temperature changes in tropical Africa that span the past 25,000 years. Tropical African temperatures gradually increased from ~21° to 25°C over the last deglaciation, which is a larger warming than estimated for the tropical Atlantic Ocean. A direct comparison with sea-surface temperature estimates from the same core revealed that the land-sea temperature difference was, through the thermal pressure gradient, an important control on central African precipitation patterns.

Experimental studies on transitional separated boundary layers

Separated transitional boundary layers appear on key aeronautical processes such as the flow around wings or turbomachinery blades. The aim of this thesis is the study of these flows in representative scenarios of technological applications, gaining knowledge about phenomenology and physical processes that occur there and, developing a simple model for scaling them. To achieve this goal, experimental measurements have been carried out in a low speed facility, ensuring the flow homogeneity and a low disturbances level such that unwanted transitional mechanisms are avoided. The studied boundary layers have been developed on a flat plate, by imposing a pressure gradient by means of contoured walls. They generate an initial acceleration region followed by a deceleration zone. The initial region is designed to obtain at the beginning of the deceleration the Blasius profile, characterized by its momentum thickness, and an edge boundary layer velocity, defining the problem characteristic velocity. The deceleration region is designed to obtain a linear evolution of the edge velocity, thereby defining the characteristic length of the problem. Several experimental techniques, both intrusive (hot wire anemometry, total pressure probes) as nonintrusive (PIV and LDV anemometry, high-speed filming), have been used in order to take advantage of each of them and allow cross-validation of the results. Once the boundary layer at the deceleration beginning has been characterized, ensuring the desired integral parameters and level of disturbance, the evolution of the laminar boundary layer up to the point of separation is studied. It has been compared with integral methods, and numerical simulations. In view of the results a new model for this evolution is proposed. Downstream from the separation, the flow near to the wall is configured as a shear layer that encloses low momentum recirculating fluid. The region where the shear layer remains laminar tends to be positioned to compensate the adverse pressure gradient associated with the imposed deceleration. Under these conditions, the momentum thickness remains almost constant. This laminar shear layer region extends up to where transitional phenomena appear, extension that scales with the momentum thickness at separation. These transitional phenomena are of inviscid type, similar to those found in free shear layers. The transitional region analysis begins with a study of the disturbances evolution in the linear growth region and the comparison of experimental results with a numerical model based on Linear Stability Theory for parallel flows and with data from other authors. The results’ coalescence for both the disturbances growth and the excited frequencies is stated. For the transition final stages the vorticity concentration into vortex blobs is found, analogously to what happens in free shear layers. Unlike these, the presence of the wall and the pressure gradient make the large scale structures to move towards the wall and quickly disappear under certain circumstances. In these cases, the recirculating flow is confined into a closed region saying the bubble is closed or the boundary layer reattaches. From the reattachment point, the fluid shows a configuration in the vicinity of the wall traditionally considered as turbulent. It has been observed that existing integral methods for turbulent boundary layers do not fit well to the experimental results, due to these methods being valid only for fully developed turbulent flow. Nevertheless, it has been found that downstream from the reattachment point the velocity profiles are self-similar, and a model has been proposed for the evolution of the integral parameters of the boundary layer in this region. Finally, the phenomenon known as bubble burst is analyzed. It has been checked the validity of existing models in literature and a new one is proposed. This phenomenon is blamed to the inability of the large scale structures formed after the transition to overcome with the adverse pressure gradient, move towards the wall and close the bubble. El estudio de capas límites transicionales con separación es de gran relevancia en distintas aplicaciones tecnológicas. Particularmente, en tecnología aeronáutica, aparecen en procesos claves, tales como el flujo alrededor de alas o álabes de turbomaquinaria. El objetivo de esta tesis es el estudio de estos flujos en situaciones representativas de las aplicaciones tecnológicas, ganando por un lado conocimiento sobre la fenomenología y los procesos físicos que aparecen y, por otra parte, desarrollando un modelo sencillo para el escalado de los mismos. Para conseguir este objetivo se han realizado ensayos en una instalación experimental de baja velocidad específicamente diseñada para asegurar un flujo homogéneo y con bajo nivel de perturbaciones, de modo que se evita el disparo de mecanismos transicionales no deseados. La capa límite bajo estudio se ha desarrollado sobre una placa plana, imponiendo un gradiente de presión a la misma por medio de paredes de geometría especificada. éstas generan una región inicial de aceleración seguida de una zona de deceleración. La región inicial se diseña para tener en al inicio de la deceleración un perfil de capa límite de Blasius, caracterizado por su espesor de cantidad de movimiento, y una cierta velocidad externa a la capa límite que se considera la velocidad característica del problema. La región de deceleración está concebida para que la variación de la velocidad externa a la capa límite sea lineal, definiendo de esta forma una longitud característica del problema. Los ensayos se han realizado explotando varias técnicas experimentales, tanto intrusivas (anemometría de hilo caliente, sondas de presión total) como no intrusivas (anemometrías láser y PIV, filmación de alta velocidad), de cara a aprovechar las ventajas de cada una de ellas y permitir validación cruzada de resultados entre las mismas. Caracterizada la capa límite al comienzo de la deceleración, y garantizados los parámetros integrales y niveles de perturbación deseados se procede al estudio de la zona de deceleración. Se presenta en la tesis un análisis de la evolución de la capa límite laminar desde el inicio de la misma hasta el punto de separación, comparando con métodos integrales, simulaciones numéricas, y proponiendo un nuevo modelo para esta evolución. Aguas abajo de la separación, el flujo en las proximidades de la pared se configura como una capa de cortadura que encierra una región de fluido recirculatorio de baja cantidad de movimiento. Se ha caracterizado la región en que dicha capa de cortadura permanece laminar, encontrando que se posiciona de modo que compensa el gradiente adverso de presión asociado a la deceleración de la corriente. En estas condiciones, el espesor de cantidad de movimiento permanece prácticamente constante y esta capa de cortadura laminar se extiende hasta que los fenómenos transicionales aparecen. Estos fenómenos son de tipo no viscoso, similares a los que aparecen en una capa de cortadura libre. El análisis de la región transicional comienza con un estudio de la evolución de las vii viii RESUMEN perturbaciones en la zona de crecimiento lineal de las mismas y la comparación de los resultados experimentales con un modelo numérico y con datos de otros autores. La coalescencia de los resultados tanto para el crecimiento de las perturbaciones como para las frecuencias excitadas queda demostrada. Para los estadios finales de la transición se observa la concentración de la vorticidad en torbellinos, de modo análogo a lo que ocurre en capas de cortadura libres. A diferencia de estas, la presencia de la pared y del gradiente de presión hace que, bajo ciertas condiciones, la gran escala se desplace hacia la pared y desaparezca rápidamente. En este caso el flujo recirculatorio queda confinado en una región cerrada y se habla de cierre de la burbuja o readherencia de la capa límite. A partir del punto de readherencia se tiene una configuración fluida en las proximidades de la pared que tradicionalmente se ha considerado turbulenta. Se ha observado que los métodos integrales existentes para capas límites turbulentas no ajustan bien a las medidas experimentales realizadas, hecho imputable a que no se obtiene en dicha región un flujo turbulento plenamente desarrollado. Se ha encontrado, sin embargo, que pasado el punto de readherencia los perfiles de velocidad próximos a la pared son autosemejantes entre sí y se ha propuesto un modelo para la evolución de los parámetros integrales de la capa límite en esta región. Finalmente, el fenómeno conocido como “estallido” de la burbuja se ha analizado. Se ha comprobado la validez de los modelos existentes en la literatura y se propone uno nuevo. Este fenómeno se achaca a la incapacidad de la gran estructura formada tras la transición para vencer el gradiente adverso de presión, desplazarse hacia la pared y cerrar la burbuja.

Multi-objective Aerodynamic Optimization of High-speed Trains in Tunnels

A genetic algorithm (GA) is employed for the multi-objective shape optimization of the nose of a high-speed train. Aerodynamic problems observed at high speeds become still more relevant when traveling along a tunnel. The objective is to minimize both the aerodynamic drag and the amplitude of the pressure gradient of the compression wave when a train enters a tunnel. The main drawback of GA is the large number of evaluations need in the optimization process. Metamodels-based optimization is considered to overcome such problem. As a result, an explicit relationship between pressure gradient and geometrical parameters is obtained.

High reynolds numbers turbulent boundary layers

Esta tesis estudia las similitudes y diferencias entre los flujos turbulentos de pared de tipo externo e interno, en régimen incompresible, y a números de Reynolds moderada¬mente altos. Para ello consideramos tanto simulaciones numéricas como experimentos de capas límites con gradiente de presiones nulo y de flujos de canal, ambos a números de Reynolds en el rango δ+ ~ 500 - 2000. Estos flujos de cortadura son objeto de numerosas investigaciones debido a la gran importancia que tienen tanto a nivel tecnológico como a nivel de física fundamental. No obstante, todavía existen muchos interrogantes sobre aspectos básicos tales como la universalidad de los perfiles medios y de fluctuación de las velocidades o de la presión, tanto en la zona cercana a la pared como en la zona logarítmica, el escalado y el efecto del número de Reynolds, o las diferencias entre los flujos internos y externos en la zona exterior. En éste estudio hemos utilizado simulaciones numéricas ya existentes de canales y capas límites a números de Reynolds δ+ ~ 2000 y δ+ ~ 700, respectivamente. Para poder comparar ambos flujos a igual número de Reynolds hemos realizado una nueva simulación directa de capa límite en el rango δ+ ~ 1000-2000. Los resultados de la misma son presentados y analizados en detalle. Los datos sin postprocesar y las estadísticas ya postprocesadas están públicamente disponibles en nuestro sitio web.162 El análisis de las estadísticas usando un único punto confirma la existencia de perfiles logarítmicos para las fluctuaciones de la velocidad transversal w'2+ y de la presión p'2+ en ambos tipos de flujos, pero no para la velocidad normal v'2+ o la velocidad longitudinal u'2+. Para aceptar o rechazar la existencia de un rango logarítmico en u'2+ se requieren números de Reynolds más altos que los considerados en éste trabajo. Una de las conse¬cuencias más importantes de poseer tales perfiles es que el valor máximo de la intensidad, que se alcanza cerca de la pared, depende explícitamente del número de Reynolds. Esto ha sido confirmado tras analizar un gran número de datos experimentales y numéricos, cor¬roborando que el máximo de u'2+, p/2+, y w'2+ aumenta proporcionalmente con el log(δ+). Por otro lado, éste máximo es más intenso en los flujos externos que en los internos. La máxima diferencia ocurre en torno a y/δ ~ 0.3-0.5, siendo esta altura prácticamente independiente del número de Reynolds considerado. Estas diferencias se originan como consecuencia del carácter intermitente de las capas límites, que es inexistente en los flujos internos. La estructura de las fluctuaciones de velocidad y de presión, junto con la de los esfuer¬zos de Reynolds, se han investigado por medio de correlaciones espaciales tridimensionales considerando dos puntos de medida. Hemos obtenido que el tamaño de las mismas es gen¬eralmente mayor en canales que en capas límites, especialmente en el caso de la correlación longitudinal Cuu en la dirección del flujo. Para esta correlación se demuestra que las es¬tructuras débilmente correladas presentan longitudes de hasta 0(75), en el caso de capas límites, y de hasta 0(185) en el caso de canales. Estas longitudes se obtienen respecti-vamente en la zona logarítmica y en la zona exterior. Las longitudes correspondientes en la dirección transversal son significativamente menores en ambos flujos, 0(5 — 25). La organización espacial de las correlaciones es compatible con la de una pareja de rollos casi paralelos con dimensiones que escalan en unidades exteriores. Esta organización se mantiene al menos hasta y ~ 0.65, altura a la cual las capas límites comienzan a organi¬zarse en rollos transversales. Este comportamiento es sin embargo más débil en canales, pudiéndose observar parcialmente a partir de y ~ 0.85. Para estudiar si estas estructuras están onduladas a lo largo de la dirección transver¬sal, hemos calculado las correlaciones condicionadas a eventos intensos de la velocidad transversal w'. Estas correlaciones revelan que la ondulación de la velocidad longitudinal aumenta conforme nos alejamos de la pared, sugiriendo que las estructuras están más alineadas en la zona cercana a la pared que en la zona lejana a ella. El por qué de esta ondulación se encuentra posiblemente en la configuración a lo largo de diagonales que presenta w'. Estas estructuras no sólo están onduladas, sino que también están inclinadas respecto a la pared con ángulos que dependen de la variable considerada, de la altura, y de el contorno de correlación seleccionado. Por encima de la zona tampón e independien¬temente del número de Reynolds y tipo de flujo, Cuu presenta una inclinación máxima de unos 10°, las correlaciones Cvv y Cm son esencialmente verticales, y Cww está inclinada a unos 35°. Summary This thesis studies the similitudes and differences between external and internal in¬compressible wall-bounded turbulent flows at moderately-high Reynolds numbers. We consider numerical and experimental zero-pressure-gradient boundary layers and chan¬nels in the range of δ+ ~ 500 — 2000. These shear flows are subjects of intensive research because of their technological importance and fundamental physical interest. However, there are still open questions regarding basic aspects such as the universality of the mean and fluctuating velocity and pressure profiles at the near-wall and logarithmic regions, their scaling and the effect of the Reynolds numbers, or the differences between internal and external flows at the outer layer, to name but a few. For this study, we made use of available direct numerical simulations of channel and boundary layers reaching δ+ ~ 2000 and δ+ ~ 700, respectively. To fill the gap in the Reynolds number, a new boundary layer simulation in the range δ+ ~ 1000-2000 is presented and discussed. The original raw data and the post-processed statistics are publicly available on our website.162 The analysis of the one-point statistic confirms the existence of logarithmic profiles for the spanwise w'2+ and pressure p'2+ fluctuations for both type of flows, but not for the wall-normal v'2+ or the streamwise u'2+ velocities. To accept or reject the existence of a logarithmic range in u'2+ requires higher Reynolds numbers than the ones considered in this work. An important consequence of having such profiles is that the maximum value of the intensities, reached near the wall, depends on the Reynolds number. This was confirmed after surveying a wide number of experimental and numerical datasets, corrob¬orating that the maximum of ul2+, p'2+, and w'2+ increases proportionally to log(δ+). On the other hand, that maximum is more intense in external flows than in internal ones, differing the most around y/δ ~ 0.3-0.5, and essentially independent of the Reynolds number. We discuss that those differences are originated as a consequence of the inter¬mittent character of boundary layers that is absent in internal flows. The structure of the velocity and pressure fluctuations, together with those of the Reynolds shear stress, were investigated using three-dimensional two-point spatial correlations. We find that the correlations extend over longer distances in channels than in boundary layers, especially in the case of the streamwise correlation Cuu in the flow direc-tion. For weakly correlated structures, the maximum streamwise length of Cuu is O(78) for boundary layers and O(188) for channels, attained at the logarithmic and outer regions respectively. The corresponding lengths for the transverse velocities and for the pressure are shorter, 0(8 — 28), and of the same order for both flows. The spatial organization of the velocity correlations is shown to be consistent with a pair of quasi-streamwise rollers that scales in outer units. That organization is observed until y ~ 0.68, from which boundary layers start to organize into spanwise rollers. This effect is weaker in channels, and it appears at y ~ 0.88. We present correlations conditioned to intense events of the transversal velocity, w', to study if these structures meander along the spanwise direction. The results indicate that the streamwise velocity streaks increase their meandering proportionally to the distance to the wall, suggesting that the structures are more aligned close to the wall than far from it. The reason behind this meandering is probably due to the characteristic organization along diagonals of w'. These structures not only meander along the spanwise direction, but they are also inclined to the wall at angles that depend on the distance from the wall, on the variable being considered, and on the correlation level used to define them. Above the buffer layer and independent of the Reynolds numbers and type of flow, the maximum inclination of Cuu is about 10°, Cvv and Cpp are roughly vertical, and Cww is inclined by 35°.

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.

A code for direct numerical simulation of turbulent boundary layers at high Reynolds numbers in BG/P supercomputers

A new high-resolution code for the direct numerical simulation of a zero pressure gradient turbulent boundary layers over a flat plate has been developed. Its purpose is to simulate a wide range of Reynolds numbers from Reθ = 300 to 6800 while showing a linear weak scaling up to 32,768 cores in the BG/P architecture. Special attention has been paid to the generation of proper inflow boundary conditions. The results are in good agreement with existing numerical and experimental data sets.

Entrainment Effects in Turbulent Boundary Layers

Esta tesis estudia el comportamiento de la región exterior de una capa límite turbulenta sin gradientes de presiones. Se ponen a prueba dos teorías relativamente bien establecidas. La teoría de semejanza para la pared supone que en el caso de haber una pared rugosa, el fluido sólo percibe el cambio en la fricción superficial que causa, y otros efectos secundarios quedarán confinados a una zona pegada a la pared. El consenso actual es que dicha teoría es aproximadamente cierta. En el extremo exterior de la capa límite existe una región producida por la interacción entre las estructuras turbulentas y el flujo irrotacional de la corriente libre llamada interfaz turbulenta/no turbulenta. La mayoría de los resultados al respecto sugieren la presencia de fuerzas de cortadura ligeramente más intensa, lo que la hace distinta al resto del flujo turbulento. Las propiedades de esa región probablemente cambien si la velocidad de crecimiento de la capa límite aumenta, algo que puede conseguirse aumentando la fricción en la pared. La rugosidad y la ingestión de masa están entonces relacionadas, y el comportamiento local de la interfaz turbulenta/no turbulenta puede explicar el motivo por el que las capas límite sobre paredes rugosas no se comportan como en el caso de tener paredes lisas precisamente en la zona exterior. Para estudiar las capas límite a números de Reynolds lo suficientemente elevados, se ha desarrollado un nuevo código de alta resolución para la simulación numérica directa de capas límite turbulentas sin gradiente de presión. Dicho código es capaz de simular capas límite en un intervalo de números de Reynolds entre ReT = 100 — 2000 manteniendo una buena escalabilidad hasta los dos millones de hilos en superordenadores de tipo Blue Gene/Q. Se ha guardado especial atención a la generación de condiciones de contorno a la entrada correctas. Los resultados obtenidos están en concordancia con los resultados previos, tanto en el caso de simulaciones como de experimentos. La interfaz turbulenta/no turbulenta de una capa límite se ha analizado usando un valor umbral del módulo de la vorticidad. Dicho umbral se considera un parámetro para analizar cada superficie obtenida de un contorno del módulo de la vorticidad. Se han encontrado dos regímenes distintos en función del umbral escogido con propiedades opuestas, separados por una transición topológica gradual. Las características geométricas de la zona escalan con o99 cuando u^/isdgg es la unidad de vorticidad. Las propiedades del íluido relativas a la posición del contorno de vorticidad han sido analizados para una serie de umbrales utilizando el campo de distancias esféricas, que puede obtenerse con independencia de la complejidad de la superficie de referencia. Las propiedades del fluido a una distancia dada del inerfaz también dependen del umbral de vorticidad, pero tienen características parecidas con independencia del número de Reynolds. La interacción entre la turbulencia y el flujo no turbulento se restringe a una zona muy fina con un espesor del orden de la escala de Kolmogorov local. Hacia el interior del flujo turbulento las propiedades son indistinguibles del resto de la capa límite. Se ha simulado una capa límite sin gradiente de presiones con una fuerza volumétrica cerca de la pared. La el forzado ha sido diseñado para aumentar la fricción en la pared sin introducir ningún efecto geométrico obvio. La simulación consta de dos dominios, un primer dominio más pequeño y a baja resolución que se encarga de generar condiciones de contorno correctas, y un segundo dominio mayor y a alta resolución donde se aplica el forzado. El estudio de los perfiles y los coeficientes de autocorrelación sugieren que los dos casos, el liso y el forzado, no colapsan más allá de la capa logarítmica por la complejidad geométrica de la zona intermitente, y por el hecho que la distancia a la pared no es una longitud característica. Los efectos causados por la geometría de la zona intermitente pueden evitarse utilizando el interfaz como referencia, y la distancia esférica para el análisis de sus propiedades. Las propiedades condicionadas del flujo escalan con 5QQ y u/uT, las dos únicas escalas contenidas en el modelo de semejanza de pared de Townsend, consistente con estos resultados. ABSTRACT This thesis studies the characteristics of the outer region of zero-pressure-gradient turbulent boundary layers at moderate Reynolds numbers. Two relatively established theories are put to test. The wall similarity theory states that with the presence of roughness, turbulent motion is mostly affected by the additional drag caused by the roughness, and that other secondary effects are restricted to a region very close to the wall. The consensus is that this theory is valid, but only as a first approximation. At the edge of the boundary layer there is a thin layer caused by the interaction between the turbulent eddies and the irroational fluid of the free stream, called turbulent/non-turbulent interface. The bulk of results about this layer suggest the presence of some localized shear, with properties that make it distinguishable from the rest of the turbulent flow. The properties of the interface are likely to change if the rate of spread of the turbulent boundary layer is amplified, an effect that is usually achieved by increasing the drag. Roughness and entrainment are therefore linked, and the local features of the turbulent/non-turbulent interface may explain the reason why rough-wall boundary layers deviate from the wall similarity theory precisely far from the wall. To study boundary layers at a higher Reynolds number, a new high-resolution code for the direct numerical simulation of a zero pressure gradient turbulent boundary layers over a flat plate has been developed. This code is able to simulate a wide range of Reynolds numbers from ReT =100 to 2000 while showing a linear weak scaling up to around two million threads in the BG/Q architecture. Special attention has been paid to the generation of proper inflow boundary conditions. The results are in good agreement with existing numerical and experimental data sets. The turbulent/non-turbulent interface of a boundary layer is analyzed by thresholding the vorticity magnitude field. The value of the threshold is considered a parameter in the analysis of the surfaces obtained from isocontours of the vorticity magnitude. Two different regimes for the surface can be distinguished depending on the threshold, with a gradual topological transition across which its geometrical properties change significantly. The width of the transition scales well with oQg when u^/udgg is used as a unit of vorticity. The properties of the flow relative to the position of the vorticity magnitude isocontour are analyzed within the same range of thresholds, using the ball distance field, which can be obtained regardless of the size of the domain and complexity of the interface. The properties of the flow at a given distance to the interface also depend on the threshold, but they are similar regardless of the Reynolds number. The interaction between the turbulent and the non-turbulent flow occurs in a thin layer with a thickness that scales with the Kolmogorov length. Deeper into the turbulent side, the properties are undistinguishable from the rest of the turbulent flow. A zero-pressure-gradient turbulent boundary layer with a volumetric near-wall forcing has been simulated. The forcing has been designed to increase the wall friction without introducing any obvious geometrical effect. The actual simulation is split in two domains, a smaller one in charge of the generation of correct inflow boundary conditions, and a second and larger one where the forcing is applied. The study of the one-point and twopoint statistics suggest that the forced and the smooth cases do not collapse beyond the logarithmic layer may be caused by the geometrical complexity of the intermittent region, and by the fact that the scaling with the wall-normal coordinate is no longer present. The geometrical effects can be avoided using the turbulent/non-turbulent interface as a reference frame, and the minimum distance respect to it. The conditional analysis of the vorticity field with the alternative reference frame recovers the scaling with 5QQ and v¡uT already present in the logarithmic layer, the only two length-scales allowed if Townsend’s wall similarity hypothesis is valid.