Three-dimensional aspects of fluid flows in channels. II. Effects of meniscus and thin film regimes on viscous fingers

We perform a three-dimensional study of steady state viscous fingers that develop in linear channels. By means of a three-dimensional lattice-Boltzmann scheme that mimics the full macroscopic equations of motion of the fluid momentum and order parameter, we study the effect of the thickness of the channel in two cases. First, for total displacement of the fluids in the channel thickness direction, we find that the steady state finger is effectively two-dimensional and that previous two-dimensional results can be recovered by taking into account the effect of a curved meniscus across the channel thickness as a contribution to surface stresses. Second, when a thin film develops in the channel thickness direction, the finger narrows with increasing channel aspect ratio in agreement with experimental results. The effect of the thin film renders the problem three-dimensional and results deviate from the two-dimensional prediction.

Three-dimensional aspects of fluid flows in channels. I. Meniscus and thin film regimes

We study the forced displacement of a fluid-fluid interface in a three-dimensional channel formed by two parallel solid plates. Using a lattice-Boltzmann method, we study situations in which a slip velocity arises from diffusion effects near the contact line. The difference between the slip and channel velocities determines whether the interface advances as a meniscus or a thin film of fluid is left adhered to the plates. We find that this effect is controlled by the capillary and Péclet numbers. We estimate the crossover from a meniscus to a thin film and find good agreement with numerical results. The penetration regime is examined in the steady state. We find that the occupation fraction of the advancing finger relative to the channel thickness is controlled by the capillary number and the viscosity contrast between the fluids. For high viscosity contrast, lattice-Boltzmann results agree with previous results. For zero viscosity contrast, we observe remarkably narrow fingers. The shape of the finger is found to be universal.

Origin of siderite mineralisation in Petrova and Trgovska Gora Mts., NW Dinarides

The Petrova and Trgovska Gora Mts. (Gora=Mountain) are Variscan basement units incorporated into the northwestern Dinarides during the Alpine orogeny. They host numerous siderite-quartz-polysulphide, siderite-chalcopyrite, siderite-galena and barite veins, as well as stratabound hydrothermal-replacement ankerite bodies within carbonates in non-metamorphosed, flysch-like Permo-Carboniferous sequences. The deposits have been mined for Cu, Pb, Ag and Fe ores since Medieval times. Fluid inclusion studies of quartz from siderite-polysulphide-quartz and barite veins of both regions have shown the presence of primary aqueous NaCl-CaCl(2)+/- MgCl(2)-H(2)O +/- CO(2) inclusions. The quartz-sulphide stage of both regions show variable salinities; 2.7-26.2 wt% NaCl eq. for the Trgovska Gora region and 3.4-23.4 wt% NaCl eq. for the Petrova gora region, and similar homogenisation temperatures (100-230A degrees C). Finally, barite is precipitated from low salinity-low temperature solutions (3.7-15.8 wt % NaCl equ. and 115-145A degrees C). P-t conditions estimated via isochore construction yield formation temperatures between 180-250A degrees C for the quartz-sulphide stage and 160-180A degrees C for the barite stage, using a maximum lithostatic pressure of 1 kbar (cc. 3 km of overburden). The sulphur isotope composition of barite from both deposits indicates the involvement of Permian seawater in ore fluids. This is supported by the elevated bromium content of the fluid inclusion leachates (120-660 ppm in quartz, 420-960 ppm in barite) with respect to the seawater, indicating evaporated seawater as the major portion of the ore-forming fluids. Variable sulphur isotope compositions of galena, pyrite and chalcopyrite, between -3.2 and +2.7aEuro degrees, are interpreted as a product of incomplete thermal reduction of the Permian marine sulphate mixed with organically- and pyrite-bound sulphur from the host sedimentary rocks. Ore-forming fluids are interpreted as deep-circulating fluids derived primarily from evaporated Permian seawater and later modified by interaction with the Variscan basement rocks. (40)Ar/(39)Ar data of the detrital mica from the host rocks yielded the Variscan age overprinted by an Early Permian tectonothermal event dated at 266-274 Ma. These ages are interpreted as those reflecting hydrothermal activity correlated with an incipient intracontinental rifting in the Tethyan domain. Nevertheless, 75 Ma recorded at a fine-grained sericite sample from the alteration zone is interpreted as a result of later resetting of white mica during Campanian opening/closure of the Sava back arc in the neighbouring Sava suture zone (Ustaszewski et al. 2008).

