Fault-controlled and stratabound dolostones in the Late Aptian-earliest Albian Benassal Formation (Maestrat Basin, E Spain) : petrology and geochemistry constrains

This study was developed under the ExxonMobil FC2 Alliance (Fundamental Controls on Flow in Carbonates). The authors wish to thank ExxonMobil Production Company and ExxonMobil Upstream Research Company for providing funding. The views in this article by Sherry L. Stafford are her own and not necessarily those of ExxonMobil. This research was supported by the Sedimentary Geology Research Group of the Generalitat de Catalunya (2014SGR251). We would like to thank Andrea Ceriani and Paola Ronchi for their critical and valuable reviews, and Associated Editor Piero Gianolla for the editorial work.

Biogeochemical investigation of sediments at the Hakon Mosby Mud Volcano (Barents Sea)

We have investigated if in a cold seep methane or sulfide is used for chemosynthetic primary production and if significant amounts of the sulfide produced by anaerobic oxidation of methane are oxidized geochemically and hence are not available for chemosynthetic production. Geochemically controlled redox reactions and biological turnover were compared in different habitats of the Håkon Mosby Mud Volcano. The center of the mud volcano is characterized by the highest fluid flow, and most primary production by the microbial community depends on oxidation of methane. The small amount of sulfide produced is oxidized geochemically with oxygen or is precipitated with dissolved iron. In the medium flow peripheral Beggiatoa habitat sulfide is largely oxidized biologically. The oxygen and nitrate supply is high enough that Beggiatoa can oxidize the sulfide completely, and chemical sulfide oxidation or precipitation is not important. An internally stored nitrate reservoir with average concentrations of 110 mmol L-1 enables the Beggiatoa to oxidize sulfide anaerobically. The pH profile indicates sequential sulfide oxidation with elemental sulfur as intermediate. Gray thiotrophic mats associated with perturbed sediments showed a high heterogeneity in sulfate turnover and high sulfide fluxes, balanced by the opposing oxygen and nitrate fluxes so that biological oxidation dominates over geochemical sulfide removal processes. The three habitats indicate substantial small-scale variability in carbon fixation pathways either through direct biological use of methane or through indirect carbon fixation of methane-derived carbon dioxide by chemolithotrophic sulfide oxidation.

Phospholipid fatty acid and quinone composition of sediments from ODP Leg 205 sites

Assessing the habitability of deep-sea sediments undergoing compaction, compression, and subduction at convergent margins adds to our understanding of the limits of the terrestrial biosphere. In this work, we report exploratory biomarker data on sediments obtained at Ocean Drilling Program (ODP) Sites 1253, 1254, and 1255 during drilling at the Costa Rica subduction trench and forearc sedimentary wedge. The samples selected for postcruise biomarker analyses were located within intervals of potentially enhanced fluid flow within the décollement and sedimentary wedge fault zones (Sites 1254 and 1255) and within basal carbonates at the reference site (Site 1253). The passage of fluids that are geochemically distinct from ambient interstitial water provides a disequilibrium setting that may enhance habitability. Biomarker data show low levels of microbial biomass in subseafloor sediments sampled at the Costa Rica convergent margin as deep as ~370 meters below seafloor.

(Table 1) Boron concentration and d11B values of interstitial water from IODP Holes 308-U1319A, U1320A, U1322B and U1324B

The Integrated Ocean Drilling Program Expedition 308 (IODP308) drilled normal-pressured sediments from the Brazos-Trinity Basin IV and over-pressured sediments from the Ursa Basin on the northern slope of the Gulf of Mexico. The interstitial water samples from the normal-pressured basin show B concentrations and B isotopic compositions ranging from 255 to 631 µM (0.6 to 1.5 times of seawater value) and from +29.1 to +42.7 per mil (relative to NIST SRM 951), respectively. A wider range is observed both for B concentrations (292 to 865 µM, 0.7 to 2.1 times of seawater value) and d11B values (+25.5 to +43.2 per mil) of the interstitial water in the over-pressured basin. The down-core distribution of B concentrations and d11B values in the interstitial waters are sensitive tracers for assessing various processes occurring in the sediment column, including boron adsorption/desorption reactions involving clay minerals and organic matter in sediments as well as fluid migration and mixing in certain horizons and in the sediment column. In the normal-pressured basin adsorption/desorption reactions in shallow sediments play the major role in controlling the B content and B isotopic composition of the interstitial water. In contrast, multiple processes affect the B content and d11B of the interstitial water in the over-pressured Ursa Basin. There, the stratigraphic level of the maxima of B and d11B correspond to seismic reflectors. The intruded fluids along the seismic reflector boundary from high to low-topography mix with local interstitial water. Fluid flow is inferred in the Blue Unit (a coarse sandstone layer, connecting the high- to low-pressured region) from the freshening of interstitial water in Ursa Basin Site U1322, and upward flow by the overpressure expels fluid from the overburden above the Blue Unit.

