Geochemistry of the sheeted dike complex at Pito Deep

Alteration of sheeted dikes exposed along submarine escarpments at the Pito Deep Rift (NE edge of the Easter microplate) provides constraints on the crustal component of axial hydrothermal systems at fast spreading mid-ocean ridges. Samples from vertical transects through the upper crust constrain the temporal and spatial scales of hydrothermal fluid flow and fluid-rock reaction. The dikes are relatively fresh (average extent of alteration is 27%), with the extent of alteration ranging from 0 to >80%. Alteration is heterogeneous on scales of tens to hundreds of meters and displays few systematic spatial trends. Background alteration is amphibole-dominated, with chlorite-rich dikes sporadically distributed throughout the dike complex, indicating that peak temperatures ranged from <300°C to >450°C and did not vary systematically with depth. Dikes locally show substantial metal mobility, with Zn and Cu depletion and Mn enrichment. Amphibole and chlorite fill fractures throughout the dike complex, whereas quartz-filled fractures and faults are only locally present. Regional variability in alteration characteristics is found on a scale of <1-2 km, illustrating the diversity of fluid-rock interaction that can be expected in fast spreading crust. We propose that much of the alteration in sheeted dike complexes develops within broad, hot upwelling zones, as the inferred conditions of alteration cannot be achieved in downwelling zones, particularly in the shallow dikes. Migration of circulating cells along rides axes and local evolution of fluid compositions produce sections of the upper crust with a distinctive character of alteration, on a scale of <1-2 km and <5-20 ka.

Chemical and mineral compositions of altered peridotites from ODP Leg 209

Abyssal peridotite from the 15°20'N area of the Mid-Atlantic Ridge show complex geochemical variations among the different sites drilled during ODP Leg 209. Major element compositions indicate variable degrees of melt depletion and refertilization as well as local hydrothermal metasomatism. Strongest evidence for melt-rock interactions are correlated Light Rare Earth Element (LREE) and High Field Strength Element (HFSE) additions at Sites 1270 and 1271. In contrast, hydrothermal alteration at Sites 1274, 1272, and 1268 causes LREE mobility associated with minor HFSE variability, reflecting the low solubility of HFSE in aqueous solutions. Site 1274 contains the least-altered, highly refractory, peridotite with strong depletion in LREE and shows a gradual increase in the intensity of isochemical serpentinization; except for the addition of H2O which causes a mass gain of up to 20 g/100 g. The formation of magnetite is reflected in decreasing Fe(2+)/Fe(3+) ratios. This style of alteration is referred to as rock-dominated serpentinization. In contrast, fluid-dominated serpentinization at Site 1268 is characterized by gains in sulfur and development of U-shaped REE pattern with strong positive Eu anomalies which are also characteristic for hot (350 to 400°C) vent-type fluids discharging from black smoker fields. Serpentinites at Site 1268 were overprinted by talc alteration under static conditions due to interaction with high a_SiO2 fluids causing the development of smooth, LREE enriched patterns with pronounced negative Eu anomalies. These results show that hydrothermal fluid-peridotite and fluid-serpentinite interaction processes are an important factor regarding the budget of exchange processes between the lithosphere and the hydrosphere in slow spreading environments.

Whole-rock geochemical analyses of serpentinized harzburgites from ODP Site 209-1270

Small-scale shear zones are present in drillcore samples of abyssal peridotites from the Mid-Atlantic ridge at 15°20'N (Ocean Drilling Program Leg 209). The shear zones act as pathways for both evolved melts and hydrothermal fluids. We examined serpentinites directly adjacent to such zones to evaluate chemical changes resulting from melt-rock and fluid-rock interaction and their influence on the mineralogy. Compared to fresh harzburgite and melt-unaffected serpentinites, serpentinites adjacent to melt-bearing veins show a marked enrichment in rare earth elements (REE), strontium and high field strength elements (HFSE) zirconium and niobium. From comparison with published chemical data of variably serpentinized and melt-unaffected harzburgites, one possible interpretation is that interaction with the adjacent melt veins caused the enrichment in HFSE, whereas the REE contents might also be enriched due to hydrothermal processes. Enrichment in alumina during serpentinization is corroborated by reaction path models for interaction of seawater with harzburgite-plagiogranite mixtures. These models explain both increased amounts of alumina in the serpentinizing fluid for increasing amounts of plagiogranitic material mixed with harzburgite, and the absence of brucite from the secondary mineralogy due to elevated silica activity. By destabilizing brucite, nearby melt veins might fundamentally influence the low-temperature alteration behaviour of serpentinites. Although observations and model results are in general agreement, due to absence of any unaltered protolith a quantification of element transport during serpentinization is not straightforward.

