Chemical composition of sediments from DSDP Hole 149

The diagenesis and geochemical evolution of deep-sea sediments are controlled by the interaction between sediments and their associated pore waters. With increasing depth, the pore water of Hole 149 (DSDP) exhibits a strong depletion in Mg and a corresponding enrichment in Ca, while the alkalinity remains relatively constant. Dissolved SiO2 is nearly constant in the upper 100 m of sediment, but is highly enriched in the deepest pore waters. The pore waters exhibit a depletion in K with increasing depth, and O18/O16 pore water ratios also decrease. The sediment section has three zones of sedimentary regimes with increasing depth in the drill hole: an upper 100 m section of detrital clays, a middle section enriched in calc-akalic volcanics which have undergone submarine weathering to a smectite phase, and a lower section of siliceous ooze which still has a diagenetic smectite phase. The quartz-feldspar ratios and O18/O16 composition of the silicate phases are in agreement with these interpretations. The submarine weathering of volcanics to a smectite can account for the observed pore water gradients. Volcanics release Ca and Mg to the pore waters causing the alkalinity values to increase. Smectite is formed, depletes the pore waters in Mg and O18 and causes the alkalinity to decrease. The net reaction allows for the observed relationship between pore water Ca and Mg gradients with little net change in alkalinity. Given the abundance of volcanics in many deep-sea sediments, especially in lower sections which often form near ridge crests, the submarine formation of smectite may be an additional oceanic Mg sink which has not yet been fully considered.

Chemical, isotopic and mineral compositions of basal metalliferous sediments from DSDP Legs 5 and 7

Sediments from near the basement of a number of Deep Sea Drilling Project (DSDP) sites, from the Bauer Deep, and from the East Pacific Rise have unusually high transition metal-to-aluminum ratios. Similarities in the chemical, isotopic, and mineralogical compositions of these deposits point to a common origin. All the sediments studied have rare-earth-element (REE) patterns strongly resembling the pattern of sea water, implying either that the REE's were coprecipitated with ferromanganese hydroxyoxides (hydroxyoxides denote a mixture of unspecified hydrated oxides and hydroxides), or that they are incorporated in small concentrations of phosphatic fish debris found in all samples. Oxygen isotopic data indicate that the metalliferous sediments are in isotopic equilibrium with sea water and are composed of varying mixtures of two end-member phases with different oxygen isotopic compositions: an iron-manganese hydroxyoxide and an iron-rich montmorillonite. A low-temperature origin for the sediments is supported by mineralogical analyses by x-ray diffraction which show that goethite, iron-rich montmorillonite, and various manganese hydroxyoxides are the dominant phases present. Sr87/Sr86 ratios for the DSDP sediments are indistinguishable from the Sr87/Sr86 ratio in modern sea water. Since these sediments were formed 30 to 90 m.y. ago, when sea water had a lower Sr87/Sr86 value, the strontium in the poorly crystalline hydroxyoxides must be exchanging with interstitial water in open contact with sea water. In contrast, uranium isotopic data indicate that the metalliferous sediments have formed a closed system for this element. The sulfur isotopic compositions suggest that sea-water sulfur dominates these sediments with little or no contribution of magmatic or bacteriologically reduced sulfur. In contrast, ratios of lead isotopes in the metalliferous deposits resemble values for oceanic tholeiite basalt, but are quite different from ratios found in authigenic marine manganese nodules. Thus, lead in the metalliferous sediments appears to be of magmatic origin. The combined mineralogical, isotopic, and chemical data for these sediments suggest that they formed from hydrothermal solutions generated by the interaction of sea water with newly formed basalt crust at mid-ocean ridges. The crystallization of solid phases took place at low temperatures and was strongly influenced by sea water, which was the source for some of the elements found in the sediments.

Decline in growth of foraminifer Marginopora rossi under eutrophication and ocean acidification scenarios

