Gas hydrate investigation at the Amsterdam mud volcano

We investigated gas hydrate in situ inventories as well as the composition and principal transport mechanisms of fluids expelled at the Amsterdam mud volcano (AMV; 2,025 m water depth) in the Eastern Mediterranean Sea. Pressure coring (the only technique preventing hydrates from decomposition during recovery) was used for the quantification of light hydrocarbons in near-surface deposits. The cores (up to 2.5 m in length) were retrieved with an autoclave piston corer, and served for analyses of gas quantities and compositions, and pore-water chemistry. For comparison, gravity cores from sites at the summit and beyond the AMV were analyzed. A prevalence of thermogenic light hydrocarbons was inferred from average C1/C2+ ratios <35 and d13C-CH4 values of -50.6 per mil. Gas venting from the seafloor indicated methane oversaturation, and volumetric gas-sediment ratios of up to 17.0 in pressure cores taken from the center demonstrated hydrate presence at the time of sampling. Relative enrichments in ethane, propane, and iso-butane in gas released from pressure cores, and from an intact hydrate piece compared to venting gas suggest incipient crystallization of hydrate structure II (sII). Nonetheless, the co-existence of sI hydrate can not be excluded from our dataset. Hydrates fill up to 16.7% of pore volume within the sediment interval between the base of the sulfate zone and the maximum sampling depth at the summit. The concave-down shapes of pore-water concentration profiles recorded in the center indicate the influence of upward-directed advection of low-salinity fluids/fluidized mud. Furthermore, the SO42- and Ba2+ pore-water profiles in the central part of the AMV demonstrate that sulfate reduction driven by the anaerobic oxidation of methane is complete at depths between 30 cm and 70 cm below seafloor. Our results indicate that methane oversaturation, high hydrostatic pressure, and elevated pore-water activity caused by low salinity promote fixing of considerable proportions of light hydrocarbons in shallow hydrates even at the summit of the AMV, and possibly also of other MVs in the region. Depending on their crystallographic structure, however, hydrates will already decompose and release hydrocarbon masses if sediment temperatures exceed ca. 19.3°C and 21.0°C, respectively. Based on observations from other mud volcanoes, the common occurrence of such temperatures induced by heat flux from below into the immediate subsurface appears likely for the AMV.

Geochemistry of basalts from the Caribbean Large Igneous Province

Oceanic flood basalts are poorly understood, short-term expressions of highly increased heat flux and mass flow within the convecting mantle. The uniqueness of the Caribbean Large Igneous Province (CLIP, 92-74 Ma) with respect to other Cretaceous oceanic plateaus is its extensive sub-aerial exposures, providing an excellent basis to investigate the temporal and compositional relationships within a starting plume head. We present major element, trace element and initial Sr-Nd-Pb isotope composition of 40 extrusive rocks from the Caribbean Plateau, including onland sections in Costa Rica, Colombia and Curaçao as well as DSDP Sites in the Central Caribbean. Even though the lavas were erupted over an area of ~3*10**6 km**2, the majority have strikingly uniform incompatible element patterns (La/Yb=0.96+/-0.16, n=64 out of 79 samples, 2sigma) and initial Nd-Pb isotopic compositions (e.g. 143Nd/144Ndin=0.51291+/-3, epsilon-Nd i=7.3+/-0.6, 206Pb/204Pbin=18.86+/-0.12, n=54 out of 66, 2sigma). Lavas with endmember compositions have only been sampled at the DSDP Sites, Gorgona Island (Colombia) and the 65-60 Ma accreted Quepos and Osa igneous complexes (Costa Rica) of the subsequent hotspot track. Despite the relatively uniform composition of most lavas, linear correlations exist between isotope ratios and between isotope and highly incompatible trace element ratios. The Sr-Nd-Pb isotope and trace element signatures of the chemically enriched lavas are compatible with derivation from recycled oceanic crust, while the depleted lavas are derived from a highly residual source. This source could represent either oceanic lithospheric mantle left after ocean crust formation or gabbros with interlayered ultramafic cumulates of the lower oceanic crust. High 3He/4He in olivines of enriched picrites at Quepos are ~12 times higher than the atmospheric ratio suggesting that the enriched component may have once resided in the lower mantle. Evaluation of the Sm-Nd and U-Pb isotope systematics on isochron diagrams suggests that the age of separation of enriched and depleted components from the depleted MORB source mantle could have been <=500 Ma before CLIP formation and interpreted to reflect the recycling time of the CLIP source. Mantle plume heads may provide a mechanism for transporting large volumes of possibly young recycled oceanic lithosphere residing in the lower mantle back into the shallow MORB source mantle.

Amino acid racemization for several DSDP samples (Table 1)

Kinetic parameters for the epimerization of isoleucine in multispecific foraminiferal asemblages were used to establish the effects of burial depth and the geothermal gradient on the extent of reaction. It was observed that with a little as thirty meters of burial in a normal thermal regime there were differences between the extent of epimerization measured and that which would have been predicted for thermal equilibrium with bottom water temperatures. As would be expected, these differences are greatest when the heat flow (the geothermal gradient) and/or the sedimentation rates are highest. These effects were observed in most of the DSDP samples studied, and have been used to estimate the average heat flux since the time of sample deposition. Occasional anomalous effects were observed which could not be related to past or present heat flux. These were determined to be due to such geologic occurrences as slumping and reworking or to recent sample contamination. Other problems emerged related to bottom water temperatures including changes over geologic time which are unknown and could not be deduced. Thus, the presence of epimerization anomalies in DSDP cores as noted above limits the effectiveness of amino acid geochronology in such cores, unless these anomalies can be recognized as ab initio.

