Calibración geométrica de cámaras no métricas. Estudio de metodologías y modelos matemáticos de distorsión

La Fotogrametría, como ciencia y técnica de obtención de información tridimensional del espacio objeto a partir de imágenes bidimensionales, requiere de medidas de precisión y en ese contexto, la calibración geométrica de cámaras ocupa un lugar importante. El conocimiento de la geometría interna de la cámara es fundamental para lograr mayor precisión en las medidas realizadas. En Fotogrametría Aérea se utilizan cámaras métricas (fabricadas exclusivamente para aplicaciones cartográficas), que incluyen objetivos fotográficos con sistemas de lentes complejos y de alta calidad. Pero en Fotogrametría de Objeto Cercano se está trabajando cada vez con más asiduidad con cámaras no métricas, con ópticas de peor calidad que exigen una calibración geométrica antes o después de cada trabajo. El proceso de calibración encierra tres conceptos fundamentales: modelo de cámara, modelo de distorsión y método de calibración. El modelo de cámara es un modelo matemático que aproxima la transformación proyectiva original a la realidad física de las lentes. Ese modelo matemático incluye una serie de parámetros entre los que se encuentran los correspondientes al modelo de distorsión, que se encarga de corregir los errores sistemáticos de la imagen. Finalmente, el método de calibración propone el método de estimación de los parámetros del modelo matemático y la técnica de optimización a emplear. En esta Tesis se propone la utilización de un patrón de calibración bidimensional que se desplaza en la dirección del eje óptico de la cámara, ofreciendo así tridimensionalidad a la escena fotografiada. El patrón incluye un número elevado de marcas, lo que permite realizar ensayos con distintas configuraciones geométricas. Tomando el modelo de proyección perspectiva (o pinhole) como modelo de cámara, se realizan ensayos con tres modelos de distorsión diferentes, el clásico de distorsión radial y tangencial propuesto por D.C. Brown, una aproximación por polinomios de Legendre y una interpolación bicúbica. De la combinación de diferentes configuraciones geométricas y del modelo de distorsión más adecuado, se llega al establecimiento de una metodología de calibración óptima. Para ayudar a la elección se realiza un estudio de las precisiones obtenidas en los distintos ensayos y un control estereoscópico de un panel test construido al efecto. ABSTRACT Photogrammetry, as science and technique for obtaining three-dimensional information of the space object from two-dimensional images, requires measurements of precision and in that context, the geometric camera calibration occupies an important place. The knowledge of the internal geometry of the camera is fundamental to achieve greater precision in measurements made. Metric cameras (manufactured exclusively for cartographic applications), including photographic lenses with complex lenses and high quality systems are used in Aerial Photogrammetry. But in Close Range Photogrammetry is working increasingly more frequently with non-metric cameras, worst quality optical components which require a geometric calibration before or after each job. The calibration process contains three fundamental concepts: camera model, distortion model and method of calibration. The camera model is a mathematical model that approximates the original projective transformation to the physical reality of the lenses. The mathematical model includes a series of parameters which include the correspondents to the model of distortion, which is in charge of correcting the systematic errors of the image. Finally, the calibration method proposes the method of estimation of the parameters of the mathematical modeling and optimization technique to employ. This Thesis is proposing the use of a pattern of two dimensional calibration that moves in the direction of the optical axis of the camera, thus offering three-dimensionality to the photographed scene. The pattern includes a large number of marks, which allows testing with different geometric configurations. Taking the projection model perspective (or pinhole) as a model of camera, tests are performed with three different models of distortion, the classical of distortion radial and tangential proposed by D.C. Brown, an approximation by Legendre polynomials and bicubic interpolation. From the combination of different geometric configurations and the most suitable distortion model, brings the establishment of a methodology for optimal calibration. To help the election, a study of the information obtained in the various tests and a purpose built test panel stereoscopic control is performed.

An engineering modification of the blade element momentum combination applicable in vertical descent flight :fundamentals and validation

An engineering modification of blade element/momentum theory is applied to describe the vertical autorotation of helicopter rotors. A full non‐linear aerodynamic model is considered for the airfoils, taking into account the dependence of lift and drag coefficients on both the angle of attack and the Reynolds number. The proposed model, which has been validated in previous work, has allowed the identification of different autorotation modes, which depend on the descent velocity and the twist of the rotor blades. These modes present different radial distributions of driven and driving blade regions, as well as different radial upwash/downwash patterns. The number of blade sections with zero tangential force, the existence of a downwash region in the rotor disk, the stability of the autorotation state, and the overall rotor autorotation efficiency, are all analyzed in terms of the flight velocity and the characteristics of the rotor. It is shown that, in vertical autorotation, larger blade twist leads to smaller values of descent velocity for a given thrust generated by the rotor in the autorotational state.

An engineering modification of the blade element momentum equation for vertical descent : an autoration case study

