On the role of global flow instability analysis in closed loop flow control

Control of linear flow instabilities has been demonstrated to be an effective theoretical flow control methodology, capable of modifying transitional flow on canonical geometries such as the plane channel and the flat-plate boundary layer.

Selected topics in the Hydrodynamics of floating offshore wind turbines: heave damping and second order surge motions

La Energía eléctrica producida mediante tecnología eólica flotante es uno de los recursos más prometedores para reducir la dependencia de energía proveniente de combustibles fósiles. Esta tecnología es de especial interés en países como España, donde la plataforma continental es estrecha y existen pocas áreas para el desarrollo de estructuras fijas. Entre los diferentes conceptos flotantes, esta tesis se ha ocupado de la tipología semisumergible. Estas plataformas pueden experimentar movimientos resonantes en largada y arfada. En largada, dado que el periodo de resonancia es largo estos puede ser inducidos por efectos de segundo orden de deriva lenta que pueden tener una influencia muy significativa en las cargas en los fondeos. En arfada las fuerzas de primer orden pueden inducir grandes movimientos y por tanto la correcta determinación del amortiguamiento es esencial para la analizar la operatividad de la plataforma. Esta tesis ha investigado estos dos efectos, para ello se ha usado como caso base el diseño de una plataforma desarrollada en el proyecto Europeo Hiprwind. La plataforma se compone de 3 columnas cilíndricas unidas mediante montantes estructurales horizontales y diagonales, Los cilindros proporcionan flotabilidad y momentos adrizante. A la base de cada columna se le ha añadido un gran “Heave Plate” o placa de cierre. El diseño es similar a otros diseños previos (Windfloat). Se ha fabricado un modelo a escala de una de las columnas para el estudio detallado del amortiguamiento mediante oscilaciones forzadas. Las dimensiones del modelo (1m diámetro en la placa de cierre) lo hacen, de los conocidos por el candidato, el mayor para el que se han publicado datos. El diseño del cilindro se ha realizado de tal manera que permite la fijación de placas de cierre planas o con refuerzo, ambos modelos se han fabricado y analizado. El modelo con refuerzos es una reproducción exacta del diseño a escala real incluyendo detalles distintivos del mismo, siendo el más importante la placa vertical perimetral. Los ensayos de oscilaciones forzadas se han realizado para un rango de frecuencias, tanto para el disco plano como el reforzado. Se han medido las fuerzas durante los ensayos y se han calculado los coeficientes de amortiguamiento y de masa añadida. Estos coeficientes son necesarios para el cálculo del fondeo mediante simulaciones en el dominio del tiempo. Los coeficientes calculados se han comparado con la literatura existente, con cálculos potenciales y por ultimo con cálculos CFD. Para disponer de información relevante para el diseño estructural de la plataforma se han medido y analizado experimentalmente las presiones en la parte superior e inferior de cada placa de cierre. Para la correcta estimación numérica de las fuerzas de deriva lenta en la plataforma se ha realizado una campaña experimental que incluye ensayos con modelo cautivo de la plataforma completa en olas bicromaticas. Pese a que estos experimentos no reproducen un escenario de oleaje realista, los mismos permiten una verificación del modelo numérico mediante la comparación de fuerzas medidas en el modelo físico y el numérico. Como resultados de esta tesis podemos enumerar las siguientes conclusiones. 1. El amortiguamiento y la masa añadida muestran una pequeña dependencia con la frecuencia pero una gran dependencia con la amplitud del movimiento. siendo coherente con investigaciones existentes. 2. Las medidas con la placa de cierre reforzada con cierre vertical en el borde, muestra un amortiguamiento significativamente menor comparada con la placa plana. Esto implica que para ensayos de canal es necesario incluir estos detalles en el modelo. 3. La masa añadida no muestra grandes variaciones comparando placa plana y placa con refuerzos. 4. Un coeficiente de amortiguamiento del 6% del crítico se puede considerar conservador para el cálculo en el dominio de la frecuencia. Este amortiguamiento es equivalente a un coeficiente de “drag” de 4 en elementos de Morison cuadráticos en las placas de cierre usadas en simulaciones en el dominio del tiempo. 5. Se han encontrado discrepancias en algunos valores de masa añadida y amortiguamiento de la placa plana al comparar con datos publicados. Se han propuesto algunas explicaciones basadas en las diferencias en la relación de espesores, en la distancia a la superficie libre y también relacionadas con efectos de escala. 6. La presión en la placa con refuerzos son similares a las de la placa plana, excepto en la zona del borde donde la placa con refuerzo vertical induce una gran diferencias de presiones entre la cara superior e inferior. 7. La máxima diferencia de presión escala coherentemente con la fuerza equivalente a la aceleración de la masa añadida distribuida sobre la placa. 8. Las masas añadidas calculadas con el código potencial (WADAM) no son suficientemente precisas, Este software no contempla el modelado de placas de pequeño espesor con dipolos, la poca precisión de los resultados aumenta la importancia de este tipo de elementos al realizar simulaciones con códigos potenciales para este tipo de plataformas que incluyen elementos de poco espesor. 9. Respecto al código CFD (Ansys CFX) la precisión de los cálculos es razonable para la placa plana, esta precisión disminuye para la placa con refuerzo vertical en el borde, como era de esperar dado la mayor complejidad del flujo. 10. Respecto al segundo orden, los resultados, en general, muestran que, aunque la tendencia en las fuerzas de segundo orden se captura bien con los códigos numéricos, se observan algunas reducciones en comparación con los datos experimentales. Las diferencias entre simulaciones y datos experimentales son mayores al usar la aproximación de Newman, que usa únicamente resultados de primer orden para el cálculo de las fuerzas de deriva media. 11. Es importante remarcar que las tendencias observadas en los resultados con modelo fijo cambiarn cuando el modelo este libre, el impacto que los errores en las estimaciones de fuerzas segundo orden tienen en el sistema de fondeo dependen de las condiciones ambientales que imponen las cargas ultimas en dichas líneas. En cualquier caso los resultados que se han obtenido en esta investigación confirman que es necesaria y deseable una detallada investigación de los métodos usados en la estimación de las fuerzas no lineales en las turbinas flotantes para que pueda servir de guía en futuros diseños de estos sistemas. Finalmente, el candidato espera que esta investigación pueda beneficiar a la industria eólica offshore en mejorar el diseño hidrodinámico del concepto semisumergible. ABSTRACT Electrical power obtained from floating offshore wind turbines is one of the promising resources which can reduce the fossil fuel energy consumption and cover worldwide energy demands. The concept is the most competitive in countries, such as Spain, where the continental shelf is narrow and does not provide space for fixed structures. Among the different floating structures concepts, this thesis has dealt with the semisubmersible one. Platforms of this kind may experience resonant motions both in surge and heave directions. In surge, since the platform natural period is long, such resonance can be excited with second order slow drift forces and may have substantial influence on mooring loads. In heave, first order forces can induce significant motion, whose damping is a crucial factor for the platform downtime. These two topics have been investigated in this thesis. To this aim, a design developed during HiPRWind EU project, has been selected as reference case study. The platform is composed of three cylindrical legs, linked together by a set of structural braces. The cylinders provide buoyancy and restoring forces and moments. Large circular heave plates have been attached to their bases. The design is similar to other documented in literature (e.g. Windfloat), which implies outcomes could have a general value. A large scale model of one of the legs has been built in order to study heave damping through forced oscillations. The final dimensions of the specimen (one meter diameter discs) make it, to the candidate’s knowledge, the largest for which data has been published. The model design allows for the fitting of either a plain solid heave plate or a flapped reinforced one; both have been built. The latter is a model scale reproduction of the prototype heave plate and includes some distinctive features, the most important being the inclusion of a vertical flap on its perimeter. The forced oscillation tests have been conducted for a range of frequencies and amplitudes, with both the solid plain model and the vertical flap one. Forces have been measured, from which added mass and damping coefficients have been obtained. These are necessary to accurately compute time-domain simulations of mooring design. The coefficients have been compared with literature, and potential flow and CFD predictions. In order to provide information for the structural design of the platform, pressure measurements on the top and bottom side of the heave discs have been recorded and pressure differences analyzed. In addition, in order to conduct a detailed investigation on the numerical estimations of the slow-drift forces of the HiPRWind platform, an experimental campaign involving captive (fixed) model tests of a model of the whole platform in bichromatic waves has been carried out. Although not reproducing the more realistic scenario, these tests allowed a preliminary verification of the numerical model based directly on the forces measured on the structure. The following outcomes can be enumerated: 1. Damping and added mass coefficients show, on one hand, a small dependence with frequency and, on the other hand, a large dependence with the motion amplitude, which is coherent with previously published research. 2. Measurements with the prototype plate, equipped with the vertical flap, show that damping drops significantly when comparing this to the plain one. This implies that, for tank tests of the whole floater and turbine, the prototype plate, equipped with the flap, should be incorporated to the model. 3. Added mass values do not suffer large alterations when comparing the plain plate and the one equipped with a vertical flap. 4. A conservative damping coefficient equal to 6% of the critical damping can be considered adequate for the prototype heave plate for frequency domain analysis. A corresponding drag coefficient equal to 4.0 can be used in time domain simulations to define Morison elements. 5. When comparing to published data, some discrepancies in added mass and damping coefficients for the solid plain plate have been found. Explanations have been suggested, focusing mainly on differences in thickness ratio and distance to the free surface, and eventual scale effects. 6. Pressures on the plate equipped with the vertical flap are similar in magnitude to those of the plain plate, even though substantial differences are present close to the edge, where the flap induces a larger pressure difference in the reinforced case. 7. The maximum pressure difference scales coherently with the force equivalent to the acceleration of the added mass, distributed over the disc surface. 8. Added mass coefficient values predicted with the potential solver (WADAM) are not accurate enough. The used solver does not contemplate modeling thin plates with doublets. The relatively low accuracy of the results highlights the importance of these elements when performing potential flow simulations of offshore platforms which include thin plates. 9. For the full CFD solver (Ansys CFX), the accuracy of the computations is found reasonable for the plain plate. Such accuracy diminishes for the disc equipped with a vertical flap, an expected result considering the greater complexity of the flow. 10. In regards to second order effects, in general, the results showed that, although the main trend in the behavior of the second-order forces is well captured by the numerical predictions, some under prediction of the experimental values is visible. The gap between experimental and numerical results is more pronounced when Newman’s approximation is considered, making use exclusively of the mean drift forces calculated in the first-order solution. 11. It should be observed that the trends observed in the fixed model test may change when the body is free to float, and the impact that eventual errors in the estimation of the second-order forces may have on the mooring system depends on the characteristics of the sea conditions that will ultimately impose the maximum loads on the mooring lines. Nevertheless, the preliminary results obtained in this research do confirm that a more detailed investigation of the methods adopted for the estimation of the nonlinear wave forces on the FOWT would be welcome and may provide some further guidance for the design of such systems. As a final remark, the candidate hopes this research can benefit the offshore wind industry in improving the hydrodynamic design of the semi-submersible concept.

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.