The combination of global and local stressors is leading to a decline in coral reef health globally. In the case of eutrophication, increased concentrations of dissolved inorganic nitrogen (DIN) and phosphorus (DIP) are largely attributed to local land use changes. From the global perspective, increased atmospheric CO2 levels are not only contributing to global warming but also ocean acidification (OA). Both eutrophication and OA have serious implications for calcium carbonate production and dissolution among calcifying organisms. In particular, benthic foraminifera precipitate the most soluble form of mineral calcium carbonate (high-Mg calcite), potentially making them more sensitive to dissolution. In this study, a manipulative orthogonal two-factor experiment was conducted to test the effects of dissolved inorganic nutrients and OA on the growth, respiration and photophysiology of the large photosymbiont-bearing benthic foraminifer, Marginopora rossi. This study found the growth rate of M. rossi was inhibited by the interaction of eutrophication and acidification. The relationship between M. rossi and its photosymbionts became destabilized due to the photosymbiont's release from nutrient limitation in the nitrate-enriched treatment, as shown by an increase in zooxanthellae cells per host surface area. Foraminifers from the OA treatments had an increased amount of Chl a per cell, suggesting a greater potential to harvest light energy, however, there was no net benefit to the foraminifer growth. Overall, this study demonstrates that the impacts of OA and eutrophication are dose dependent and interactive. This research indicates an OA threshold at pH 7.6, alone or in combination with eutrophication, will lead to a decline in M. rossi calcification. The decline in foraminifera calcification associated with pollution and OA will have broad ecological implications across their ubiquitous range and suggests that without mitigation it could have serious implications for the future of coral reefs.

Numerical simulation of tangling in jet engine turbines

The numerical analysis of certain safety related problems presents serious difficulties, since the large number of components present leads to huge finite elementmodels that can only be solved by using large and expensive computers or by making rough approaches to the problem. Tangling, or clashing, in the turbine of a jet engine airplane is an example of such problems. This is caused by the crash and friction between rotor and stator blades in the turbine after an eventual shaft failure. When facing the study of an event through numerical modelling, the accurate simulation of this problem would require the engineer to model all the rotor and stator blades existing in the turbine stage, using a small element size in all pieces. Given that the number of stator and rotor blades is usually around 200, such simulations would require millions of elements. This work presents a new numerical methodology, specifically developed for the accurate modelling of the tangling problem that, depending on the turbine configuration, is able to reduce the number of nodes up to an order of magnitude without losing accuracy. The methodology, which benefits from the cyclic configuration of turbines, is successfully applied to the numerical analysis of a hypothetical tangling event in a turbine, providing valuable data such as the rotating velocity decrease of the turbine, the braking torque and the damage suffered by the blades. The methodology is somewhat general and can be applied to any problem in which damage caused by the interaction between a rotating and static piece is to be analysed.

Consolidation problems

The analysis of deformation in soils is of paramount importance in geotechnical engineering. For a long time the complex behaviour of natural deposits defied the ingenuity of engineers. The time has come that, with the aid of computers, numerical methods will allow the solution of every problem if the material law can be specified with a certain accuracy. Boundary Techniques (B.E.) have recently exploded in a splendid flowering of methods and applications that compare advantegeously with other well-established procedures like the finite element method (F.E.). Its application to soil mechanics problems (Brebbia 1981) has started and will grow in the future. This paper tries to present a simple formulation to a classical problem. In fact, there is already a large amount of application of B.E. to diffusion problems (Rizzo et al, Shaw, Chang et al, Combescure et al, Wrobel et al, Roures et al, Onishi et al) and very recently the first specific application to consolidation problems has been published by Bnishi et al. Here we develop an alternative formulation to that presented in the last reference. Fundamentally the idea is to introduce a finite difference discretization in the time domain in order to use the fundamental solution of a Helmholtz type equation governing the neutral pressure distribution. Although this procedure seems to have been unappreciated in the previous technical literature it is nevertheless effective and straightforward to implement. Indeed for the special problem in study it is perfectly suited, because a step by step interaction between the elastic and flow problems is needed. It allows also the introduction of non-linear elastic properties and time dependent conditions very easily as will be shown and compares well with performances of other approaches.

Convective instabilities in two superposed horizontal liquid layers heated laterally.