Chemistry of secondary minerals from the deep sheeted dike complex of ODP Holes 137-504B and 140-504B

Dolerites sampled from the lower sheeted dikes from Hole 504B during Ocean Drilling Program Legs 137 and 140, between 1562.4 and 2000.4 mbsf, were examined to document the mineralogy, petrography, and mineral parageneses associated with secondary alteration, to constrain the thermal history and composition of hydrothermal fluids. The main methods used were mineral chemical analyses by electron microprobe, X-ray diffraction, and cathodoluminescence microscopy. Temperatures of alteration were estimated on the basis of single and/or coexisting mineral chemistry. Permeability is important in controlling the type and extent of alteration in the studied dike section. At the meter-scale, intervals of weakly altered dolerites containing fresh olivine are interpreted as having experienced restricted exposure to hydrothermal fluids. At the centimeter- or millimeter-scale, alteration patches and extensively altered halos adjacent to veins reflect the permeability related to intergranular primary porosity and cracks. Most of the sheeted dike alteration in this case resulted from non-focused, pervasive fluid-rock interaction. This study confirms and extends the previous model for hydrothermal alteration at Hole 504B: hydrothermal alteration at the ridge axis followed by seawater recharge and off-axis alteration. The major new discoveries, all related to higher temperatures of alteration, are: (1) the presence of hydrothermal plagioclase (An80-95), (2) the presence of deuteric and/or hydrothermal diopside, and (3) the general increasing proportion of amphiboles, and particularly magnesio-hornblende with depth. We propose that the dolerites at Hole 504B were altered in five stages. Stage 1 occurred at high temperatures (less than 500° to 700°C) and involved late-magmatic formation of Na- and Ti-rich diopside, the hydrothermal formation of Na, Ti-poor diopside and the hydrothermal formation of an assemblage of An-rich plagioclase + hornblende. Stage 2 occurred at lower temperatures (250°-320°C) and is characterized by the appearance of actinolite, chlorite, chlorite-smectite, and/or talc (in low permeability zones) and albite. During Stage 3, quartz and epidote precipitated from evolved hydrothermal fluids at temperatures between 310° and 320°C. Anhydrite appeared during Stage 4 and likely precipitated directly from heated seawater. Stage 5 occurred off-axis at low temperatures (250°C) with laumontite and prehnite from evolved fluids.

Geochemistry of Atlantis Massif serpentinites

The Lost City hydrothermal system at the southern Atlantis Massif (Mid-Atlantic Ridge, 30°N) provides a natural laboratory for studying serpentinization processes, the temporal evolution of ultramafic-hosted hydrothermal systems, and alteration conditions during formation and emplacement of an oceanic core complex. Here we present B, O, and Sr isotope data to investigate fluid/rock interaction and mass transfer during detachment faulting and exhumation of lithospheric sequences within the Atlantis Massif. Our data indicate that extensive serpentinization was a seawater-dominated process that occurred predominately at temperatures of 150-250 °C and at high integrated W/R ratios that led to a marked boron enrichment (34-91 ppm). Boron removal from seawater during serpentinization is positively correlated with changes in d11B (11-16 per mil) but shows no correlation with O-isotope composition. Modeling indicates that B concentrations and isotope values of the serpentinites are controlled by transient temperature-pH conditions. In contrast to prior studies, we conclude that low-temperature marine weathering processes are insignificant for boron geochemistry of the Atlantis Massif serpentinites. Talc- and amphibole-rich fault rocks formed within a zone of detachment faulting at temperatures of approximately 270-350 °C and at low W/R ratios. Talc formation in ultramafic domains in the massif was subsequent to an early stage of serpentinization and was controlled by the access of Si-rich fluids derived through seawater-gabbro interactions. Replacement of serpentine by talc resulted in boron loss and significant lowering of d11B values (9-10 per mil), which we model as the product of progressive extraction of boron. Our study provides new constraints on the boron geochemical cycle at oceanic spreading ridges and suggests that serpentinization associated with ultramafic-hosted hydrothermal systems may have important implications for the behavior of boron in subduction zone settings.

Experimental and statistical investigation of canonical problems in sloshing