Origin of CO2 and carbonate veins in mantle-derived xenoliths in the Pannonian Basin

The origin and evolution of CO2 inclusions and calcite veins in peridotite xenoliths of the Pannonian Basin, Hungary, were investigated by means of petrographic investigation and stable isotope analyses. The fluid inclusions recovered in paragenetic olivine and clinopyroxene belong to distinct populations: type A (texturally early) inclusions with regular shapes (often with negative crystal forms) forming intragranular trails, type B (texturally late) inclusions defining randomly oriented trails that reach grain boundaries Type B inclusions are often associated with silicate melt (type C) inclusions Stable carbon isotope compositions in inclusion-hosted CO2 were obtained by vacuum crushing followed by conventional dual inlet as well as continuous flow mass spectrometry in order to eliminate possible lab artifacts. Olivines, clino- and orthopyroxenes of the host peridotite have oxygen isotope compositions from 5.3 to 6.0 parts per thousand (relative to V-SMOW). without any relationship with xenolith texture. Some of the xenoliths contained calcite in various forms veins and infillings in silicate globules in veins, secondary carbonate veins filling cracks and metasomatic veins with diffuse margins The former two carbonate types have delta C-13 values around -13 parts per thousand (relative to V-PDB) and low Sr contents (<05 wt %), whereas the third type,veins with high-temperature metasomatic features have a delta C-13 value of -5 0 parts per thousand and high Sr contents up to 34 wt.% In spite of the mantle-like delta C-13 value and the unusually high Sr content typical for mantle-derived carbonate, trace element compositions have proven a crustal origin. This observation supports the conclusions of earlier studies that the carbonate melt droplets found on peridotite xenoliths in the alkaline basalts represent mobilized sedimentary carbonate. The large delta C-13 range and the C-12-enrichment in the carbonates can be attributed to devolanlization of the migrating carbonate or infiltration of surficial fluids containing C-12-rich dissolved carbon Carbon isotope compositions of inclusion-hosted CO2 range from -17 8 to -4.8 parts per thousand (relative to V-PDB) with no relation to the amount of CO2 released by vacuum crushing. Low-delta C-13 values measured by stepwise heating under vacuum suggest that the carbon component is pristine and not related to surficial contamination, and that primary mantle fluids with delta C-13 values around -5 parts per thousand were at least partly preserved in the xenoliths Tectonic reworking and heating by the basaltic magma resulted in partial CO2 release and local C-13-depletion. (C) 2010 Elsevier B V All rights reserved

Magmatic fluids in the Breccia-hosted epithermal Au-Ag deposit of Rosia Montana, Romania