Application of the homogenization techniques to the determination of the effective flexure constants of plate structural systems

A Mindlin plate with periodically distributed ribs patterns is analyzed by using homogenization techniques based on asymptotic expansion methods. The stiffness matrix of the homogenized plate is found to be dependent on the geometrical characteristics of the periodical cell, i.e. its skewness, plan shape, thickness variation etc. and on the plate material elastic constants. The computation of this plate stiffness matrix is carried out by averaging over the cell domain some solutions of different periodical boundary value problems. These boundary value problems are defined in variational form by linear first order differential operators on the cell domain and the boundary conditions of the variational equation correspond to a periodic structural problem. The elements of the stiffness matrix of homogenized plate are obtained by linear combinations of the averaged solution functions of the above mentioned boundary value problems. Finally, an illustrative example of application of this homogenization technique to hollowed plates and plate structures with ribs patterns regularly arranged over its area is shown. The possibility of using in the profesional practice the present procedure to the actual analysis of floors of typical buildings is also emphasized.

DNA analysis and diagnostics on oligonucleotide microchips.

We present a further development in the technology of sequencing by hybridization to oligonucleotide microchips (SHOM) and its application to diagnostics for genetic diseases. A robot has been constructed to manufacture sequencing "microchips." The microchip is an array of oligonucleotides immobilized into gel elements fixed on a glass plate. Hybridization of the microchip with fluorescently labeled DNA was monitored in real time simultaneously for all microchip elements with a two-wavelength fluorescent microscope equipped with a charge-coupled device camera. SHOM has been used to detect beta-thalassemia mutations in patients by hybridizing PCR-amplified DNA with the microchips. A contiguous stacking hybridization technique has been applied for the detection of mutations; it can simplify medical diagnostics and enhance its reliability. The use of multicolor monitoring of contiguous stacking hybridization is suggested for large-scale diagnostics and gene polymorphism studies. Other applications of the SHOM technology are discussed.

