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.

(Table 1) Lithium contents and isotopic composition in sediments from DSDP Holes 64-477 and 64-477A

Lithium isotopic compositions of hydrothermally altered sediments of Deep Sea Drilling Project (DSDP) site 477/477A, as well as high temperature vent fluids of the Guaymas Basin, have been determined to gain an understanding of lithium exchange during fluid-sediment interaction at this sediment-covered spreading center. Unaltered turbidite of the basin has a d6Li value of -10%, 5-7% heavier than fresh oceanic basalts. Contact metamorphism induced by a shallow sill intrusion results in a decrease of the lithium content of the adjacent sediments and a lighter isotopic value (-8%). Below the sill, sediments altered by a deep-seated hydrothermal system show strong depletions in lithium, while lithium isotopic compositions vary greatly, ranging from -11 to +1%. The shift to lighter composition is the result of preferential retention of the lighter isotope in recrystallized phases after destruction of the primary minerals. The complexity of the isotope profile is attributed to inhomogeneity in mineral composition, the tortuous pathway of fluids and the temperature effect on isotopic fractionation. The range of lithium concentration and d6Li values for the vent fluids sampled in 1982 and 1985 overlaps with that of the sediment-free mid-ocean ridge systems. The lack of a distinct expression of sediment input is explained in terms of a flow-through system with continuous water recharge. The observations on the natural system agree well with the results of laboratory hydrothermal experiments. The experimental study demonstrates the importance of temperature, pressure, water/rock ratio, substrate composition and reaction time on the lithium isotopic composition of the reacted fluid. High temperature authigenic phases do not seem to constitute an important sink for lithium and sediments of a hydrothermal system such as Guaymas are a source of lithium to the ocean. The ready mobility of lithium in the sediment under elevated temperature and pressure conditions also has important implications for lithium cycling in subduction zones.

Lithologic, mineral, chemical and isotopic compositions of altered rocks with increasing distance across the eastern flank of the Juan de Fuca Ridge, ODP Leg 168 data

During ODP Leg 168, 10 sites were drilled across the eastern flank of the Juan de Fuca Ridge (JdFR), to examine the conditions of fluid-rock interaction in three distinct hydrothermal regimes (referred to as the Hydrothermal Transition (HT), Buried Basement (BB) and Rough Basement (RB) transects), extending over a ~120 km linear transect perpendicular to the spreading ridge. This was carried out in an attempt to constrain the conditions and processes that control the location, style and magnitude of low temperature (<150°C) fluid-rock interaction within this setting. This paper presents new data on the petrology, mineral chemistry and whole rock strontium and oxygen isotopic compositions of basalts from the eastern flank of the JdFR, in order to investigate the extent, style and sequence of low-temperature hydrothermal alteration and to establish how the hydrothermal regime evolved with time. Throughout the flank, a progressive sequence of low-temperature hydrothermal alteration has been identified, marked by changes in the dominant secondary mineral assemblage, changing from: chlorite+chlorite/smectite; to iron oyxhydroxide+celadonite; to saponite+/-pyrite; culminating at present with Ca- to CaMg(+/-Fe,Mn)-carbonate. The changes in secondary mineralogy have been used to infer a series of systematic shifts in the conditions of alteration that occurred as the basement moved off-axis and was progressively buried by sediment. In general, hydrothermal alteration of the uppermost oceanic crust commenced under open, oxidative conditions, with interaction between unmodified to slightly modified seawater and basaltic crust, to a regime in which circulation of a strongly modified seawater-derived fluid was more restricted, and alteration occurred under non-oxidative conditions. Across the flank, petrological observations and microprobe analyses indicate that the observed ranges in secondary mineral composition are directly related to changes in the geochemical and textural characteristics of the basement, as well as to interaction between fluids and phases from the four stages of alteration. This is suggestive of an increase in fluid-rock increased with time. Whole rock 87Sr/86Sr and d18O analyses of basalts from across the eastern flank of the JdFR reinforce petrological observations, with 87Sr/86Sr and d18O values slightly elevated above accepted pristine MORB values for this region. These results are consistent with an increase in the amount of fluid-rock interaction with time. Across the flank, enrichment in the 87Sr/86Sr and d18O relative to MORB, is influenced by a number of factors, including: local and regional variations in the crustal lithology and structure; the age of the crust; the extent of bulk rock alteration; and theoretically, the relative abundance of different isotopically-enriched secondary mineral phases in the crust.

Effects of potential models on the adsorption of ethane and ethylene on graphitized thermal carbon black. Study of two-dimensional critical temperature and isosteric heat versus loading

Adsorption of ethylene and ethane on graphitized thermal carbon black and in slit pores whose walls are composed of graphene layers is studied in detail to investigate the packing efficiency, the two-dimensional critical temperature, and the variation of the isosteric heat of adsorption with loading and temperature. Here we used a Monte Carlo simulation method with a grand canonical Monte Carlo ensemble. A number of two-center Lennard-Jones (LJ) potential models are investigated to study the impact of the choice of potential models in the description of adsorption behavior. We chose two 2C-LJ potential models in our investigation of the (i) UA-TraPPE-LJ model of Martin and Siepmann (J. Phys. Chem. B 1998,102, 25692577) for ethane and Wick et al. (J. Phys. Chem. B 2000,104, 8008-8016) for ethylene and (ii) AUA4-LJ model of Ungerer et al. (J. Chem. Phys. 2000,112, 5499-5510) for ethane and Bourasseau et al. (J. Chem. Phys. 2003, 118, 3020-3034) for ethylene. These models are used to study the adsorption of ethane and ethylene on graphitized thermal carbon black. It is found that the solid-fluid binary interaction parameter is a function of adsorbate and temperature, and the adsorption isotherms and heat of adsorption are well described by both the UA-TraPPE and AUA models, although the UA-TraPPE model performs slightly better. However, the local distributions predicted by these two models are slightly different. These two models are used to explore the two-dimensional condensation for the graphitized thermal carbon black, and these values are 110 K for ethylene and 120 K for ethane.