The breccia-hosted epithermal Au-Ag deposit of Rosia Montana is located 7 kin northeast of Abrud, in the northern part of the South Apuseni Mountains, Romania. Estimated total reserves of 214.91 million metric toils (Mt) of ore at 1.46 g/t An and 6.9 g/t Ag (10.1 Moz of An and 47.6 Moz of Ag) make Rosia Montana one of the largest gold deposits in Europe. At this location, Miocene calc-alkaline magmatic and hydrothermal activity was associated with local extensional tectonics within a strike-slip regime related to the indentation of the Adriatic microplate into the European plate during the Carpathian orogenesis. The host rocks of the magmatic complex consist of pre-Mesozoic metamorphosed continental crust covered by Cretaceous turbiditic sediment (flysch). Magmatic activity at Rosia Montana and its surroundings occurred in several pulses and lasted about 7 m.y, Rosia Montana is a breccia-hosted epithermal system related to strong phreatomagmatic activity due to the shallow emplacement of the Montana dacite. The Montana dacite intruded Miocene volcaniclastic material (volcaniclastic breccias) and crops out at Cetate and Carnic Hills. Current mining is focused primarily on the Cetate open pit, which was mapped in detail, leading to the recognition of three distinct breccia bodies: the dacite breccia with a dominantly hydrothermal matrix, the gray polymict breccia with a greater proportion of sand-sized matrix support, and the black polymict breccia, which reached to the surface, contains carbonized tree trunks and has a dominantly barren elastic matrix. The hydrothermal alteration is pervasive. Adularia alteration with a phyllic overprint is ubiquitous; silicification and argillic alteration occur locally. Mineralization consists of quartz, adularia, carbonates (commonly Mn-rich), pyrite, Fe-poor sphalerite, galena, chalcopyrite, tetrahedrite, and native gold and occurs as disseminations, as well as in veins and filling vugs within the Montana dacite and the different breccias. The age of mineralization (12.85 +/- 0.07 Ma) was determined by Ar-40- Ar-39 dating on hydrothermal adularia crystals from vugs in the dacite breccia in the Cetate open pit. Microthermometric measurements of fluid inclusions in quartz phenocrysts from the Montana dacite revealed two fluid types that are absent from the hydrothermal breccia and must have been trapped at depth prior to dacite dome emplacement: brine inclusions (32-55 -wt % NaCl equiv, homogenizing at T-h > 460 degrees C) and intermediate density fluids (4.9-15.6 wt % NaCl equiv, T, between 345 degrees-430 degrees C). Secondary aqueous fluid inclusion assemblages in the phenocrysts have salinities of 0.2 to 2.2 wt percent NaCl equiv and T-h of 200 degrees to 280 degrees C. Fluid inclusion assemblages in hydrothermal quartz from breccias and veins have salinities of 0.2 to 3.4 wt percent NaCl equiv and T-h, from 200 degrees to 270 degrees C. The oxygen isotope composition of several zones of an ore-related epithermal quartz crystal indicate a very constant delta O-18 of 4.5 to 5.0 per mil for the mineralizing fluid, despite significant salinity and temperature variation over time. Following microthermometry, selected fluid inclusion assemblages were analyzed by laser ablation-inductively coupled-plasma mass spectrometry (LA-ICMS). Despite systematic differences in salinity between phenocryst-hosted fluids trapped at depth and fluids from quartz in the epithermal breccias, all fluids have overlapping major and trace cation ratios, including identical Na/K/Rb/Sr/Cs/Ba. Consistent with the constant near-magmatic oxygen isotope composition of the hydrothermal fluids, these data strongly indicate a common magmatic component of these chemically conservative solutes in all fluids. Cu, Pb, Zn, and Mn show variations in concentration relative to the relatively non-reactive alkalis, reflecting the precipitation of sulfide minerals together with An in the epithermal breccia, and possibly of Cu in an inferred subjacent porphyry environment. The magmatic-hydrothermal processes responsible for epithermal Au-Ag mineralization at Rosia Montana are, however, not directly related to the formation of the spatially associated porphyry Cu-Au deposit of Rosia Poieni, which occurred lout 3 m.y. later.

Mineral zoning and geochemistry of epithermal polymetallic Zn-Pb-Ag-Cu-Bi mineralization at Cerro de Pasco, Peru