Tuning the functional properties of Silicon Carbide Thin Film and Nanowires

The present thesis has been devoted to the synthesis and investigation of functional properties of silicon carbide thin films and nanowires. The work took profit from the experience of the research group in the synthesis of 3C-SiC from vapour phase. 3C-SiC thin films Thin films heteroepitaxy on silicon substrates was carried out in a vapour phase epitaxy reactor. The initial efforts were committed to the process development in order to enhance the crystal quality of the epi-layer. The carbonization process and a buffer layer procedure were optimized in order to obtain good quality monocrystalline 3C-SiC layers. The films characterization was used not only to improve the entire process, but also to assess the crystalline quality and to identify the defects. Methyltrichlorosilane (MTS) was introduced during the synthesis to increase the growth rate and enhance crystalline quality. The effect of synthesis parameters such as MTS flow and process temperature was studied in order to promote defect density reduction and the release of the strain due to lattice mismatch between 3C-SiC and silicon substrate. In-growth n-type doping was implemented using a nitrogen gas line and the effect of different synthesis parameters on doping level was studied. Raman measurements allowed a contactless characterization and evaluation of electrically active dopant. The effect of MTS on nitrogen incorporation was investigated and a promotion of dopant concentration together with a higher growth rate were demonstrated. This result allows to obtain higher doping concentrations without deteriorating crystal quality in 3C-SiC and, to the best of our knowledge, it has never been demonstrated before. 3C-SiC nanowires Core-shell SiC-SiO2 nanowires were synthesized using a chemical vapour deposition technique in an open tube configuration reactor on silicon substrates. Metal catalyst were used to promote a uniaxial growth and a dense bundle of nanowires 100 µm long and 60 nm thick was obtained. Substrate preparation was found to be fundamental in order to obtain a uniform nanowire density. Morphological characterization was carried out using scanning electron microscopy and the analysis of structural, compositional, optical properties is reported.