An engineering modification of blade element/momentum theory is applied to describe the vertical autorotation of helicopter rotors. A full non-linear aerodynamic model is considered for the airfoils, taking into account the dependence of lift and drag coefficients on both the angle of attack and the Reynolds number. The proposed model, which has been validated in previous work, has allowed the identification of different autorotation modes, which depend on the descent velocity and the twist of the rotor blades. These modes present different radial distributions of driven and driving blade regions, as well as different radial upwash/downwash patterns. The number of blade sections with zero tangential force, the existence of a downwash region in the rotor disk, the stability of the autorotation state, and the overall rotor autorotation efficiency, are all analyzed in terms of the flight velocity and the characteristics of the rotor. It is shown that, in vertical autorotation, larger blade twist leads to smaller values of descent velocity for a given thrust generated by the rotor in the autorotational state.

Anatomía comparada de la madera de Cupressaceae y su correspondencia con los estudios de filogenia

La familia Cupressaceae incluye un total de 133 especies agrupadas en 30 géneros, 17 de los cuales son monospecíficos. Esta familia se encuentra representada en todos los continentes salvo en la Antártida. Sus especies se distribuyen en distintas regiones climáticas, y en altitudes que varían desde el nivel del mar hasta los 5.000 m. La falta de descripción anatómica de muchos de los géneros y especies de Cupressaceae es notable, así como la contradicción que aparece entre distintas investigaciones sobre las características anatómicas de la madera descritas para cada especie. Este estudio describe la anatomía de la madera de Cupressaceae y analiza las características que podrían representar sinapomorfías de los clados delimitados en los estudios filogenéticos. Siguiendo los métodos tradicionales de preparación y descripción de la madera a nivel microscópico, se ha estudiado la madera de 113 especies de los 30 géneros de Cupressaceae. Para ello se han empleado muestras de madera de origen trazable, procedentes de colecciones de madera de distintas instituciones internacionales. Se ha empleado una robusta filogenia molecular para la reconstrucción de los caracteres ancestrales. La anatomía de la madera de los 30 géneros de Cupressaceae, pone de manifiesto la gran homogeneidad de la familia, caracterizada por la presencia de traqueidas axiales sin engrosamientos helicoidales, parénquima radial con paredes horizontales lisas, punteaduras del campo de cruce de tipo cupresoide y la carencia de canales resiníferos fisiológicos. Además, todos presentan parénquima axial (salvo Neocallitropsis, Thuja y Xanthocyparis), punteaduras radiales areoladas con toro definido (salvo Thuja y Thujopsis), siendo habitual la presencia de punteaduras areoladas en las paredes tangenciales de la madera tardía, y verrugosidades en la cara interna de las traqueidas (salvo Ca. macleayana, Libocedrus, Papuacedrus y Neocallitropsis). Los radios leñosos son homogéneos y están compuestos de parénquima radial (con la presencia de traqueidas radiales en algunas especies de Cupressus, Sequoia, Thujopsis y X. nootkatensis) con paredes finales lisas o lisas y noduladas (exclusivamente noduladas en Cal. macrolepis, C. bakeri y en la mayoría de especies de Juniperus), y el rango de altura de los radios leñosos se encuentra entre 5 y 15 células. Se consideran posibles sinapomorfismos de Cupressaceae la presencia de verrugosidades en la cara interna de las traqueidas, la presencia de traqueidas axiales sin engrosamientos helicoidales, la presencia de parénquima axial, la presencia de radios leñosos homogéneos (compuestos únicamente de parénquima radial), la tipología de las paredes horizontales del parénquima radial, las punteaduras del campo de cruce de tipo cupresoide y la ausencia de canales resiníferos fisiológicos, pero lo que realmente diferencia a este grupo de coníferas es la simultaneidad de todos estos caracteres en sus maderas. Como sinapomorfías específicas por clados se proponen: la ausencia de toro definido y muescas en el borde de las punteaduras en Thuja-Thujopsis, la existencia de extensiones de toro en Diselma-Fitzroya-Widdringtonia; la presencia de engrosamientos callitroides en Callitris-Actinostrobus; la presencia de espacios intercelulares y las muescas en el borde de las punteaduras en el clado formado por el género Juniperus y las especies de Cupressus en la región oriental; la presencia de paredes finales del parénquima radial tanto lisas como noduladas en los clados formados por el género Xanthocyparis y las especies de Cupressus en la región occidental y en Fitzroya-Diselma; y por último, la presencia de punteaduras del campo de cruce de tipo taxodioide en los clados taxodioid y sequoioid. ABSTRACT The Cupressaceae family comprises 133 species grouped into 30 genera, 17 of which are monotypic. The family is represented in all continents except Antarctica. Its species are distributed in various climate zones and at altitudes from sea level to 5,000 m. There is a considerable lack of anatomical descriptions for many genera and species of Cupressaceae and much contradiction between studies about the wood anatomical features described for each species. This study describes the wood anatomy of Cupressaceae and analyses the features that could represent synapomorphies of the clades recovered in phylogenetic studies. Following the traditional methods of preparation and description of wood at microscopic level, a study was made of the wood of 113 species of the 30 Cupressaceae genera. The study samples had traceable origins and came from wood collections held at various international institutions. A robust molecular phylogeny was used for ancestral state reconstruction. The wood anatomy of the 30 genera of the Cupressaceae shows the high homogeneity of the family, which is characterised by the presence of axial tracheids without helical thickenings, smooth horizontal walls of ray parenchyma cells, cupressoid cross-field pits, and the absence of physiological resin canals. In addition, they all have axial parenchyma (except Neocallitropsis, Thuja and Xanthocyparis), a warty layer on the inner wall of the tracheids (except Ca. macleayana, Libocedrus, Papuacedrus and Neocallitropsis) and tracheid pitting in radial walls with a well defined torus (except Thuja and Thujopsis); tracheid pitting in the tangential walls of the latewood is common. Rays are homogeneous and are composed of ray parenchyma (with the presence of ray tracheids in some species of Cupressus, Sequoia, Thujopsis and X. nootkatensis), with smooth end walls or both smooth and nodular end walls (exclusively nodular in Cal. macrolepis, C. bakeri and most Juniperus species), and ray height range is 5 to 15 cells. Possible synapomorphies of Cupressaceae are the presence of a warty layer on the inner layer of the tracheids, axial tracheids without helical thickenings, the presence of axial parenchyma, homogeneous rays (composed exclusively of ray parenchyma), the typology of the horizontal walls of ray parenchyma cells, cupressoid cross-field pits and the absence of physiological resin canals, but what truly differentiates this group of softwoods is the co-occurrence of all these features in their wood. The following are proposed as clade-specific synapomorphies: absence of a well-defined torus and presence of pits with notched borders in Thuja-Thujopsis, torus extensions in Diselma-Fitzroya-Widdringtonia; callitroid thickenings in Callitris-Actinostrobus; intercellular spaces and pits with notched borders in the clade formed by the genus Juniperus and the species of Cupressus in the eastern region; smooth and nodular ray parenchyma end walls in the clades formed by the genus Xanthocyparis and the species of Cupressus in the western region and in Fitzroya-Diselma, and taxodioid cross-field pits in the taxodioid and sequoioid clades.