The large Cerro de Pasco Cordilleran base metal deposit in central Peru is located on the eastern margin of a middle Miocene diatreme-dome complex and comprises two mineralization stages. The first stage consists of a large pyrite-quartz body replacing Lower Mesozoic Pucara carbonate rocks and, to a lesser extent, diatreme breccia. This body is composed of pyrite with pyrrhotite inclusions, quartz, and black and red chalcedony (containing hypogene hematite). At the contact with the pyrite-quartz body, the diatreme breccia is altered to pyrite-quartz-sericite-pyrite. This body was, in part, replaced by pipelike pyrrhotite bodies zoned outward to carbonate-replacement Zn-Pb ores hearing Fe-rich sphalerite (up to 24 mol % Fes). The second mineralization stage is partly superimposed on the first and consists of zoned east-west-trending Cu-Ag-(Au-Zn-Pb) enargite-pyrite veins hosted in the diatreme breccia in the western part of the deposit and well-zoned Zn-Pb-(Bi-Ag-Cu) carbonate-replacement orebodies; in both cases, sphalerite is Fe poor and the inner parts of the orebodies show typically advanced argillic alteration assemblages, including aluminum phosphate Sulfate (APS) minerals. The zoned enargite-pyrite veins display mineral zoning, from a core of enargite-pyrite +/- alunite with traces of Au, through an intermediate zone of tennantite, chalcopyrite, and Bi minerals to a poorly developed Outer zone hearing sphalerite-galena +/- kaolinite. The carbonate-hosted replacement ores are controlled along N 35 degrees E, N 90 degrees E, N 120 degrees E, and N 170 degrees E faults. They form well-zoned upward-flaring pipelike orebodies with a core of famatinite-pyrite and alunite, an intermediate zone with tetrahedrite-pyrite, chalcopyrite, matildite, cuprobismutite, emplectite, and other Bi minerals accompanied by APS minerals, kaolinite, and dickite, and an outer zone composed of Fe-poor sphalerite (in the range of 0.05-3.5 mol % Fes) and galena. The outermost zone consists of hematite, magnetite, and Fe-Mn-Zn-Ca-Mg carbonates. Most of the second-stage carbonate-replacement orebodies plunge between 25 degrees and 60 degrees to the west, suggesting that the hydrothermal fluids ascended from deeper levels and that no lateral feeding from the veins to the carbonate-replacement orebodies took place. In the Venencocha and Santa Rosa areas, located 2.5 km northwest of the Cerro de Pasco open pit and in the southern part of the deposit, respectively, advanced argillic altered dacitic domes and oxidized veins with advanced argillic alteration halos occur. The latter veins are possibly the oxidized equivalent of the second-stage enargite-pyrite veins located in the western part of the deposit. The alteration assemblage quartz-muscovite-pyrite associated with the pyrite-quartz body suggests that the first stage precipitated at slightly, acidic fin. The sulfide mineral assemblages define an evolutionary path close to the pyrite-pyrrhotite boundary and are characteristic of low-sulfidation states; they suggest that the oxidizing slightly acidic hydrothermal fluid was buffered by phyllite, shale, and carbonate host rock. However, the presence in the pyrite-quartz body of hematite within quartz suggests that, locally, the fluids were less buffered by the host rock. The mineral assemblages of the second mineralization stage are characteristic of high- to intermediate-sulfidation states. High-sulfidation states and oxidizing conditions were achieved and maintained in the cores of the second-stage orebodies, even in those replacing carbonate rocks. The observation that, in places, second-stage mineral assemblages are found in the inner and outer zones is explained in terms of the hydrothermal fluid advancing and waning. Microthermometric data from fluid inclusions in quartz indicate that the different ores of the first mineralization stage formed at similar temperatures and moderate salinities (200 degrees-275 degrees C and 0.2-6.8 wt % NaCl equiv in the pyrite-quartz body; 192 degrees-250 degrees C and 1.1-4.3 wt % NaCl equiv in the pyrrhotite bodies; and 183 degrees-212 degrees C and 3.2-4.0 wt % NaCl equiv in the Zn-Pb ores). These values are similar to those obtained for fluid inclusions in quartz and sphalerite from the second-stage ores (187 degrees-293 degrees C and 0.2-5.2 wt % NaCl equiv in the enargite-pyrite veins: 178 degrees-265 degrees C and 0.2-7.5 wt % NaCl equiv in quartz of carbonate-replacement orebodies; 168 degrees-999 degrees C and 3-11.8 wt % NaCl equiv in sphalerite of carbonate-replacement orebodies; and 245 degrees-261 degrees C and 3.2-7.7 wt % NaCl equiv in quartz from Venencocha). Oxygen and hydrogen isotope compositions oil kaolinite from carbonate-replacement orebodies (delta(18)O = 5.3-11.5%o, delta D = -82 to -114%o) and on alunite from the Venencocha and Santa Rosa areas (delta(18)O = 1.9-6.9%o, delta D = -56 to -73%o). Oxygen isotope compositions of quartz from the first and second stages have 6180 values from 9.1 to 1.7.8 per mil. Calculated fluids in equilibrium with kaolinite have delta(18)O values of 2.0 to 8.2 and delta D values of -69 to -97 per mil; values in equilibrium with alunite are -1.4 to -6.4 and -62 to -79 per mil. Sulfur isotope compositions of sulfides from both stages have a narrow range of delta(34)S values, between -3.7 and +4.2 per mil; values for sulfates from the second stage are between 4.2 and 31.2 per mil. These results define two mixing trends for the ore-forming fluids. The first trend reflects mixing between a moderately saline (similar to 10 wt % NaCl equiv) magmatic end member that had degassed (as indicated by the low delta D values) and meteoric water. The second mixing indicates condensation of magmatic vapor with HCl and SO(2) into meteoric water, which formed alunite. The hydrothermal system at Cerro de Pasco was emplaced at a shallow depth (similar to 500 m) in the epithermal and upper part of a porphyry environment. The similar temperatures and salinities obtained for the first stage and second stages, together with the stable isotope data, indicate that both stages are linked and represent successive stages of epithermal polymetallic mineralization in the upper part of a porphyry system.