O elemento finito T6-3i na análise de placas e dinâmica de cascas.

O método dos elementos finitos é o método numérico mais difundido na análise de estruturas. Ao longo das últimas décadas foram formulados inúmeros elementos finitos para análise de cascas e placas. As formulações de elementos finitos lidam bem com o campo de deslocamentos, mas geralmente faltam testes que possam validar os resultados obtidos para o campo das tensões. Este trabalho analisa o elemento finito T6-3i, um elemento finito triangular de seis nós proposto dentro de uma formulação geometricamente exata, em relação aos seus resultados de tensões, comparando-os com as teorias analíticas de placas, resultados de tabelas para o cálculo de momentos em placas retangulares e do ANSYSr, um software comercial para análise estrutural, mostrando que o T6-3i pode apresentar resultados insatisfatórios. Na segunda parte deste trabalho, as potencialidades do T6-3i são expandidas, sendo proposta uma formulação dinâmica para análise não linear de cascas. Utiliza-se um modelo Lagrangiano atualizado e a forma fraca é obtida do Teorema dos Trabalhos Virtuais. São feitas simulações numéricas da deformação de domos finos que apresentam vários snap-throughs e snap-backs, incluindo domos com vincos curvos, mostrando a robustez, simplicidade e versatilidade do elemento na sua formulação e na geração das malhas não estruturadas necessárias para as simulações.