(Table 2) Carbonate chemistry and isotopes from northern Cascadia margin samples

Most authigenic carbonates previously recovered from the Cascadia slope have 87Sr/86Sr signatures that reflect shallow precipitation in equilibrium with coeval seawater. There is also evidence for carbonate formation supported by fluids that have been modified by reactions with the incoming Juan de Fuca plate (87Sr/86Sr = 0.7071; Teichert et al., 2005, doi:10.1016/j.epsl.2005.08.002) or with terrigenous turbidites (87Sr/86Sr = 0.70975 to 0.71279; Sample et al., 1993, doi:10.1130/0091-7613(1993)021<0507:CCICFF>2.3.CO;2). We report on the strontium isotopic composition of carbonates and fluids from IODP Site U1329 and nearby Barkley Canyon (offshore Vancouver Island), which have strontium isotope ratios as low as 0.70539. Whereas the strontium and oxygen isotopic compositions of carbonates from paleoseeps in the uplifted Coast Range forearc indicate formation in ambient bottom seawater, several samples from the Pysht/Sooke Fm. show a 87Sr-depleted signal (87Sr/86Sr = 0.70494 and 0.70511) similar to that of the anomalous Site U1329 and Barkley Canyon carbonates. Our data, when analyzed in the context of published elemental and isotopic composition of these carbonates (Joseph et al., 2012, doi:10.1016/j.palaeo.2013.01.012 ), point to two formation mechanisms: 1) shallow precipitation driven by the anaerobic oxidation of methane (AOM) with d13C values as low as -50 per mil and contemporaneous 87Sr/86Sr seawater ratios, and 2) carbonate precipitation driven by fluids that have circulated through the oceanic crust, which are depleted in 87Sr. Carbonates formed from the second mechanism precipitate both at depth and at sites of deep-sourced fluid seepage on the seafloor. The 87Sr-depleted carbonates and pore fluids found at Barkley Canyon represent migration of a deep, exotic fluid similar to that found in high permeability conglomerate layers at 188 mbsf of Site U1329, and which may have fed paleoseeps in the Pysht/Sooke Fm. These exotic fluids likely reflect interaction with the 52-57 Ma igneous Crescent Terrane, which supplies fluids with high calcium, manganese and strontium enriched in the non-radiogenic nucleide. Tectonic compression and dehydration reactions then force these fluids updip, where they pick up the thermogenic hydrocarbons and 13C-enriched dissolved inorganic carbon that are manifested in fluids and carbonates sampled at Barkley Canyon and at Site U1329. The Crescent Terrane may have sourced cold seeps in this margin since at least the late Oligocene.

On the feasibility of time-lapse superconducting gravimetry for reservoir monitoring