Stable isotope (C, O, S) systematics of the mercury mineralization at Idrija, Slovenia: constraints on fluid source and alteration processes

The world-class Idrija mercury deposit (western Slovenia) is hosted by highly deformed Permocarboniferous to Middle Triassic sedimentary rocks within a complex tectonic structure at the transition between the External Dinarides and the Southern Alps. Concordant and discordant mineralization formed concomitant with Middle Triassic bimodal volcanism in an aborted rift. A multiple isotopic (C, O, S) investigation of host rocks and ore minerals was performed to put constraints on the source and composition of the fluid, and the hydrothermal alteration. The distributions of the delta(13)C and delta(18)O values of host and gangue carbonates are indicative of a fracture-controlled hydrothermal system, with locally high fluid-rock ratios. Quantitative modeling of the delta(13)C and delta(18)O covariation for host carbonates during temperature dependent fluid-rock interaction, and concomitant precipitation of void-filling dolomites points to a slightly acidic hydrothermal fluid (delta(13)Capproximate to-4parts per thousand and delta(18)Oapproximate to+10parts per thousand), which most likely evolved during isotopic exchange with carbonates under low fluid/rock ratios. The delta(34)S values of hydrothermal and sedimentary sulfur minerals were used to re-evaluate the previously proposed magmatic and evaporitic sulfur sources for the mineralization, and to assess the importance of other possible sulfur sources such as the contemporaneous seawater sulfate, sedimentary pyrite, and organic sulfur compounds. The delta(34)S values of the sulfides show a large variation at deposit down to hand-specimen scale. They range for cinnabar and pyrite from -19.1 to +22.8parts per thousand, and from -22.4 to +59.6parts per thousand, respectively, suggesting mixing of sulfur from different sources. The peak of delta(34)S values of cinnabar and pyrite close to 0parts per thousand is compatible with ore sulfur derived dominantly from a magmatic fluid and/or from hydrothermal leaching of basement rocks. The similar stratigraphic trends of the delta(34)S values of both cinnabar and pyrite suggest a minor contribution of sedimentary sulfur (pyrite and organic sulfur) to the ore formation. Some of the positive delta(34)S values are probably derived from thermochemical reduction of evaporitic and contemporaneous seawater sulfates.

The Imprint of Geologic History on Within‐Island Diversification of Woodlouse-­Hunter Spiders (Araneae, Dysderidae) in the Canary Islands.

Geological processes and ecological adaptation are major drivers of diversification on oceanic islands. Although diversification in these islands is often interpreted as resulting from dispersal or island hopping rather than vicariance, this may not be the case in islands with complex geological histories. The island of Tenerife, in the Canary Islands, emerged in the late Miocene as 3 precursor islands that were subsequently connected and reisolated by volcanic cycles. The spider Dysdera verneaui is endemic to the island of Tenerife, where it is widely distributed throughout most island habitats, providing an excellent model to investigate the role of physical barriers and ecological adaptation in shaping within-island diversity. Here, we present evidence that the phylogeographic patterns of this species trace back to the independent emergence of the protoislands. Molecular markers (mitochondrial genes cox1, 16S, and nad1 and the nuclear genes ITS-2 and 28S) analyzed from 100 specimens (including a thorough sampling of D. verneaui populations and additional outgroups) identify 2 distinct evolutionary lineages that correspond to 2 precursor islands, each with diagnostic genital characters indicative of separate species status. Episodic introgression events between these 2 main evolutionary lineages explain the observed incongruence between mitochondrial and nuclear markers, probably as a result of the homogenization of their ITS-2 sequence types. The most widespread lineage exhibits a complex population structure, which is compatible with either secondary contact, following connection of deeply divergent lineages, or alternatively, a back colonization from 1 precursor island to another.