Particel size distribution and microstructures of sediments from ODP Leg 205 sites (Table 1)

Ocean Drilling Program Legs 170 and 205 offshore Costa Rica provide structural observations which support a new model for the geometry and deformation response to the seismic cycle of the frontal sedimentary prism and decollement. The model is based on drillcore, thin section, and electron microscope observations. The decollement damage zone is a few tens of meters in width, it develops mainly within the frontal prism. A clear cm-thick fault core is observed 1.6 km from the trench. The lower boundary of the fault core is coincident with the lithological boundary between the frontal prism and the hemipelagic and pelagic sediment of the Cocos plate. Breccia clast distributions in the upper portion of the decollement damage zone were studied through fractal analysis. This analysis shows that the fractal dimension changes with brecciated fragment size, implying that deformation was not accommodated by self-similar fracturing. A higher fractal dimensionality correlates with smaller particle size, which indicates that different or additional grain-size reduction processes operated during shearing. The co-existence of two distinct fracturing processes is also confirmed by microscopic analysis in which extension fracturing in the upper part of the damage zone farthest from the fault core is frequent, while both extension and shear fracturing occur approaching the fault core. The coexistence of extensional and shear fracturing seems to be best explained by fluid pressure variations in response to variations of the compressional regime during the seismic cycle. During the co-seismic event, sub-horizontal compression and fluid pressure increase, triggering shear fracturing and fluid expulsion. Fractures migrate upward with fluids, contributing to the asymmetric shape of the decollement, while slip propagates. In the inter-seismic interval the frontal prismrelaxes and fluid pressure drops. The frontal prismgoes into diffuse extension during the intervalwhen plate convergence is accommodated by creep along the ductile fault core. The fault core is typically a barrier to deformation, which is explained by its weak, but impermeable, nature. The localized development of a damage zone beneath the fault core is characterized by shear fracturing that appears as the result of local strengthening of the detachment.

Rock eval analysis of organic matter in sediment core 159-959D

The Cretaceous Equatorial Atlantic Gateway between the Central and South Atlantic basins is of interest not only for paleoceanographic and paleoclimatic studies, but also because it provided particularly favourable conditions for the accumulation and preservation of organic-rich sediments. Deposition of carbonaceous sediments along the Côte d'Ivoire-Ghana Transform Margin (Ocean Drilling Program Leg 159) was intimately linked to the plate tectonic and paleoceanographic evolution of this gateway. Notably, the formation of a marginal basement ridge on the southeastern border of the transform margin provided an efficient shelter of the landward Deep Ivorian Basin against erosive and potentially oxidizing currents. Different subsidence histories across the transform margin were responsible for the development of distinct depositional settings on the crest and on both sides of the basement ridge. Whereas the southern, oceanward flank of the basement ridge was characterized by rapid, continuous deepening since last Albian-early Cenomanian, marine sedimentation on the northern, landward flank was interrupted by a period of uplift and erosion in the late Albian, and rapid subsidence started after the early Coniacian. Organic-rich sediments occur throughout almost the entire Cretaceous section, but hydrogen-rich marine black shales were exclusively recovered from core sections above an uplift-related unconformity. These black shales formed when separation of Africa and South America was sufficient to allow permanent oceanic midwater exchange after the late Albian. Four periods of black shale accumulation are recovered, some of them are correlated with the global oceanic anoxic events: in the last Albian-earliest Cenomanian, at the Cenomanian-Turronian boundary, during the middle Coniacian-early Campanian, and in the mid-Maastrichtian. These periods were characterized by increasing carbon flux to the seafloor, induced by enhanced palaeoproductivity and intensified supply of terrestrial organic matter. Black shale depostion appears to be intimately linked to periods of rising or maximum eustatic sea level and to the expansion of the oxygen minimum zone, as indicated by foraminiferal biofacies. Intervals between black shales units, in contrast, indicate a shrinking oxygen minimum zone and enhanced detrital flux rates, probably related to lowering sea level. Upper Cretaceous detritral limestones with high porosities may provide excellent hydrocarbon reservoirs, alsthough their areal extent appears to be limited. Palaeogene porcellanites, capped by Neogene pelagic marls and clays, extend over a wider area and max provide another target for hydrocarbon exploration.