Coherent structures in statistically-stationary homogeneous shear turbulence

La región cerca de la pared de flujos turbulentos de pared ya está bien conocido debido a su bajo número de Reynolds local y la separación escala estrecha. La región lejos de la pared (capa externa) no es tan interesante tampoco, ya que las estadísticas allí se escalan bien por las unidades exteriores. La región intermedia (capa logarítmica), sin embargo, ha estado recibiendo cada vez más atención debido a su propiedad auto-similares. Además, de acuerdo a Flores et al. (2007) y Flores & Jiménez (2010), la capa logarítmica es más o menos independiente de otras capas, lo que implica que podría ser inspeccionado mediante el aislamiento de otras dos capas, lo que reduciría significativamente los costes computacionales para la simulación de flujos turbulentos de pared. Algunos intentos se trataron después por Mizuno & Jiménez (2013), quien simulan la capa logarítmica sin la región cerca de la pared con estadísticas obtenidas de acuerdo razonablemente bien con los de las simulaciones completas. Lo que más, la capa logarítmica podría ser imitado por otra turbulencia sencillo de cizallamiento de motor. Por ejemplo, Pumir (1996) encontró que la turbulencia de cizallamiento homogéneo estadísticamente estacionario (SS-HST) también irrumpe, de una manera muy similar al proceso de auto-sostenible en flujos turbulentos de pared. Según los consideraciones arriba, esta tesis trata de desvelar en qué medida es la capa logarítmica de canales similares a la turbulencia de cizalla más sencillo, SS-HST, mediante la comparación de ambos cinemática y la dinámica de las estructuras coherentes en los dos flujos. Resultados sobre el canal se muestran mediante Lozano-Durán et al. (2012) y Lozano-Durán & Jiménez (2014b). La hoja de ruta de esta tarea se divide en tres etapas. En primer lugar, SS-HST es investigada por medio de un código nuevo de simulación numérica directa, espectral en las dos direcciones horizontales y compacto-diferencias finitas en la dirección de la cizalla. Sin utiliza remallado para imponer la condición de borde cortante periódica. La influencia de la geometría de la caja computacional se explora. Ya que el HST no tiene ninguna longitud característica externa y tiende a llenar el dominio computacional, las simulaciopnes a largo plazo del HST son ’mínimos’ en el sentido de que contiene sólo unas pocas estructuras media a gran escala. Se ha encontrado que el límite principal es el ancho de la caja de la envergadura, Lz, que establece las escalas de longitud y velocidad de la turbulencia, y que las otras dos dimensiones de la caja debe ser suficientemente grande (Lx > 2LZ, Ly > Lz) para evitar que otras direcciones estando limitado también. También se encontró que las cajas de gran longitud, Lx > 2Ly, par con el paso del tiempo la condición de borde cortante periódica, y desarrollar fuertes ráfagas linealizadas no físicos. Dentro de estos límites, el flujo muestra similitudes y diferencias interesantes con otros flujos de cizalla, y, en particular, con la capa logarítmica de flujos turbulentos de pared. Ellos son exploradas con cierto detalle. Incluyen un proceso autosostenido de rayas a gran escala y con una explosión cuasi-periódica. La escala de tiempo de ruptura es de aproximadamente universales, ~20S~l(S es la velocidad de cizallamiento media), y la disponibilidad de dos sistemas de ruptura diferentes permite el crecimiento de las ráfagas a estar relacionado con algo de confianza a la cizalladura de turbulencia inicialmente isotrópico. Se concluye que la SS-HST, llevado a cabo dentro de los parámetros de cílculo apropiados, es un sistema muy prometedor para estudiar la turbulencia de cizallamiento en general. En segundo lugar, las mismas estructuras coherentes como en los canales estudiados por Lozano-Durán et al. (2012), es decir, grupos de vórticidad (fuerte disipación) y Qs (fuerte tensión de Reynolds tangencial, -uv) tridimensionales, se estudia mediante simulación numérica directa de SS-HST con relaciones de aspecto de cuadro aceptables y número de Reynolds hasta Rex ~ 250 (basado en Taylor-microescala). Se discute la influencia de la intermitencia de umbral independiente del tiempo. Estas estructuras tienen alargamientos similares en la dirección sentido de la corriente a las familias separadas en los canales hasta que son de tamaño comparable a la caja. Sus dimensiones fractales, longitudes interior y exterior como una función del volumen concuerdan bien con sus homólogos de canales. El estudio sobre sus organizaciones espaciales encontró que Qs del mismo tipo están alineados aproximadamente en la dirección del vector de velocidad en el cuadrante al que pertenecen, mientras Qs de diferentes tipos están restringidos por el hecho de que no debe haber ningún choque de velocidad, lo que hace Q2s (eyecciones, u < 0,v > 0) y Q4s (sweeps, u > 0,v < 0) emparejado en la dirección de la envergadura. Esto se verifica mediante la inspección de estructuras de velocidad, otros cuadrantes como la uw y vw en SS-HST y las familias separadas en el canal. La alineación sentido de la corriente de Qs ligada a la pared con el mismo tipo en los canales se debe a la modulación de la pared. El campo de flujo medio condicionado a pares Q2-Q4 encontró que los grupos de vórticidad están en el medio de los dos, pero prefieren los dos cizalla capas alojamiento en la parte superior e inferior de Q2s y Q4s respectivamente, lo que hace que la vorticidad envergadura dentro de las grupos de vórticidad hace no cancele. La pared amplifica la diferencia entre los tamaños de baja- y alta-velocidad rayas asociados con parejas de Q2-Q4 se adjuntan como los pares alcanzan cerca de la pared, el cual es verificado por la correlación de la velocidad del sentido de la corriente condicionado a Q2s adjuntos y Q4s con diferentes alturas. Grupos de vórticidad en SS-HST asociados con Q2s o Q4s también están flanqueadas por un contador de rotación de los vórtices sentido de la corriente en la dirección de la envergadura como en el canal. La larga ’despertar’ cónica se origina a partir de los altos grupos de vórticidad ligada a la pared han encontrado los del Álamo et al. (2006) y Flores et al. (2007), que desaparece en SS-HST, sólo es cierto para altos grupos de vórticidad ligada a la pared asociados con Q2s pero no para aquellos asociados con Q4s, cuyo campo de flujo promedio es en realidad muy similar a la de SS-HST. En tercer lugar, las evoluciones temporales de Qs y grupos de vórticidad se estudian mediante el uso de la método inventado por Lozano-Durán & Jiménez (2014b). Las estructuras se clasifican en las ramas, que se organizan más en los gráficos. Ambas resoluciones espaciales y temporales se eligen para ser capaz de capturar el longitud y el tiempo de Kolmogorov puntual más probable en el momento más extrema. Debido al efecto caja mínima, sólo hay un gráfico principal consiste en casi todas las ramas, con su volumen y el número de estructuras instantáneo seguien la energía cinética y enstrofía intermitente. La vida de las ramas, lo que tiene más sentido para las ramas primarias, pierde su significado en el SS-HST debido a las aportaciones de ramas primarias al total de Reynolds estrés o enstrofía son casi insignificantes. Esto también es cierto en la capa exterior de los canales. En cambio, la vida de los gráficos en los canales se compara con el tiempo de ruptura en SS-HST. Grupos de vórticidad están asociados con casi el mismo cuadrante en términos de sus velocidades medias durante su tiempo de vida, especialmente para los relacionados con las eyecciones y sweeps. Al igual que en los canales, las eyecciones de SS-HST se mueven hacia arriba con una velocidad promedio vertical uT (velocidad de fricción) mientras que lo contrario es cierto para los barridos. Grupos de vórticidad, por otra parte, son casi inmóvil en la dirección vertical. En la dirección de sentido de la corriente, que están advección por la velocidad media local y por lo tanto deforman por la diferencia de velocidad media. Sweeps y eyecciones se mueven más rápido y más lento que la velocidad media, respectivamente, tanto por 1.5uT. Grupos de vórticidad se mueven con la misma velocidad que la velocidad media. Se verifica que las estructuras incoherentes cerca de la pared se debe a la pared en vez de pequeño tamaño. Los resultados sugieren fuertemente que las estructuras coherentes en canales no son especialmente asociado con la pared, o incluso con un perfil de cizalladura dado. ABSTRACT Since the wall-bounded turbulence was first recognized more than one century ago, its near wall region (buffer layer) has been studied extensively and becomes relatively well understood due to the low local Reynolds number and narrow scale separation. The region just above the buffer layer, i.e., the logarithmic layer, is receiving increasingly more attention nowadays due to its self-similar property. Flores et al. (20076) and Flores & Jim´enez (2010) show that the statistics of logarithmic layer is kind of independent of other layers, implying that it might be possible to study it separately, which would reduce significantly the computational costs for simulations of the logarithmic layer. Some attempts were tried later by Mizuno & Jimenez (2013), who simulated the logarithmic layer without the buffer layer with obtained statistics agree reasonably well with those of full simulations. Besides, the logarithmic layer might be mimicked by other simpler sheardriven turbulence. For example, Pumir (1996) found that the statistically-stationary homogeneous shear turbulence (SS-HST) also bursts, in a manner strikingly similar to the self-sustaining process in wall-bounded turbulence. Based on these considerations, this thesis tries to reveal to what extent is the logarithmic layer of channels similar to the simplest shear-driven turbulence, SS-HST, by comparing both kinematics and dynamics of coherent structures in the two flows. Results about the channel are shown by Lozano-Dur´an et al. (2012) and Lozano-Dur´an & Jim´enez (20146). The roadmap of this task is divided into three stages. First, SS-HST is investigated by means of a new direct numerical simulation code, spectral in the two horizontal directions and compact-finite-differences in the direction of the shear. No remeshing is used to impose the shear-periodic boundary condition. The influence of the geometry of the computational box is explored. Since HST has no characteristic outer length scale and tends to fill the computational domain, longterm simulations of HST are ‘minimal’ in the sense of containing on average only a few large-scale structures. It is found that the main limit is the spanwise box width, Lz, which sets the length and velocity scales of the turbulence, and that the two other box dimensions should be sufficiently large (Lx > 2LZ, Ly > Lz) to prevent other directions to be constrained as well. It is also found that very long boxes, Lx > 2Ly, couple with the passing period of the shear-periodic boundary condition, and develop strong unphysical linearized bursts. Within those limits, the flow shows interesting similarities and differences with other shear flows, and in particular with the logarithmic layer of wallbounded turbulence. They are explored in some detail. They include a self-sustaining process for large-scale streaks and quasi-periodic bursting. The bursting time scale is approximately universal, ~ 20S~l (S is the mean shear rate), and the availability of two different bursting systems allows the growth of the bursts to be related with some confidence to the shearing of initially isotropic turbulence. It is concluded that SS-HST, conducted within the proper computational parameters, is a very promising system to study shear turbulence in general. Second, the same coherent structures as in channels studied by Lozano-Dur´an et al. (2012), namely three-dimensional vortex clusters (strong dissipation) and Qs (strong tangential Reynolds stress, -uv), are studied by direct numerical simulation of SS-HST with acceptable box aspect ratios and Reynolds number up to Rex ~ 250 (based on Taylor-microscale). The influence of the intermittency to time-independent threshold is discussed. These structures have similar elongations in the streamwise direction to detached families in channels until they are of comparable size to the box. Their fractal dimensions, inner and outer lengths as a function of volume agree well with their counterparts in channels. The study about their spatial organizations found that Qs of the same type are aligned roughly in the direction of the velocity vector in the quadrant they belong to, while Qs of different types are restricted by the fact that there should be no velocity clash, which makes Q2s (ejections, u < 0, v > 0) and Q4s (sweeps, u > 0, v < 0) paired in the spanwise direction. This is verified by inspecting velocity structures, other quadrants such as u-w and v-w in SS-HST and also detached families in the channel. The streamwise alignment of attached Qs with the same type in channels is due to the modulation of the wall. The average flow field conditioned to Q2-Q4 pairs found that vortex clusters are in the middle of the pair, but prefer to the two shear layers lodging at the top and bottom of Q2s and Q4s respectively, which makes the spanwise vorticity inside vortex clusters does not cancel. The wall amplifies the difference between the sizes of low- and high-speed streaks associated with attached Q2-Q4 pairs as the pairs reach closer to the wall, which is verified by the correlation of streamwise velocity conditioned to attached Q2s and Q4s with different heights. Vortex clusters in SS-HST associated with Q2s or Q4s are also flanked by a counter rotating streamwise vortices in the spanwise direction as in the channel. The long conical ‘wake’ originates from tall attached vortex clusters found by del A´ lamo et al. (2006) and Flores et al. (2007b), which disappears in SS-HST, is only true for tall attached vortices associated with Q2s but not for those associated with Q4s, whose averaged flow field is actually quite similar to that in SS-HST. Third, the temporal evolutions of Qs and vortex clusters are studied by using the method invented by Lozano-Dur´an & Jim´enez (2014b). Structures are sorted into branches, which are further organized into graphs. Both spatial and temporal resolutions are chosen to be able to capture the most probable pointwise Kolmogorov length and time at the most extreme moment. Due to the minimal box effect, there is only one main graph consist by almost all the branches, with its instantaneous volume and number of structures follow the intermittent kinetic energy and enstrophy. The lifetime of branches, which makes more sense for primary branches, loses its meaning in SS-HST because the contributions of primary branches to total Reynolds stress or enstrophy are almost negligible. This is also true in the outer layer of channels. Instead, the lifetime of graphs in channels are compared with the bursting time in SS-HST. Vortex clusters are associated with almost the same quadrant in terms of their mean velocities during their life time, especially for those related with ejections and sweeps. As in channels, ejections in SS-HST move upwards with an average vertical velocity uτ (friction velocity) while the opposite is true for sweeps. Vortex clusters, on the other hand, are almost still in the vertical direction. In the streamwise direction, they are advected by the local mean velocity and thus deformed by the mean velocity difference. Sweeps and ejections move faster and slower than the mean velocity respectively, both by 1.5uτ . Vortex clusters move with the same speed as the mean velocity. It is verified that the incoherent structures near the wall is due to the wall instead of small size. The results suggest that coherent structures in channels are not particularly associated with the wall, or even with a given shear profile.