The Imprint of Geologic History on Within‐Island Diversification of Woodlouse-­Hunter Spiders (Araneae, Dysderidae) in the Canary Islands.

Geological processes and ecological adaptation are major drivers of diversification on oceanic islands. Although diversification in these islands is often interpreted as resulting from dispersal or island hopping rather than vicariance, this may not be the case in islands with complex geological histories. The island of Tenerife, in the Canary Islands, emerged in the late Miocene as 3 precursor islands that were subsequently connected and reisolated by volcanic cycles. The spider Dysdera verneaui is endemic to the island of Tenerife, where it is widely distributed throughout most island habitats, providing an excellent model to investigate the role of physical barriers and ecological adaptation in shaping within-island diversity. Here, we present evidence that the phylogeographic patterns of this species trace back to the independent emergence of the protoislands. Molecular markers (mitochondrial genes cox1, 16S, and nad1 and the nuclear genes ITS-2 and 28S) analyzed from 100 specimens (including a thorough sampling of D. verneaui populations and additional outgroups) identify 2 distinct evolutionary lineages that correspond to 2 precursor islands, each with diagnostic genital characters indicative of separate species status. Episodic introgression events between these 2 main evolutionary lineages explain the observed incongruence between mitochondrial and nuclear markers, probably as a result of the homogenization of their ITS-2 sequence types. The most widespread lineage exhibits a complex population structure, which is compatible with either secondary contact, following connection of deeply divergent lineages, or alternatively, a back colonization from 1 precursor island to another.

The Imprint of Geologic History on Within‐Island Diversification of Woodlouse-­Hunter Spiders (Araneae, Dysderidae) in the Canary Islands.

Geological processes and ecological adaptation are major drivers of diversification on oceanic islands. Although diversification in these islands is often interpreted as resulting from dispersal or island hopping rather than vicariance, this may not be the case in islands with complex geological histories. The island of Tenerife, in the Canary Islands, emerged in the late Miocene as 3 precursor islands that were subsequently connected and reisolated by volcanic cycles. The spider Dysdera verneaui is endemic to the island of Tenerife, where it is widely distributed throughout most island habitats, providing an excellent model to investigate the role of physical barriers and ecological adaptation in shaping within-island diversity. Here, we present evidence that the phylogeographic patterns of this species trace back to the independent emergence of the protoislands. Molecular markers (mitochondrial genes cox1, 16S, and nad1 and the nuclear genes ITS-2 and 28S) analyzed from 100 specimens (including a thorough sampling of D. verneaui populations and additional outgroups) identify 2 distinct evolutionary lineages that correspond to 2 precursor islands, each with diagnostic genital characters indicative of separate species status. Episodic introgression events between these 2 main evolutionary lineages explain the observed incongruence between mitochondrial and nuclear markers, probably as a result of the homogenization of their ITS-2 sequence types. The most widespread lineage exhibits a complex population structure, which is compatible with either secondary contact, following connection of deeply divergent lineages, or alternatively, a back colonization from 1 precursor island to another.

Fluid flow processes at basin scale

S u b s u r face fluid flow plays a significant role in many geologic processes and is increasingly being studied in the scale of sedimentary basins and geologic time perspective. Many economic resources such as petroleum and mineral deposits are products of basin scale fluid flow operating over large periods of time. Such ancient flow systems can be studied through analysis of diagenetic alterations and fluid inclusions to constrain physical and chemical conditions of fluids and rocks during their paleohy d r og e o l ogic evolution. Basin simulation models are useful to complement the paleohy d r og e o l ogic record preserved in the rocks and to derive conceptual models on hydraulic basin evolution and generation of economic resources. Different types of fluid flow regimes may evo l ve during basin evolution. The most important with respect to flow rates and capacity for transport of solutes and thermal energy is gr avitational fluid flow driven by the topographic configuration of a basin. Such flow systems require the basin to be elevated above sea level. Consolidational fluid flow is the principal fluid migration process in basins below sea level, caused by loading of compressible rocks. Flow rates of such systems are several orders of magnitude below topogr a p hy driven flow. Howeve r, consolidation may create significant fluid ove rpressure. Episodic dewatering of ove rpressured compart m e n t s m ay cause sudden fluid release with elevated flow velocities and may cause a transient local thermal and chemical disequilibrium betwe e n fluid and rock. This paper gives an ove rv i ew on subsurface fluid flow processes at basin scale and presents examples related to the Pe n e d è s basin in the central Catalan continental margin including the offshore Barcelona half-graben and the compressive South-Pyrenean basin.