(Table 1) Results from XRD analysis as well as sediment type and selected results from ring shear experiments

The location of the seaward tip of a subduction thrust controls material transfer at convergent plate margins, and hence global mass balances. At approximately half of those margins, the material of the subducting plate is completely underthrust so that no accretion or even subduction erosion takes place. Along the remaining margins, material is scraped off the subducting plate and added to the upper plate by frontal accretion. We here examine the physical properties of subducting sediments off Costa Rica and Nankai, type examples for an erosional and an accretionary margin, to investigate which parameters control the level where the frontal thrust cuts into the incoming sediment pile. A series of rotary-shear experiments to measure the frictional strength of the various lithologies entering the two subduction zones were carried out. Results include the following findings: (1) At Costa Rica, clay-rich strata at the top of the incoming succession have the lowest strength (µres = 0.19) while underlying calcareous ooze, chalk and diatomite are strong (up to µres = 0.43; µpeak = 0.56). Hence the entire sediment package is underthrust. (2) Off Japan, clay-rich deposits within the lower Shikoku Basin inventory are weakest (µres = 0.13-0.19) and favour the frontal proto-thrust to migrate into one particular horizon between sandy, competent turbidites below and ash-bearing mud above. (3) Taking in situ data and earlier geotechnical testing into account, it is suggested that mineralogical composition rather than pore-pressure defines the position of the frontal thrust, which locates in the weakest, clay mineral-rich (up to 85 wt.%) materials. (4) Smectite, the dominant clay mineral phase at either margin, shows rate strengthening and stable sliding in the frontal 50 km of the subduction thrust (0.0001-0.1 mm/s, 0.5-25 MPa effective normal stress). (5) Progressive illitization of smectite cannot explain seismogenesis, because illite-rich samples also show velocity strengthening at the conditions tested.

Grain size distribution analysis of prominent beds on Southern Central Chilean seamounts

Gravity cores obtained from isolated seamounts located within, and rising up to 300 m from the sediment-filled Peru-Chile Trench off Southern Central Chile (36°S-39°S) contain numerous turbidite layers which are much coarser than the hemipelagic background sedimentation. The mineralogical composition of some of the beds indicates a mixed origin from various source terrains while the faunal assemblage of benthic foraminifera in one of the turbidite layers shows a mixed origin from upper shelfal to middle-lower bathyal depths which could indicate a multi-source origin and therefore indicate an earthquake triggering of the causing turbidity currents. The bathymetric setting and the grain size distribution of the sampled layers, together with swath echosounder and sediment echosounder data which monitor the distribution of turbidites on the elevated Nazca Plate allow some estimates on the flow direction, flow velocity and height of the causing turbidity currents. We discuss two alternative models of deposition, both of which imply high (175-450 m) turbidity currents and we suggest a channelized transport process as the general mode of turbidite deposition. Whether these turbidites are suspension fallout products of thick turbiditic flows or bedload deposits from sheet-like turbidity currents overwhelming elevated structures cannot be decided upon using our sedimentological data, but the specific morphology of the seamounts rather argues for the first option. Oxygen isotope stratigraphy of one of the cores indicates that the turbiditic sequences were deposited during the last Glacial period and during the following transition period and turbiditic deposition stopped during the Holocene. This climatic coupling seems to be dominant, while the occurrence of megathrust earthquakes provides a trigger mechanism. This seismic triggering takes effect only during times of very high sediment supply to the shelf and slope.