The feasibility of monitoring fluid flow subsurface processes that result in density changes, using the iGrav superconducting gravimeter, is investigated. Practical targets include steam-assisted gravity drainage (SAGD) bitumen depletion and water pumping from aquifers, for which there is currently a void in low-impact, inexpensive monitoring techniques. This study demonstrates that the iGrav has the potential to be applied to multi-scale and diverse reservoirs. Gravity and gravity gradient signals are forward modeled for a real SAGD reservoir at two time steps, and for surface-fed and groundwater-fed aquifer pumping models, to estimate signal strength and directional dependency of water flow. Time-lapse gravimetry on small-scale reservoirs exhibits two obstacles, namely, a µgal sensitivity requirement and high noise levels in the vicinity of the reservoir. In this study, both limitations are overcome by proposing (i) a portable superconducting gravimeter, and (ii) a pair of instruments under various baseline geometries. This results in improved spatial resolution for locating depletion zones, as well as the cancellation of noise common in both instruments. Results indicate that a pair of iGrav superconducting gravimeters meet the sensitivity requirements and the spatial focusing desired to monitor SAGD bitumen migration at the reservoir scales. For SAGD reservoirs, the well pair separation, reservoir depth, and survey sampling determine the resolvability of individual well pair depletion patterns during the steam chamber rising phase, and general reservoir depletion patterns during the steam chamber spreading phase. Results show that monitoring water table elevation changes due to pumping and tracking whether groundwater or surface water is being extracted are feasible.

Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico

LINS, Filipe C. A. et al. Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico. In: CONGRESSO BRASILEIRO DE PESQUISA E DESENVOLVIMENTO EM PETRÓLEO E GÁS, 5. 2009, Fortaleza, CE. Anais... Fortaleza: CBPDPetro, 2009.

Modelling the mixing of lead bullion during the refining process

Abstract not available

Investigating wetting characteristics on microstructured surfaces for superhydrophobicity and metal microcasting

The engineering of liquid behavior on surfaces is important for infrastructure, transportation, manufacturing, and sensing. Surfaces can be rendered superhydrophobic by microstructuring, and superhydrophobic devices could lead to practical corrosion inhibition, self-cleaning, fluid flow control, and surface drag reduction. To more fully understand how liquid interacts with microstructured surfaces, this dissertation introduces a direct method for determining droplet solid-liquid-vapor interfacial geometry on microstructured surfaces. The technique performs metrology on molten metal droplets deposited onto microstructured surfaces and then frozen. Unlike other techniques, this visualization technique can be used on large areas of curved and opaque microstructured surfaces to determine contact line. This dissertation also presents measurements and models for how curvature and flexing of microstructured polymers affects hydrophobicity. Increasing curvature of microstructured surfaces leads to decreased slide angle for liquid droplets suspended on the surface asperities. For a surface with regularly spaced asperities, as curvature becomes more positive, droplets suspended on the tops of asperities are suspended on fewer asperities. Curvature affects superhydrophobicity because microscopic curvature changes solid-liquid interaction, pitch is altered, and curvature changes the shape of the three phase contact line. This dissertation presents a model of droplet interactions with curved microstructured surfaces that can be used to design microstructure geometries that maintain the suspension of a droplet when curved surfaces are covered with microstructured polymers. Controlling droplet dynamics could improve microfluidic devices and the shedding of liquids from expensive equipment, preventing corrosion and detrimental performance. This dissertation demonstrates redirection of dynamic droplet spray with anisotropic microstructures. Superhydrophobic microstructured surfaces can be economically fabricated using metal embossing masters, so this dissertation describes casting-based microfabrication of metal microstructures and nanostructures. Low melting temperature metal was cast into flexible silicone molds which were themselves cast from microfabricated silicon templates. The flexibility of the silicone mold permits casting of curved surfaces, which this dissertation demonstrates by fabricating a cylindrical metal roller with microstructures. The metal microstructures can be in turn used as a reusable molding tool. This dissertation also describes an industrial investment casting process to produce aluminum molds having integrated microstructures. Unlike conventional micromolding tools, the aluminum mold was large and had complex curved surfaces. The aluminum was cast into curved microstructured ceramic molds which were themselves cast from curved microstructured rubber. Many structures were successfully cast into the aluminum with excellent replication fidelity, including circular, square, and triangular holes. This dissertation demonstrates molding of large, curved surfaces having surface microstructures using the aluminum mold. This work contributes a more full understanding of the phenomenon of superhydrophobicity and techniques for the economic fabrication of superhydrophobic microstructures.

Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico

LINS, Filipe C. A. et al. Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico. In: CONGRESSO BRASILEIRO DE PESQUISA E DESENVOLVIMENTO EM PETRÓLEO E GÁS, 5. 2009, Fortaleza, CE. Anais... Fortaleza: CBPDPetro, 2009.