The influence of a certain volume of fluid in the lower tear meniscus on intraocular pressure measurements using the non-contact tonometer

Background: Intraocular pressure (IOP) is the pressure inside the eye that helps to maintain the integrity and the suitable form of the ocular globe. Precise and accurate measures of IOP are needed for the diagnosis as well as follow-up of glaucoma. In daily clinical practice, Goldmann applanation tonometer (GAT) and Non-contact tonometer (NCT) are the most common devices for measuring IOP. A close agreement between these methods has been showed, particularly in normotensive patients and a poor agreement, especially when IOP levels are above the normal range. Ophthalmologists have noticed a poor agreement between NCT and GAT, observing that by using NCT and after comparing with GAT, there is an overestimation of IOP readings, and particularly it occurs when the eyes are tearful. Previous studies investigate the effect of tears in Non-contact tonometer readings by the instillation of artificial tears, concluding in one of the studies that the variation was less than 1mmHg and not clinically significant, in contrast with another study which the increases were sadistically significant. Tear menisci are a thin strip of tear fluid located between the bulbar conjunctiva and the eyelid margins. We think that the overestimation of IOP readings using NCT could be due to the presence of a higher volume of tear in the lower tear meniscus which might cause an optical interference in the optoelectronic applanation monitoring system of this deviceObjectives: To research the influence of a certain volume of fluid in the lower tear meniscus on IOP measurements using the NCT in healthy eyes. Moreover, to investigate the agreement between IOP readings obtained by NCT and GAT in the presence and absence of this volume of fluidMethods: The study design will be transversal for diagnostic tests of repeated measures. We will study patients with no ocular pathology and IOP<21mmHg. It will consist in the measurement of IOP using NCT before and after the instillation of COLIRCUSÍ FLUOTEST, used as a volume of fluid in the lower tear meniscus, to observe if there will be differences using the paired t-test. Moreover, we will take IOP measures by GAT in order to know the agreement between these methods after and before the application of these eyedrops, using the ICC (intraclass correlation coefficient) and the Bland-Altmann method

HPLC-ESI-MS/MS analysis of oxidized di-caffeoylquinic acids generated by metalloporphyrin-catalyzed reactions

This paper reports an HPLC-ESI-MS/MS investigation on the oxidation of 3,5- and 4,5- dicaffeoylquinic acid using iron(III) tetraphenylporphyrin chloride as catalyst. Two major mono-oxidised products of the quinic acid moiety have been identified for both compounds. However, only the 4,5-derivative afforded two different tri-oxo products. Thus, it seems that the oxidation pattern depends on the number and positions of the caffeic acid moieties present in caffeoylquinic acid molecules.

Methods and Models for Accelerating Dynamic Simulation of Fluid Power Circuits