Beauty illustrated by an analysis and classification of beauty in woman : with a critical view of the hypotheses of Hume, Hogarth, Burke, Knight, Alison, etc., and the hypotheses of beauty in sculpture and painting, by Leonardo da Vinci, Winckelmann, Mengs, Bossi, etc. /

The plates include one frontispiece, 22 plates numbered I-XXII and one addtional plate numbered XVIb.

Non-linear analysis of the thermo-electro-mechanical behaviour of shear deformable FGM plates with piezoelectric actuators

This paper investigates the non-linear bending behaviour of functionally graded plates that are bonded with piezoelectric actuator layers and subjected to transverse loads and a temperature gradient based on Reddy's higher-order shear deformation plate theory. The von Karman-type geometric non-linearity, piezoelectric and thermal effects are included in mathematical formulations. The temperature change is due to a steady-state heat conduction through the plate thickness. The material properties are assumed to be graded in the thickness direction according to a power-law distribution in terms of the volume fractions of the constituents. The plate is clamped at two opposite edges, while the remaining edges can be free, simply supported or clamped. Differential quadrature approximation in the X-axis is employed to convert the partial differential governing equations and the associated boundary conditions into a set of ordinary differential equations. By choosing the appropriate functions as the displacement and stress functions on each nodal line and then applying the Galerkin procedure, a system of non-linear algebraic equations is obtained, from which the non-linear bending response of the plate is determined through a Picard iteration scheme. Numerical results for zirconia/aluminium rectangular plates are given in dimensionless graphical form. The effects of the applied actuator voltage, the volume fraction exponent, the temperature gradient, as well as the characteristics of the boundary conditions are also studied in detail. Copyright (C) 2004 John Wiley Sons, Ltd.

Transient bending of a piezoelectric circular plate

Thin, piezoelectric circular plates are frequently used as active components in transducer and smart materials applications. This paper reports on the exact, explicit solution for the transient motion of a piezoelectric circular plate, built-in or simply supported on the edge and electrically grounded over the entire surface. Expressed by elementary Bessel functions and obtained via exact inverse Laplace transforms, the solution enables the efficient calculation of accurate system parameters. (C) 2004 Elsevier Ltd. All rights reserved.

Formation of mesostructured titania thin films using iospropoxide precursors

Mesostructured titania thin films were prepared by an evaporation-induced self-assembly process. The highly acidic sot precursors contained titanium(IV) tetraisopropoxide (TTIP) as a titanium source, a tri-block copolymer Pluronic P123 as a template, and acetylacetonate and HCl as hydrolysis inhibitors. Characteristics of the resultant titania thin films were studied using X-ray diffraction (XRD) analysis, N-2-adsorption/desorption analysis, and transmission electron microscopy (TEM). XRD and TEM investigations on the as-synthesised films revealed the appearance of cubic-like, pseudohexagonal, and lamellar mesophases; depending on the amount of water in the sols of film precursors. Template removal by a calcination process yields high surface area (320-360 m(2)/g) mesoporous materials with crystalline anatase frameworks. Water content also influences the degree of anatase crystallinity of the calcined films. Higher water content resulted in improved anatase crystallinity. These nanostructured materials are of interest for photocatalysts, pbotoelectrochemical solar cells and other photonic devices. (C) 2003 Elsevier B.V. All rights reserved.