Energetic considerations of ciliary beating and the advantage of metachronal coordination

The internal mechanism of cilia is among the most ancient biological motors on an evolutionary scale. It produces beat patterns that consist of two phases: during the effective stroke, the cilium moves approximately as a straight rod, and during the recovery stroke, it rolls close to the surface in a tangential motion. It is commonly agreed that these two phases are designed for efficient functioning: the effective stroke encounters strong viscous resistance and generates thrust, whereas the recovery stroke returns the cilium to starting position while avoiding viscous resistance. Metachronal coordination between cilia, which occurs when many of them beat close to each other, is believed to be an autonomous result of the hydrodynamical interactions in the system. Qualitatively, metachronism is perceived as a way for reducing the energy expenditure required for beating. This paper presents a quantitative study of the energy expenditure of beating cilia, and of the energetic significance of metachronism. We develop a method for computing the work done by model cilia that beat in a viscous fluid. We demonstrate that for a single cilium, beating in water, the mechanical work done during the effective stroke is approximately five times the amount of work done during the recovery stroke. Investigation of multicilia configurations shows that having neighboring cilia beat metachronally is energetically advantageous and perhaps even crucial for multiciliary functioning. Finally, the model is used to approximate the number of dynein arm attachments that are likely to occur during the effective and recovery strokes of a beat cycle, predicting that almost all of the available dynein arms should participate in generating the motion.