This work is devoted to the theoretical study of the stability of two superposed horizontal liquid layers bounded by two solid planes and subjected to a horizontal temperature gradient. The liquids are supposed to be immiscible with a nondeformable interface. The forces acting on the system are buoyancy and interfacial tension. Four different flow patterns and temperature profiles are found for the basic state. A linear perturbative analysis with respect to two- and three-dimensional perturbations reveals the existence of three kinds of patterns. Depending on the relative height of both liquids several situations are predicted: either wave propa- gation from cold to the hot regions, or waves propagating in the opposite direction or still stationary longitu- dinal rolls. The behavior of three different pairs of liquids which have been used in experiments on bilayers under vertical gradient by other authors have been examined. The instability mechanisms are discussed and a qualitative interpretation of the different behaviors exhibited by the system is provided. In some configurations it is possible to find a codimension-two point created by the interaction of two Hopf modes with different frequencies and wave numbers. These results suggest to consider two liquid layers as an interesting prototype ? nard-Marangoni problem.

Inertial behavior of offshore devices

In the last years, there has been a continued growth in the number of offshore operations for handling large equipment and objects, with emphasis on installation and maintenance of devices for exploiting marine renewable energy like generators for harnessing wind, waves and currents energy. Considering the behaviour of these devices during manoeuvrings, and due to their size and by the interaction with the surrounding fluid, the effect of inertial forces and torques is very important, which requires a specific modelling. This paper especially discusses the masses and moments of inertia modelling problem, with the aim to use it in the simulation of the complex manoeuvres of these devices and in the automatic control systems designed for their offshore operations. Given the importance and complexity of the added mass modelling, a method for its early design identification, developed by the R&D Group on Marine Renewable Energy Technologies of the UPM (GITERM) and its use on special cases like emersion manoeuvres of devices from underwater to the surface will be presented.

Raman spectroscopy in Group IV nanowires and nanowire axial heterostructures

The control of the SiGe NW composition is fundamental for the fabrication of high quality heterostructures. Raman spectroscopy has been used to analyse the composition of SiGe alloys. We present a study of the Raman spectrum of SiGe nanowires and SiGe/Si heterostructures. The inhomogeneity of the Ge composition deduced from the Raman spectrum is explained by the existence of a Ge-rich outer shell and by the interaction of the NW with the electromagnetic field associated with the laser beam.

Flamelet structures in spray ignition

In typical liquid-fueled burners the fuel is injected as a high-velocity liquid jet that breaks up to form the spray. The initial heating and vaporization of the liquid fuel rely on the relatively large temperatures of the sourrounding gas, which may include hot combustion products and preheated air. The heat exchange between the liquid and the gas phases is enhanced by droplet dispersion arising from the turbulent motion. Chemical reaction takes place once molecular mixing between the fuel vapor and the oxidizer has occurred in mixing layers separating the spray flow from the hot air stream. Since in most applications the injection velocities are much larger than the premixed-flame propagation velocity, combustion stabilization relies on autoignition of the fuel-oxygen mixture, with the combustion stand-off distance being controlled by the interaction of turbulent transport, droplet heating and vaporization, and gas-phase chemical reactions. In this study, conditions are identified under which analyses of laminar flamelets canshed light on aspects of turbulent spray ignition. This study extends earlier fundamental work by Liñan & Crespo (1976) on ignition in gaseous mixing layers to ignition of sprays. Studies of laminar mixing layers have been found to be instrumental in developing un-derstanding of turbulent combustion (Peters 2000), including the ignition of turbulent gaseous diffusion flames (Mastorakos 2009). For the spray problem at hand, the configuration selected, shown in Figure 1, involves a coflow mixing layer formed between a stream of hot air moving at velocity UA and a monodisperse spray moving at velocity USUA. The boundary-layer approximation will be used below to describe the resulting sl ender flow, which exhibits different igniting behaviors depending on the characteristics of t he fuel. In this approximation, consideration of the case U A = U S enables laminar ignition distances to be related to ignition times of unstrained spray flamelets, thereby pro viding quantitative information of direct applicability in regions of low scala r dissipation-rate in turbulent reactive flows (see the discussion in pp. 181–186 of Peters (2000)) . This report is organized as follows. Effects of droplet dispersion dynamics on ignition of sprays in turbulent mixing layers are discussed in Section 2. The formulation f or ignition in laminar mixing layers is outlined in Sections 3 and 4. The results are presented in Section 5. In Section 6, the mixture-fraction field and associated scalar dissipat ion rates for spray ignition are discussed. Finally, some brief conclusions are drawn in Section 7.

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.