Studies of turbulent flows in continuous casting of steel with and without magnetic field

This thesis develops and tests various transient and steady-state computational models such as direct numerical simulation (DNS), large eddy simulation (LES), filtered unsteady Reynolds-averaged Navier-Stokes (URANS) and steady Reynolds-averaged Navier-Stokes (RANS) with and without magnetic field to investigate turbulent flows in canonical as well as in the nozzle and mold geometries of the continuous casting process. The direct numerical simulations are first performed in channel, square and 2:1 aspect rectangular ducts to investigate the effect of magnetic field on turbulent flows. The rectangular duct is a more practical geometry for continuous casting nozzle and mold and has the option of applying magnetic field either perpendicular to broader side or shorter side. This work forms the part of a graphic processing unit (GPU) based CFD code (CU-FLOW) development for magnetohydrodynamic (MHD) turbulent flows. The DNS results revealed interesting effects of the magnetic field and its orientation on primary, secondary flows (instantaneous and mean), Reynolds stresses, turbulent kinetic energy (TKE) budgets, momentum budgets and frictional losses, besides providing DNS database for two-wall bounded square and rectangular duct MHD turbulent flows. Further, the low- and high-Reynolds number RANS models (k-ε and Reynolds stress models) are developed and tested with DNS databases for channel and square duct flows with and without magnetic field. The MHD sink terms in k- and ε-equations are implemented as proposed by Kenjereš and Hanjalić using a user defined function (UDF) in FLUENT. This work revealed varying accuracies of different RANS models at different levels. This work is useful for industry to understand the accuracies of these models, including continuous casting. After realizing the accuracy and computational cost of RANS models, the steady-state k-ε model is then combined with the particle image velocimetry (PIV) and impeller probe velocity measurements in a 1/3rd scale water model to study the flow quality coming out of the well- and mountain-bottom nozzles and the effect of stopper-rod misalignment on fluid flow. The mountain-bottom nozzle was found more prone to the longtime asymmetries and higher surface velocities. The left misalignment of stopper gave higher surface velocity on the right leading to significantly large number of vortices forming behind the nozzle on the left. Later, the transient and steady-state models such as LES, filtered URANS and steady RANS models are combined with ultrasonic Doppler velocimetry (UDV) measurements in a GaInSn model of typical continuous casting process. LES-CU-LOW is the fastest and the most accurate model owing to much finer mesh and a smaller timestep. This work provided a good understanding on the performance of these models. The behavior of instantaneous flows, Reynolds stresses and proper orthogonal decomposition (POD) analysis quantified the nozzle bottom swirl and its importance on the turbulent flow in the mold. Afterwards, the aforementioned work in GaInSn model is extended with electromagnetic braking (EMBr) to help optimize a ruler-type brake and its location for the continuous casting process. The magnetic field suppressed turbulence and promoted vortical structures with their axis aligned with the magnetic field suggesting tendency towards 2-d turbulence. The stronger magnetic field at the nozzle well and around the jet region created large scale and lower frequency flow behavior by suppressing nozzle bottom swirl and its front-back alternation. Based on this work, it is advised to avoid stronger magnetic field around jet and nozzle bottom to get more stable and less defect prone flow.

Chaotic mixing of a competitive-consecutive reaction

The evolution of a competitive-consecutive chemical reaction is computed numerically in a two-dimensional chaotic fluid flow with initially segregated reactants. Results from numerical simulations are used to evaluate a variety of reduced models commonly adopted to model the full advection-reaction-diffusion problem. Particular emphasis is placed upon fast reactions, where the yield varies most significantly with Peclet number (the ratio of diffusive to advective time scales). When effects of the fluid mechanical mixing are strongest, we find that the yield of the reaction is underestimated by a one-dimensional lamellar model that ignores the effects of fluid mixing, but overestimated by two other lamellar models that include fluid mixing.

Three-dimensional kinetic simulations of active suspensions: effect of chemotaxis and steric interaction