En hidrodinámica, el fenómeno de Sloshing se puede definir como el movimiento de la superficie libre de un fluido dentro de un contenedor sometido a fuerzas y perturbaciones externas. El fluido en cuestión experimenta violentos movimientos con importantes deformaciones de su superficie libre. La dinámica del fluido puede llegar a generar cargas hidrodinámicas considerables las cuales pueden afectar la integridad estructural y/o comprometer la estabilidad del vehículo que transporta dicho contenedor. El fenómeno de Sloshing ha sido extensivamente investigado matemática, numérica y experimentalmente, siendo el enfoque experimental el más usado debido a la complejidad del problema, para el cual los modelos matemáticos y de simulación son aun incapaces de predecir con suficiente rapidez y precisión las cargas debidas a dicho fenómeno. El flujo generado por el Sloshing usualmente se caracteriza por la presencia de un fluido multifase (gas-liquido) y turbulencia. Reducir al máximo posible la complejidad del fenómeno de Sloshing sin perder la esencia del problema es el principal reto de esta tesis doctoral, donde un trabajo experimental enfocado en casos canónicos de Sloshing es presentado y documentado con el objetivo de aumentar la comprensión de dicho fenómeno y por tanto intentar proveer información valiosa para validaciones de códigos numéricos. El fenómeno de Sloshing juega un papel importante en la industria del transporte marítimo de gas licuado (LNG). El mercado de LNG en los últimos años ha reportado un crecimiento hasta tres veces mayor al de los mercados de petróleo y gas convencionales. Ingenieros en laboratorios de investigación e ingenieros adscritos a la industria del LNG trabajan continuamente buscando soluciones económicas y seguras para contener, transferir y transportar grandes volúmenes de LNG. Los buques transportadores de LNG (LNGC) han pasado de ser unos pocos buques con capacidad de 75000 m3 hace unos treinta años, a una amplia flota con una capacidad de 140000 m3 actualmente. En creciente número, hoy día se construyen buques con capacidades que oscilan entre 175000 m3 y 250000 m3. Recientemente un nuevo concepto de buque LNG ha salido al mercado y se le conoce como FLNG. Un FLNG es un buque de gran valor añadido que solventa los problemas de extracción, licuefacción y almacenamiento del LNG, ya que cuenta con equipos de extracción y licuefacción a bordo, eliminando por tanto las tareas de transvase de las estaciones de licuefacción en tierra hacia los buques LNGC. EL LNG por tanto puede ser transferido directamente desde el FLNG hacia los buques LNGC en mar abierto. Niveles de llenado intermedios en combinación con oleaje durante las operaciones de trasvase inducen movimientos en los buques que generan por tanto el fenómeno de Sloshing dentro de los tanques de los FLNG y los LNGC. El trabajo de esta tesis doctoral lidia con algunos de los problemas del Sloshing desde un punto de vista experimental y estadístico, para ello una serie de tareas, descritas a continuación, se han llevado a cabo : 1. Un dispositivo experimental de Sloshing ha sido configurado. Dicho dispositivo ha permitido ensayar secciones rectangulares de tanques LNGC a escala con movimientos angulares de un grado de libertad. El dispositivo experimental ha sido instrumentado para realizar mediciones de movimiento, presiones, vibraciones y temperatura, así como la grabación de imágenes y videos. 2. Los impactos de olas generadas dentro de una sección rectangular de un LNGC sujeto a movimientos regulares forzados han sido estudiados mediante la caracterización del fenómeno desde un punto de vista estadístico enfocado en la repetitividad y la ergodicidad del problema. 3. El estudio de los impactos provocados por movimientos regulares ha sido extendido a un escenario más realístico mediante el uso de movimientos irregulares forzados. 4. El acoplamiento del Sloshing generado por el fluido en movimiento dentro del tanque LNGC y la disipación de la energía mecánica de un sistema no forzado de un grado de libertad (movimiento angular) sujeto a una excitación externa ha sido investigado. 5. En la última sección de esta tesis doctoral, la interacción entre el Sloshing generado dentro en una sección rectangular de un tanque LNGC sujeto a una excitación regular y un cuerpo elástico solidario al tanque ha sido estudiado. Dicho estudio corresponde a un problema de interacción fluido-estructura. Abstract In hydrodynamics, we refer to sloshing as the motion of liquids in containers subjected to external forces with large free-surface deformations. The liquid motion dynamics can generate loads which may affect the structural integrity of the container and the stability of the vehicle that carries such container. The prediction of these dynamic loads is a major challenge for engineers around the world working on the design of both the container and the vehicle. The sloshing phenomenon has been extensively investigated mathematically, numerically and experimentally. The latter has been the most fruitful so far, due to the complexity of the problem, for which the numerical and mathematical models are still incapable of accurately predicting the sloshing loads. The sloshing flows are usually characterised by the presence of multiphase interaction and turbulence. Reducing as much as possible the complexity of the sloshing problem without losing its essence is the main challenge of this phd thesis, where experimental work on selected canonical cases are presented and documented in order to better understand the phenomenon and to serve, in some cases, as an useful information for numerical validations. Liquid sloshing plays a key roll in the liquified natural gas (LNG) maritime transportation. The LNG market growth is more than three times the rated growth of the oil and traditional gas markets. Engineers working in research laboratories and companies are continuously looking for efficient and safe ways for containing, transferring and transporting the liquified gas. LNG carrying vessels (LNGC) have evolved from a few 75000 m3 vessels thirty years ago to a huge fleet of ships with a capacity of 140000 m3 nowadays and increasing number of 175000 m3 and 250000 m3 units. The concept of FLNG (Floating Liquified Natural Gas) has appeared recently. A FLNG unit is a high value-added vessel which can solve the problems of production, treatment, liquefaction and storage of the LNG because the vessel is equipped with a extraction and liquefaction facility. The LNG is transferred from the FLNG to the LNGC in open sea. The combination of partial fillings and wave induced motions may generate sloshing flows inside both the LNGC and the FLNG tanks. This work has dealt with sloshing problems from a experimental and statistical point of view. A series of tasks have been carried out: 1. A sloshing rig has been set up. It allows for testing tanks with one degree of freedom angular motion. The rig has been instrumented to measure motions, pressure and conduct video and image recording. 2. Regular motion impacts inside a rectangular section LNGC tank model have been studied, with forced motion tests, in order to characterise the phenomenon from a statistical point of view by assessing the repeatability and practical ergodicity of the problem. 3. The regular motion analysis has been extended to an irregular motion framework in order to reproduce more realistic scenarios. 4. The coupled motion of a single degree of freedom angular motion system excited by an external moment and affected by the fluid moment and the mechanical energy dissipation induced by sloshing inside the tank has been investigated. 5. The last task of the thesis has been to conduct an experimental investigation focused on the strong interaction between a sloshing flow in a rectangular section of a LNGC tank subjected to regular excitation and an elastic body clamped to the tank. It is thus a fluid structure interaction problem.