Regulation of Growth Anisotropy in Well-Watered and Water-Stressed Maize Roots. II. Role of Cortical Microtubules and Cellulose Microfibrils1

We tested the hypothesis that the degree of anisotropic expansion of plant tissues is controlled by the degree of alignment of cortical microtubules or cellulose microfibrils. Previously, for the primary root of maize (Zea mays L.), we quantified spatial profiles of expansion rate in length, radius, and circumference and the degree of growth anisotropy separately for the stele and cortex, as roots became thinner with time from germination or in response to low water potential (B.M. Liang, A.M. Dennings, R.E. Sharp, T.I. Baskin [1997] Plant Physiol 115:101–111). Here, for the same material, we quantified microtubule alignment with indirect immunofluorescence microscopy and microfibril alignment throughout the cell wall with polarized-light microscopy and from the innermost cell wall layer with electron microscopy. Throughout much of the growth zone, mean orientations of microtubules and microfibrils were transverse, consistent with their parallel alignment specifying the direction of maximal expansion rate (i.e. elongation). However, where microtubule alignment became helical, microfibrils often made helices of opposite handedness, showing that parallelism between these elements was not required for helical orientations. Finally, contrary to the hypothesis, the degree of growth anisotropy was not correlated with the degree of alignment of either microtubules or microfibrils. The mechanisms plants use to specify radial and tangential expansion rates remain uncharacterized.