Polymer Fluid Dynamics: Continuum and Molecular Appoaches

To solve problems in polymer fluid dynamics, one needs the equation of continuity, motion, and energy. The last two equations contain the stress tensor and the heat-flux vector for the material. There are two ways to formulate the stress tensor: (1) one can write a continuum expression for the stress tensor in terms of kinematic tensors, or (2) one can select a molecular model that represents the polymer molecule, and then develop an expression for the stress tensor from kinetic theory. The advantage of the kinetic theory approach is that one gets information about the relation between the molecular structure of the polymers and the rheological properties. In this review, we restrict the discussion primarily to the simplest stress tensor expressions or “constitutive equations” containing from two to four adjustable parameters, although we do indicate how these formulations may be extended to give more complicated expressions. We also explore how these simplest expressions are recovered as special cases of a more general framework, the Oldroyd 8-constant model. The virtue of studying the simplest models is that we can discover some general notions as to which types of empiricisms or which types of molecular models seem to be worth investigating further. We also explore equivalences between continuum and molecular approaches. We restrict the discussion to several types of simple flows, such as shearing flows and extensional flows. These are the flows that are of greatest importance in industrial operations. Furthermore, if these simple flows cannot be well described by continuum or molecular models, then it is not necessary to lavish time and energy to apply them to more complex flow problems.

Temperature fields in a channel partially filled with a porous material under local thermal non-equilibrium condition-an exact solution

This work examines analytically the forced convection in a channel partially filled with a porous material and subjected to constant wall heat flux. The Darcy–Brinkman–Forchheimer model is used to represent the fluid transport through the porous material. The local thermal non-equilibrium, two-equation model is further employed as the solid and fluid heat transport equations. Two fundamental models (models A and B) represent the thermal boundary conditions at the interface between the porous medium and the clear region. The governing equations of the problem are manipulated, and for each interface model, exact solutions, for the solid and fluid temperature fields, are developed. These solutions incorporate the porous material thickness, Biot number, fluid to solid thermal conductivity ratio and Darcy number as parameters. The results can be readily used to validate numerical simulations. They are, further, applicable to the analysis of enhanced heat transfer, using porous materials, in heat exchangers.

Multiphase modeling of melting / solidification with high density variations using XFEM

La modélisation de la cryolite, utilisée dans la fabrication de l’aluminium, implique plusieurs défis, notament la présence de discontinuités dans la solution et l’inclusion de la difference de densité entre les phases solide et liquide. Pour surmonter ces défis, plusieurs éléments novateurs ont été développés dans cette thèse. En premier lieu, le problème du changement de phase, communément appelé problème de Stefan, a été résolu en deux dimensions en utilisant la méthode des éléments finis étendue. Une formulation utilisant un multiplicateur de Lagrange stable spécialement développée et une interpolation enrichie a été utilisée pour imposer la température de fusion à l’interface. La vitesse de l’interface est déterminée par le saut dans le flux de chaleur à travers l’interface et a été calculée en utilisant la solution du multiplicateur de Lagrange. En second lieu, les effets convectifs ont été inclus par la résolution des équations de Stokes dans la phase liquide en utilisant la méthode des éléments finis étendue aussi. Troisièmement, le changement de densité entre les phases solide et liquide, généralement négligé dans la littérature, a été pris en compte par l’ajout d’une condition aux limites de vitesse non nulle à l’interface solide-liquide pour respecter la conservation de la masse dans le système. Des problèmes analytiques et numériques ont été résolus pour valider les divers composants du modèle et le système d’équations couplés. Les solutions aux problèmes numériques ont été comparées aux solutions obtenues avec l’algorithme de déplacement de maillage de Comsol. Ces comparaisons démontrent que le modèle par éléments finis étendue reproduit correctement le problème de changement phase avec densités variables.

Prediction of Upward Flame Spread over Polymers

In this work, the existing understanding of flame spread dynamics is enhanced through an extensive study of the heat transfer from flames spreading vertically upwards across 5 cm wide, 20 cm tall samples of extruded Poly (Methyl Methacrylate) (PMMA). These experiments have provided highly spatially resolved measurements of flame to surface heat flux and material burning rate at the critical length scale of interest, with a level of accuracy and detail unmatched by previous empirical or computational studies. Using these measurements, a wall flame model was developed that describes a flame’s heat feedback profile (both in the continuous flame region and the thermal plume above) solely as a function of material burning rate. Additional experiments were conducted to measure flame heat flux and sample mass loss rate as flames spread vertically upwards over the surface of seven other commonly used polymers, two of which are glass reinforced composite materials. Using these measurements, our wall flame model has been generalized such that it can predict heat feedback from flames supported by a wide range of materials. For the seven materials tested here – which present a varied range of burning behaviors including dripping, polymer melt flow, sample burnout, and heavy soot formation – model-predicted flame heat flux has been shown to match experimental measurements (taken across the full length of the flame) with an average accuracy of 3.9 kW m-2 (approximately 10 – 15 % of peak measured flame heat flux). This flame model has since been coupled with a powerful solid phase pyrolysis solver, ThermaKin2D, which computes the transient rate of gaseous fuel production of constituents of a pyrolyzing solid in response to an external heat flux, based on fundamental physical and chemical properties. Together, this unified model captures the two fundamental controlling mechanisms of upward flame spread – gas phase flame heat transfer and solid phase material degradation. This has enabled simulations of flame spread dynamics with a reasonable computational cost and accuracy beyond that of current models. This unified model of material degradation provides the framework to quantitatively study material burning behavior in response to a wide range of common fire scenarios.