A time-domain finite element method for seakeeping and wave resistance problems

El objetivo de la tesis es la investigación de algoritmos numéricos para el desarrollo de herramientas numéricas para la simulación de problemas tanto de comportamiento en la mar como de resistencia al avance de buques y estructuras flotantes. La primera herramienta desarrollada resuelve el problema de difracción y radiación de olas. Se basan en el método de los elementos finitos (MEF) para la resolución de la ecuación de Laplace, así como en esquemas basados en MEF, integración a lo largo de líneas de corriente, y en diferencias finitas desarrollados para la condición de superficie libre. Se han desarrollado herramientas numéricas para la resolución de la dinámica de sólido rígido en sistemas multicuerpos con ligaduras. Estas herramientas han sido integradas junto con la herramienta de resolución de olas difractadas y radiadas para la resolución de problemas de interacción de cuerpos con olas. También se han diseñado algoritmos de acoplamientos con otras herramientas numéricas para la resolución de problemas multifísica. En particular, se han realizado acoplamientos con una herramienta numérica basada de cálculo de estructuras con MEF para problemas de interacción fluido-estructura, otra de cálculo de líneas de fondeo, y con una herramienta numérica de cálculo de flujos en tanques internos para problemas acoplados de comportamiento en la mar con “sloshing”. Se han realizado simulaciones numéricas para la validación y verificación de los algoritmos desarrollados, así como para el análisis de diferentes casos de estudio con aplicaciones diversas en los campos de la ingeniería naval, oceánica, y energías renovables marinas. ABSTRACT The objective of this thesis is the research on numerical algorithms to develop numerical tools to simulate seakeeping problems as well as wave resistance problems of ships and floating structures. The first tool developed is a wave diffraction-radiation solver. It is based on the finite element method (FEM) in order to solve the Laplace equation, as well as numerical schemes based on FEM, streamline integration, and finite difference method tailored for solving the free surface boundary condition. It has been developed numerical tools to solve solid body dynamics of multibody systems with body links across them. This tool has been integrated with the wave diffraction-radiation solver to solve wave-body interaction problems. Also it has been tailored coupling algorithms with other numerical tools in order to solve multi-physics problems. In particular, it has been performed coupling with a MEF structural solver to solve fluid-structure interaction problems, with a mooring solver, and with a solver capable of simulating internal flows in tanks to solve couple seakeeping-sloshing problems. Numerical simulations have been carried out to validate and verify the developed algorithms, as well as to analyze case studies in the areas of marine engineering, offshore engineering, and offshore renewable energy.