Influência do acabamento superficial no desempenho de lubrificantes de motor novos e usados em automóveis abastecidos com E22 e E100.

Superfícies anisotrópicas lisas e rugosas foram usadas para avaliar o efeito da rugosidade e da direção de acabamento na formação de MoS2 a partir de MoDTC em ensaios tribologicos lubrificados com óleos de motor completamente formulados. Igualmente foi avaliada a resposta de atrito de lubrificantes de motor usados em carros de passageiros e em testes de dinamômetro abastecidos com etanol (E100) e gasolina (E22). Encontrou-se que tanto a direção de acabamento quanto a rugosidade foram fundamentais na reação MoDTC - MoS2. A direção de acabamento influenciou na medida que carregamentos tangenciais geram respostas diferentes nos ensaios quando são realizados paralelos e perpendiculares às linhas de acabamento, dado que para os últimos apresenta-se maior deformação plástica das asperezas, o qual favorece a obtenção de superfícies livres de óxidos, que tem sido indicada como uma condição necessário para que aconteça a reação MoDTC - MoS2. Por esta razão os valores de coeficiente de atrito próprios da formação de MoS2 foram obtidos somente nas superfícies rugosas ensaiadas perpendiculares às marcas de acabamento. Para superfícies com valores de índice de plasticidade superiores a 1 e nos quais não são formados filmes com boas capacidades redutoras de atrito, como é o caso de ensaios realizados com óleos base (livres de aditivos), o coeficiente de atrito não depende da rugosidade e da direção de acabamento. Nos ensaios lubrificados com óleos usado, encontraram-se valores de coeficiente de atrito similares aos obtidos nas condições de lubrificação com óleo livres de aditivos, devido provavelmente à redução do MoDTC no lubrificante como tem sido identificado por diferentes autores. Quando foram comparados os óleos usados contaminados com etanol com os óleos usados contaminados com gasolina, encontrou-se maior oxidação nestes últimos. Mesmo que estas diferenças de oxidação dos óleos não significaram diferenças em termos de atrito, estas podem ser importantes na medida em que óleos mais oxidados podem favorecer o desgaste oxidativo.