In this work, the effects of chemotaxis and steric interactions in active suspensions are analyzed by extending the kinetic model proposed by Saintillan and Shelley [1, 2]. In this model, a conservation equation for the active particle configuration is coupled to the Stokes equation for the flow arising from the force dipole exerted by the particles on the fluid. The fluid flow equations are solved spectrally and the conservation equation is solved by second-order finite differencing in space and second-order Adams-Bashforth time marching. First, the dynamics in suspensions of oxytactic run-and-tumble bacteria confined in thin liquid films surrounded by air is investigated. These bacteria modify their tumbling behavior by making temporal comparisons of the oxygen concentration, and, on average, swim towards high concentrations of oxygen. The kinetic model proposed by Saintillan and Shelley [1, 2] is modified to include run-and-tumble effects and oxygentaxis. The spatio-temporal dynamics of the oxygen and bacterial concentration are analyzed. For small film thicknesses, there is a weak migration of bacteria to the boundaries, and the oxygen concentration is high inside the film as a result of diffusion; both bacterial and oxygen concentrations quickly reach steady states. Above a critical film thickness (approximately 200 micron), a transition to chaotic dynamics is observed and is characterized by turbulent-like 3D motion, the formation of bacterial plumes, enhanced oxygen mixing and transport into the film, and hydrodynamic velocities of magnitudes up to 7 times the single bacterial swimming speed. The simulations demonstrate that the combined effects of hydrodynamic interactions and oxygentaxis create collective three-dimensional instabilities which enhances oxygen availability for the bacteria. Our simulation results are consistent with the experimental findings of Sokolov et al. [3], who also observed a similar transition with increasing film thickness. Next, the dynamics in concentrated suspensions of active self-propelled particles in a 3D periodic domain are analyzed. We modify the kinetic model of Saintillan and Shelley [1, 2] by including an additional nematic alignment torque proportional to the local concentration in the equation for the rotational velocity of the particles, causing them to align locally with their neighbors (Doi and Edwards [4]). Large-scale three- dimensional simulations show that, in the presence of such a torque both pusher and puller suspensions are unstable to random fluctuations and are characterized by highly nematic structures. Detailed measures are defined to quantify the degree and direction of alignment, and the effects of steric interactions on pattern formation will be presented. Our analysis shows that steric interactions have a destabilizing effect in active suspensions.

Adaptive coatings based on polyaniline for direct 2D observation of diffusion processes in microfluidic systems

[EN] Therefore the understanding and proper evaluation of the flow and mixing behaviour at microscale becomes a very important issue. In this study, the diffusion behaviour of two reacting solutions of HCI and NaOH were directly observed in a glass/polydimethylsiloxane microfluidic device using adaptive coatings based on the conductive polymer polyaniline that are covalently attached to the microchannel walls. The two liquid streams were combined at the junction of a Y-shaped microchannel, and allowed to diffuse into each other and react. The results showed excellent correlation between optical observation of the diffusion process and the numerical results. A numerical model which is based on finite volume method (FVM) discretisation of steady Navier-Stokes (fluid flow) equations and mass transport equations without reactions was used to calculate the flow variables at discrete points in the finite volume mesh element. The high correlation between theory and practical data indicates the potential of such coatings to monitor diffusion processes and mixing behaviour inside microfluidic channels in a dye free environment.

Inverse Spectral Methods in Acoustic Normal Mode Velocimetry of High Reynolds Number Spherical Couette Flows

Experimental geophysical fluid dynamics often examines regimes of fluid flow infeasible for computer simulations. Velocimetry of zonal flows present in these regimes brings many challenges when the fluid is opaque and vigorously rotating; spherical Couette flows with molten metals are one such example. The fine structure of the acoustic spectrum can be related to the fluid’s velocity field, and inverse spectral methods can be used to predict and, with sufficient acoustic data, mathematically reconstruct the velocity field. The methods are to some extent inherited from helioseismology. This work develops a Finite Element Method suitable to matching the geometries of experimental setups, as well as modelling the acoustics based on that geometry and zonal flows therein. As an application, this work uses the 60-cm setup Dynamo 3.5 at the University of Maryland Nonlinear Dynamics Laboratory. Additionally, results obtained using a small acoustic data set from recent experiments in air are provided.