Desenvolvimento da modelagem de turbulência e interação fluido-estrutura para as vibrações induzidas por vórtices de cilindro rígido.

Esta tese apresenta o desenvolvimento e aplicação de modelos de turbulência, transição laminar-turbulenta e de interações fluido-estrutura ao escoamento externo em cilindro rígido estacionário e em vibrações induzidas por vórtices. Tais desenvolvimentos foram realizados no código ReFRESCO, baseado em técnicas de dinâmica de fluidos computacional (CFD). Realizou-se um estudo quanto ao desempenho do modelo k- SST em extensa faixa de números de Reynolds, segundo o qual se identificaram as deficiências de modelagem para este escoamento. A modelagem adaptativa das escalas (SAS) e o modelo de transição por correlações locais (LCTM), ambos combinados ao SST, melhoraram a aderência aos resultados experimentais para este escoamento, em uma contribuição original deste trabalho. A aplicação de técnicas de verificação e validação possibilitou a estimação de incertezas e erros para os modelos e números de Reynolds e também de identificada como outra contribuição deste trabalho. A combinação da modelagem em SST, SAS e LCTM com movimentos impostos de realizada para números de Reynolds moderados, diferentes frequências e amplitudes de vibração, algo que poucas publicações abordam em detalhes. Com relação aos movimentos livres, este trabalho traz contribuições com a aplicação dos modelos SST e SAS ao estudo de vibrações induzidas por vórtices em dois graus de liberdade, baixa razão de massa e números de Reynolds moderados, mais altos do que normalmente observados na literatura. Por fim, a investigação da importância relativa de efeitos da turbulência aos casos de movimentos livres e impostos, com relação ao caso de cilindro estacionário, comprovou a conjetura formulada na parte inicial deste trabalho, no que tange à escolha do modelo de turbulência em determinadas aplicações. Tal escolha mostrou-se menos decisiva no caso do cilindro em movimento imposto e ainda menos nos movimentos livres, em comparação ao caso estacionário, uma vez que a resposta em movimentos do corpo filtra grande parte dos efeitos turbulentos de ordem superior. Esta observação mostra-se relevante, uma vez que pode permitir simplificações na modelagem e aplicação de ferramentas de CFD em uma classe importante de projetos de engenharia.

Método para análise da interação fluido-estrutura em travessas do pré-distribuidor de turbinas hidráulicas.

Um dos grandes desafios enfrentados pelos fabricantes de turbinas hidráulicas é prevenir o aparecimento de vibrações induzidas pelo escoamento nas travessas do pré-distribuidor e pás do rotor. Considerando apenas as travessas, e atribuídos a tais vibrações, foram relatados 28 casos de trincas ou ruídos anormais nas últimas décadas, que acarretaram enormes prejuízos associados a reparos, atrasos e perda de geração. O estado da arte na prevenção destes problemas baseia-se na utilização de sofisticados, e caros, programas comerciais de dinâmica dos fluidos computacional para o cálculo transiente do fenômeno. Este trabalho faz uma ampla revisão bibliográfica e levantamento de eventos de trincas ou ruídos ocorridos em travessas nos últimos 50 anos. Propõe, então, um enfoque alternativo, baseado exclusivamente em ferramentas de código aberto. A partir de hipóteses simplificadoras devidamente justificadas, o problema é formulado matematicamente de forma bidimensional, no plano da seção transversal da travessa, levando em conta a interação fluido-estrutura. Nesta estratégia, as equações de Navier-Stokes são resolvidas pelo método dos elementos finitos por meio da biblioteca gratuita oomph-lib. Um código especial em C++ é desenvolvido para o problema de interação fluido-estrutura, no qual o fenômeno de turbulência é levado em consideração por meio de um algoritmo baseado no modelo de Baldwin-Lomax. O método proposto é validado por meio da comparação dos resultados obtidos com referências e medições disponíveis na literatura, que tratam de problemas de barras retangulares suportadas elasticamente. O trabalho finaliza com a aplicação do método a um estudo de caso envolvendo uma travessa particular.