Tratamiento de lixiviados de residuos de origen urbano mediante MBR

Se ha utilizado una planta de tratamiento a escala laboratorio consiste en un biorreactor de membrana (MBR). Esta planta está compuesta por un reactor biológico de 25 L de capacidad. Se utilizó una membrana plana de micro filtración marca Kubota de polietileno clorado, tamaño de poro 0,1 μm y área de filtración 0.116 m2. Se utilizaron como condiciones de operación: tiempo de residencia hidráulico 3 días, caudal de permeado 0.35 L/h y LMH 3 L/m2h. Se ha podido comprobar que es posible adaptar una población microbiológica a las particulares características químicas del lixiviado procedente de la planta y tratar estos lixiviados en un reactor biológico de membrana sumergida operando en condiciones habituales de sólidos en suspensión en el reactor entre 8-12 g/L durante un periodo de 6 meses. El proceso utilizado permite reducir la materia orgánica (97% DBO5 y 40% DQO) presente en estas corrientes residuales, agotando prácticamente toda la materia biodegradable. Respecto a los contenidos de nutrientes, el tratamiento MBR ensayado permite reducir de 35-40% el nitrógeno total, 45-50% el nitrógeno amoniacal y un 65-70% el fósforo total. Los sólidos en suspensión se han reducido en el efluente tratado en más de un 99%.

Influence of sludge retention time on filtration performance and biomass characteristics in a hollow fiber membrane bioreactor

In this study, the filtration process and the biomass characteristics in a laboratory-scale submerged membrane bioreactor (MBR) equipped with a hollow fiber (HF) microfiltration membrane were studied at different solid retention times (SRT). The MBR was fed by synthetic wastewater and the organic loading rate (OLR) was 0.5, 0.2, 0.1, and 0.08 kg COD kg VSS−1 d−1 for 10, 30, 60, and 90 days of SRT, respectively. The hydraulic retention time was 8.4 h and the permeate flux was 6 L m−2 h−1(LMH). Data analysis confirmed that at all the studied SRTs, the HF-MBR operated very good obtaining of high quality permeates. Chemical Oxygen Demand (COD) removal efficiencies were higher than 95%. The best filtration performance was reached at SRT of 30 d. On the other hand, the respirometric analysis showed that biomass was more active and there was more biomass production at low SRTs. The concentration of soluble extracellular polymeric substances (EPS) decreased with increasing SRT. A decrease of soluble EPS caused a decrease of membrane fouling rate, decreasing the frequency of chemical cleanings. The floc size decreased with SRT increasing. At high SRTs, there was more friction among particles due to the increase of the cellular density and the flocs broke decreasing their size.

Tight-Binding Approximations in 1D and 2D Coupled-Cavity Photonic Crystal Structures

Light confinement and controlling an optical field has numerous applications in the field of telecommunications for optical signals processing. When the wavelength of the electromagnetic field is on the order of the period of a photonic microstructure, the field undergoes reflection, refraction, and coherent scattering. This produces photonic bandgaps, forbidden frequency regions or spectral stop bands where light cannot exist. Dielectric perturbations that break the perfect periodicity of these structures produce what is analogous to an impurity state in the bandgap of a semiconductor. The defect modes that exist at discrete frequencies within the photonic bandgap are spatially localized about the cavity-defects in the photonic crystal. In this thesis the properties of two tight-binding approximations (TBAs) are investigated in one-dimensional and two-dimensional coupled-cavity photonic crystal structures We require an efficient and simple approach that ensures the continuity of the electromagnetic field across dielectric interfaces in complex structures. In this thesis we develop \textrm{E} -- and \textrm{D} --TBAs to calculate the modes in finite 1D and 2D two-defect coupled-cavity photonic crystal structures. In the \textrm{E} -- and \textrm{D} --TBAs we expand the coupled-cavity \overrightarrow{E} --modes in terms of the individual \overrightarrow{E} -- and \overrightarrow{D} --modes, respectively. We investigate the dependence of the defect modes, their frequencies and quality factors on the relative placement of the defects in the photonic crystal structures. We then elucidate the differences between the two TBA formulations, and describe the conditions under which these formulations may be more robust when encountering a dielectric perturbation. Our 1D analysis showed that the 1D modes were sensitive to the structure geometry. The antisymmetric \textrm{D} mode amplitudes show that the \textrm{D} --TBA did not capture the correct (tangential \overrightarrow{E} --field) boundary conditions. However, the \textrm{D} --TBA did not yield significantly poorer results compared to the \textrm{E} --TBA. Our 2D analysis reveals that the \textrm{E} -- and \textrm{D} --TBAs produced nearly identical mode profiles for every structure. Plots of the relative difference between the \textrm{E} and \textrm{D} mode amplitudes show that the \textrm{D} --TBA did capture the correct (normal \overrightarrow{E} --field) boundary conditions. We found that the 2D TBA CC mode calculations were 125-150 times faster than an FDTD calculation for the same two-defect PCS. Notwithstanding this efficiency, the appropriateness of either TBA was found to depend on the geometry of the structure and the mode(s), i.e. whether or not the mode has a large normal or tangential component.

Investigation du rôle de Myo9b dans la migration des interneurones GABAergiques corticaux