The objective of this dissertation is to improve the dynamic simulation of fluid power circuits. A fluid power circuit is a typical way to implement power transmission in mobile working machines, e.g. cranes, excavators etc. Dynamic simulation is an essential tool in developing controllability and energy-efficient solutions for mobile machines. Efficient dynamic simulation is the basic requirement for the real-time simulation. In the real-time simulation of fluid power circuits there exist numerical problems due to the software and methods used for modelling and integration. A simulation model of a fluid power circuit is typically created using differential and algebraic equations. Efficient numerical methods are required since differential equations must be solved in real time. Unfortunately, simulation software packages offer only a limited selection of numerical solvers. Numerical problems cause noise to the results, which in many cases leads the simulation run to fail. Mathematically the fluid power circuit models are stiff systems of ordinary differential equations. Numerical solution of the stiff systems can be improved by two alternative approaches. The first is to develop numerical solvers suitable for solving stiff systems. The second is to decrease the model stiffness itself by introducing models and algorithms that either decrease the highest eigenvalues or neglect them by introducing steady-state solutions of the stiff parts of the models. The thesis proposes novel methods using the latter approach. The study aims to develop practical methods usable in dynamic simulation of fluid power circuits using explicit fixed-step integration algorithms. In this thesis, twomechanisms whichmake the systemstiff are studied. These are the pressure drop approaching zero in the turbulent orifice model and the volume approaching zero in the equation of pressure build-up. These are the critical areas to which alternative methods for modelling and numerical simulation are proposed. Generally, in hydraulic power transmission systems the orifice flow is clearly in the turbulent area. The flow becomes laminar as the pressure drop over the orifice approaches zero only in rare situations. These are e.g. when a valve is closed, or an actuator is driven against an end stopper, or external force makes actuator to switch its direction during operation. This means that in terms of accuracy, the description of laminar flow is not necessary. But, unfortunately, when a purely turbulent description of the orifice is used, numerical problems occur when the pressure drop comes close to zero since the first derivative of flow with respect to the pressure drop approaches infinity when the pressure drop approaches zero. Furthermore, the second derivative becomes discontinuous, which causes numerical noise and an infinitely small integration step when a variable step integrator is used. A numerically efficient model for the orifice flow is proposed using a cubic spline function to describe the flow in the laminar and transition areas. Parameters for the cubic spline function are selected such that its first derivative is equal to the first derivative of the pure turbulent orifice flow model in the boundary condition. In the dynamic simulation of fluid power circuits, a tradeoff exists between accuracy and calculation speed. This investigation is made for the two-regime flow orifice model. Especially inside of many types of valves, as well as between them, there exist very small volumes. The integration of pressures in small fluid volumes causes numerical problems in fluid power circuit simulation. Particularly in realtime simulation, these numerical problems are a great weakness. The system stiffness approaches infinity as the fluid volume approaches zero. If fixed step explicit algorithms for solving ordinary differential equations (ODE) are used, the system stability would easily be lost when integrating pressures in small volumes. To solve the problem caused by small fluid volumes, a pseudo-dynamic solver is proposed. Instead of integration of the pressure in a small volume, the pressure is solved as a steady-state pressure created in a separate cascade loop by numerical integration. The hydraulic capacitance V/Be of the parts of the circuit whose pressures are solved by the pseudo-dynamic method should be orders of magnitude smaller than that of those partswhose pressures are integrated. The key advantage of this novel method is that the numerical problems caused by the small volumes are completely avoided. Also, the method is freely applicable regardless of the integration routine applied. The superiority of both above-mentioned methods is that they are suited for use together with the semi-empirical modelling method which necessarily does not require any geometrical data of the valves and actuators to be modelled. In this modelling method, most of the needed component information can be taken from the manufacturer’s nominal graphs. This thesis introduces the methods and shows several numerical examples to demonstrate how the proposed methods improve the dynamic simulation of various hydraulic circuits.

Fluid dynamic modelling and simulation of circulating fluidized bed reactors: importance of the interface turbulence transfer

The main objective of this work is to analyze the importance of the gas-solid interface transfer of the kinetic energy of the turbulent motion on the accuracy of prediction of the fluid dynamic of Circulating Fluidized Bed (CFB) reactors. CFB reactors are used in a variety of industrial applications related to combustion, incineration and catalytic cracking. In this work a two-dimensional fluid dynamic model for gas-particle flow has been used to compute the porosity, the pressure, and the velocity fields of both phases in 2-D axisymmetrical cylindrical co-ordinates. The fluid dynamic model is based on the two fluid model approach in which both phases are considered to be continuous and fully interpenetrating. CFB processes are essentially turbulent. The model of effective stress on each phase is that of a Newtonian fluid, where the effective gas viscosity was calculated from the standard k-epsilon turbulence model and the transport coefficients of the particulate phase were calculated from the kinetic theory of granular flow (KTGF). This work shows that the turbulence transfer between the phases is very important for a better representation of the fluid dynamics of CFB reactors, especially for systems with internal recirculation and high gradients of particle concentration. Two systems with different characteristics were analyzed. The results were compared with experimental data available in the literature. The results were obtained by using a computer code developed by the authors. The finite volume method with collocated grid, the hybrid interpolation scheme, the false time step strategy and SIMPLEC (Semi-Implicit Method for Pressure Linked Equations - Consistent) algorithm were used to obtain the numerical solution.