Les encéphalopathies épileptogènes sont des maladies graves de l’enfance associant une épilepsie, souvent réfractaire, et un retard de développement. Les mécanismes sous-tendant ces maladies sont peu connus. Cependant, nous postulons que ces épilepsies puissent être causées par une dysfonction du réseau inhibiteur. En effet, des défauts de migration ou de maturation des interneurones GABAergiques (INs) corticaux induisent l’épilepsie, tant chez l’humain que chez la souris. Dans le but d’étudier les causes génétiques des encéphalopathies épileptogènes sporadiques inexpliquées, le laboratoire de la Dre Rossignol a procédé au séquençage d’exome entier d’une cohorte d’enfants atteints. Cela a permis d’identifier, chez un patient, une nouvelle mutation de novo, possiblement pathogène, dans le gène MYO9b. MYO9b est impliqué dans la migration de cellules immunitaires et cancéreuses et est exprimée durant le développement cérébral. Nous émettons l’hypothèse voulant que MYO9b puisse être importante pour la migration des INs corticaux. Les résultats présentés dans ce mémoire démontrent que Myo9b est exprimé dès le stade embryonnaire par les progéniteurs des INs corticaux et que son expression se restreint aux INs dans le cortex mature. De plus, nous démontrons que la répression ex vivo de Myo9b sélectivement dans les INs au sein de tranches corticales organotypiques embryonnaires mène à des défauts morphologiques majeurs de ces cellules en migration. En effet, ces cellules présentent une morphologie multipolaire et des neurites rostraux plus longs et plus complexes. Ces changements morphologiques pourraient avoir un impact majeur sur la migration des INs et ainsi perturber le développement des réseaux inhibiteurs.

Quantifying the effect of stress relaxation on excavation

Stress relaxation is relevant to the design of both civil and mining excavations. While many authors refer to the adverse effect of stress relaxation on excavation stability, some present compelling empirical evidence indicating that stress relaxation does not have a significant effect. Establishing clear definitions of stress relaxation was critical to understanding and quantifying stress relaxation of the various types that have been referred to in the literature. This paper defines three types of stress relaxation – partial relaxation, full relaxation and tangential relaxation. Once clear definitions were determined, it became clear that the theoretical arguments and empirical evidence presented by various authors to support their respective cases are not contradictory; rather, the different conclusions can be attributed to different types of stress relaxation. In particular, when the minor principal stress is negative the intermediate principal stress has been identified as significantly affecting jointed rock mass behaviour. The aim of the study was to review and evaluate existing methods of quantifying the effect of stress relaxation around underground excavations and, if necessary, propose a new set of recommendations. An empirical stope stability model, that has been termed the Extended Mathews stability chart, was considered to be the most appropriate method of quantifying the effects of stress relaxation. A new set of guidelines to account for the effect of stress relaxation on excavation stability in the Extended Mathews stability chart has been proposed from a back-analysis of 55 case histories of stress relaxation.

Mode II delamination testing in uniaxially oriented PVC pipes

The mode II fracture toughness of an oriented PVC pipe was measured using an End Notched Flexure test geometry. A relatively low value of G(IIC) was found of 1.07 kJ m(-2) and this indicates that it is energetically more favorable for a crack to propagate in the tangential direction rather than radially through the wall of the pipe. Examination of the mechanism of crack advanced showed that although the crack was propagating globally in mode II, micro-cracks were opening ahead of the crack in mode I or in mixed mode. Growth of the crack occurred by linking up of these micro-cracks. This is similar to the mechanism found for mode II cracking in carbon fibre epoxy composites. (C) 2004 Kluwer Academic Publishers.

Comparative analysis of CNS populations in knockout mice with altered growth hormone responsiveness

Recently we have shown that growth hormone (GH) inhibits neuronal differentiation and that this process is blocked by suppressor of cytokine signalling-2 (SOCS2). Here we examine several cortical and subcortical neuronal populations in GH hyper-responsive SOCS2 null (-/-) mice and GH non-responsive GH receptor null (GHR-/-) mice. While SOCS2-/- mice showed a 30% decrease in density of NeuN positive neurons in cortex compared to wildtype, GHR-/- mice showed a 25% increase even though brain size was decreased. Interneuron sub-populations were variably affected, with a slight decrease in cortical parvalbumin expressing interneurons in SOCS2-/- mice and an increase in cortical calbindin and calretinin and striatal cholinergic neuron density in GHR-/- mice. Analysis of glial cell numbers in cresyl violet or glial fibrillary acidic protein (GFAP) stained sections of cortex showed that the neuron: glia ratio was increased in GHR-/- mice and decreased in SOCS2-/- mice. The astrocytes in GHR-/- mice appeared smaller, while they were larger in SOCS2-/- mice. Neuronal soma size also varied in the different genotypes, with smaller striatal cholinergic neurons in GHR-/- mice. While the size of layer 5 pyramidal neurons was not significantly different from wildtype, SOCS2-/- neurons were larger than GHR-/- neurons. In addition, primary dendritic length was similar in all genotypes but dendritic branching of pyramidal neurons in the cortex appeared sparser in GHR-/- and SOCS2-/- mice. These results suggest that GH, possibly regulated by SOCS2, has multiple effects on central nervous system (CNS) development and maturation, regulating the number and size of multiple neuronal and glial cell types.