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.

Dynamics of driven three-dimensional thin films: from hdrophilic to superhydrophobic substrates

We study the forced displacement of a thin film of fluid in contact with vertical and inclined substrates of different wetting properties, that range from hydrophilic to hydrophobic, using the lattice-Boltzmann method. We study the stability and pattern formation of the contact line in the hydrophilic and superhydrophobic regimes, which correspond to wedge-shaped and nose-shaped fronts, respectively. We find that contact lines are considerably more stable for hydrophilic substrates and small inclination angles. The qualitative behavior of the front in the linear regime remains independent of the wetting properties of the substrate as a single dispersion relation describes the stability of both wedges and noses. Nonlinear patterns show a clear dependence on wetting properties and substrate inclination angle. The effect is quantified in terms of the pattern growth rate, which vanishes for the sawtooth pattern and is finite for the finger pattern. Sawtooth shaped patterns are observed for hydrophilic substrates and low inclination angles, while finger-shaped patterns arise for hydrophobic substrates and large inclination angles. Finger dynamics show a transient in which neighboring fingers interact, followed by a steady state where each finger grows independently.

Fluid-rock-strain feedbacks in extensional shear zones

La présence de fluide météorique synchrone à l'activité du détachement (Farmin, 2003 ; Mulch et al., 2007 ; Gébelin et al., 2011), implique que les zones de cisaillement sont des systèmes ouverts avec des cellules de convections à l'échelle crustale et un intense gradient géothermique au sein du détachement (Morrison et Anderson, 1998, Gottardi et al., 2011). De plus, les réactions métamorphiques liées à des infiltrations fluides dans les zones de cisaillement extensionnel peuvent influencer les paramètres rhéologiques du système (White and Knipe, 1978), et impliquer la localisation de la déformation dans la croûte. Dans ce manuscrit, deux zones de cisaillement infiltrées par des fluides météoriques sont étudiées, l'une étant largement quartzitique, et l'autre de nature granitique ; les relations entre déformation, fluides, et roches s'appuient sur des approches structurales, microstructurales, chimiques et isotopiques. L'étude du détachement du Columbia river (WA, USA) met en évidence que la déformation mylonitique se développe en un million d'années. La phase de cisaillement principal s'effectue à 365± 30°C d'après les compositions isotopiques en oxygène du quartz et de la muscovite. Ces minéraux atteignent l'équilibre isotopique lors de leur recristallisation dynamique contemporaine à la déformation. La zone de cisaillement enregistre une baisse de température, remplaçant le mécanisme de glissement par dislocation par celui de dissolution- précipitation dans les derniers stades de l'activité du détachement. La dynamique de circulation fluide bascule d'une circulation pervasive à chenalisée, ce qui engendre localement la rupture des équilibres d'échange isotopiques. La zone de cisaillement de Bitterroot (MT, USA) présente une zone mylonitique de 600m d'épaisseur, progressant des protomylonites aux ultramylonites. L'intensité de la localisation de la déformation se reflète directement sur l'hydratation des feldspaths, réaction métamorphique majeure dite de « rock softening ». Une étude sur roche totale indique des transferts de masse latéraux au sein des mylonites, et d'importantes pertes de volume dans les ultramylonites. La composition isotopique en hydrogène des phyllosilicates met en évidence la présence (1) d'une source magmatique/métamorphique originelle, caractérisée par les granodiorites ayant conservé leur foliation magmatique, jusqu'aux protomylonites, et (2) une source météorique qui tamponne les valeurs des phyllosilicates des fabriques mylonitiques jusqu'aux veines de quartz non-déformées. Les compositions isotopiques en oxygène des minéraux illustrent le tamponnement de la composition du fluide météorique par l'encaissant. Ce phénomène cesse lors du processus de chloritisation de la biotite, puisque les valeurs des chlorites sont extrêmement négatives (-10 per mil). La thermométrie isotopique indique une température d'équilibre isotopique de la granodiorite entre 600-500°C, entre 500-300°C dans les mylonites, et entre 300 et 200°C dans les fabriques cassantes (cataclasites et veines de quartz). Basé sur les résultats issus de ce travail, nous proposons un modèle général d'interactions fluide-roches-déformation dans les zones de détachements infiltrées par des fluides météoriques. Les zones de détachements évoluent rapidement (en quelques millions d'années) au travers de la transition fragile-ductile ; celle-ci étant partiellement contrôlée par l'effet thermique des circulations de fluide météoriques. Les systèmes de détachements sont des lieux où la déformation et les circulations fluides sont couplées ; évoluant rapidement vers une localisation de la déformation, et de ce fait, une exhumation efficace. - The presence of meteoric fluids synchronous with the activity of extensional detachment zones (Famin, 2004; Mulch et al., 2007; Gébelin et al., 2011) implies that extensional systems involve fluid convection at a crustal scale, which results in high geothermal gradients within active detachment zones (Morrison and Anderson, 1998, Gottardi et al., 2011). In addition, the metamorphic reactions related to fluid infiltration in extensional shear zones can influence the rheology of the system (White and Knipe, 1978) and ultimately how strain localizes in the crust. In this thesis, two shear zones that were permeated by meteoric fluids are studied, one quartzite-dominated, and the other of granitic composition; the relations between strain, fluid, and evolving rock composition are addressed using structural, microstructural, and chemical/isotopic measurements. The study of the Columbia River detachment that bounds the Kettle core complex (Washington, USA) demonstrates that the mylonitic fabrics in the 100 m thick quartzite- dominated detachment footwall developed within one million years. The main shearing stage occurred at 365 ± 30°C when oxygen isotopes of quartz and muscovite equilibrated owing to coeval deformation and dynamic recrystallization of these minerals. The detachment shear zone records a decrease in temperature, and dislocation creep during detachment shearing gave way to dissolution-precipitation and fracturing in the later stages of detachment activity. Fluid flow switched from pervasive to channelized, leading to isotopic disequilibrium between different minerals. The Bitterroot shear zone detachment (Montana, USA) developed a 600 m thick mylonite zone, with well-developed transitions from protomylonite to ultramylonite. The localization of deformation relates directly to the intensity of feldspar hydration, a major rock- softening metamorphic reaction. Bulk-rock analyses of the mylonitic series indicate lateral mass transfer in the mylonite (no volume change), and significant volume loss in ultramylonite. The hydrogen isotope composition of phyllosilicates shows (1) the presence of an initial magmatic/metamorphic source characterized by the granodiorite in which a magmatic, and gneissic (protomylonite) foliation developed, and (2) a meteoric source that buffers the values of phyllosilicates in mylonite, ultramylonite, cataclasite, and deformed and undeformed quartz veins. The mineral oxygen isotope compositions were buffered by the host-rock compositions until chloritization of biotite started; the chlorite oxygen isotope values are negative (-10 per mil). Isotope thermometry indicates a temperature of isotopic equilibrium of the granodiorite between 600-500°C, between 500-300°C in the mylonite, and between 300 and 200°C for brittle fabrics (cataclasite and quartz veins). Results from this work suggest a general model for fluid-rock-strain feedbacks in detachment systems that are permeated by meteoric fluids. Phyllosilicates have preserved in their hydrogen isotope values evidence for the interaction between rock and meteoric fluids during mylonite development. Fluid flow generates mass transfer along the tectonic anisotropy, and mylonites do not undergo significant volume change, except locally in ultramylonite zones. Hydration of detachment shear zones attends mechanical grain size reduction and enhances strain softening and localization. Self-exhuming detachment shear zones evolve rapidly (a few million years) through the transition from ductile to brittle, which is partly controlled by the thermal effect of circulating surface fluids. Detachment systems are zones in the crust where strain and fluid flow are coupled; these systems. evolve rapidly toward strain localization and therefore efficient exhumation.

Manifolds on the verge of a hyperbolicity breakdown

We study numerically the disappearance of normally hyperbolic invariant tori in quasiperiodic systems and identify a scenario for their breakdown. In this scenario, the breakdown happens because two invariant directions of the transversal dynamics come close to each other, losing their regularity. On the other hand, the Lyapunov multipliers associated with the invariant directions remain more or less constant. We identify notable quantitative regularities in this scenario, namely that the minimum angle between the two invariant directions and the Lyapunov multipliers have power law dependence with the parameters. The exponents of the power laws seem 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).

Amyloid arthropathy in patients undergoing periodical haemodialysis for chronic renal failure: a new complication.

Seven patients (five male and two female) with chronic renal failure (CRF) treated by periodical haemodialysis presented with swelling and effusion of more than three months' duration in knees (four bilateral), shoulders (two, one of them bilateral), elbow (one), and ankle (one). Four had a carpal tunnel syndrome both clinically and electromyographically (three bilateral). All patients had hyperparathyroidism secondary to their CRF, which was not due to amyloidosis in any of them. The dialysis duration period varied from five to 14 years, with an average of 8.6 years. Amyloid deposits (Congo red positive areas with green birefringence under polarising microscopy) were shown in six of the seven synovial biopsy specimens of the knee, in five of the sediments of the synovial fluids, and in specimens removed during carpal tunnel syndrome surgery. No amyloid was found in the biopsy specimen of abdominal fat of six of the patients. The finding of amyloid only in the synovial membrane and fluid, and carpal tunnel, its absence in abdominal fat, and the lack of other manifestations of generalised amyloidosis (cardiomyopathy, malabsorption syndrome, macroglossia, etc.) and of Bence Jones myeloma (protein immunoelectrophoresis normal) raises the possibility that this is a form of amyloidosis which is peculiar to CRF treated by periodical haemodialysis.

Dynamics of Fréedericksz transition in a fluctuating magnetic field

Fréedericksz transition under twist deformation in a nematic layer is discussed when the magnetic field has a random component. A dynamical model which includes the thermal fluctuations of the system is presented. The randomness of the field produces a shift of the instability point. Beyond this instability point the time constant characteristic of the approach to the stationary stable state decreases because of the field fluctuations. The opposite happens for fields smaller than the critical one. The decay time of an unstable state, calculated as a mean first-passage time, is also decreased by the field fluctuations.

"Bottleneck effect" in two-Dimensional microfluidics

An anomalously long transient is needed to achieve a steady pressurization of a fluid when forced to flow through micronarrowed channels under constant mechanical driving. This phenomenon, known as the "bottleneck effect" is here revisited from a different perspective, by using confined displacements of interfacial fluids. Compared to standard microfluidics, such effect admits in this case a neat quantitative characterization, which reveals intrinsic material characteristics of flowing monolayers and permits to envisage strategies for their controlled micromanipulation.

Quasi-two-dimensional electrodeposition under forced fluid flow

Experimental quasi-two-dimensional Zn electrodeposits are grown under forced convection conditions. Large-scale effects, with preferential growth towards the impinging flow, together with small-scale roughness suppression effects are evidenced and separately analyzed by using two different radial cell configurations. Interpretations are given in terms of primary concepts concerning current and concentration distributions.

Diffusion in spatially and temporarily inhomogeneous media: Effects of turbulent mixing

We consider diffusion of a passive substance C in a phase-separating nonmiscible binary alloy under turbulent mixing. The substance is assumed to have different diffusion coefficients in the pure phases A and B, leading to a spatially and temporarily dependent diffusion ¿coefficient¿ in the diffusion equation plus convective term. In this paper we consider especially the effects of a turbulent flow field coupled to both the Cahn-Hilliard type evolution equation of the medium and the diffusion equation (both, therefore, supplemented by a convective term). It is shown that the formerly observed prolonged anomalous diffusion [H. Lehr, F. Sagués, and J.M. Sancho, Phys. Rev. E 54, 5028 (1996)] is no longer seen if a flow of sufficient intensity is supplied.

Stochastic generation of homogeneous isotropic turbulence with well-defined-spectra

A precise and simple computational model to generate well-behaved two-dimensional turbulent flows is presented. The whole approach rests on the use of stochastic differential equations and is general enough to reproduce a variety of energy spectra and spatiotemporal correlation functions. Analytical expressions for both the continuous and the discrete versions, together with simulation algorithms, are derived. Results for two relevant spectra, covering distinct ranges of wave numbers, are given.

Front propagation in spatially modulated media

A two-dimensional reaction-diffusion front which propagates in a modulated medium is studied. The modulation consists of a spatial variation of the local front velocity in the transverse direction to that of the front propagation. We study analytically and numerically the final steady-state velocity and shape of the front, resulting from a nontrivial interplay between the local curvature effects and the global competition process between different maxima of the control parameter. The transient dynamics of the process is also studied numerically and analytically by means of singular perturbation techniques.

Inertial effects on reactive particles advected by turbulence

We study the problem of the advection of passive particles with inertia in a two-dimensional, synthetic, and stationary turbulent flow. The asymptotic analytical result and numerical simulations show the importance of inertial bias in collecting the particles preferentially in certain regions of the flow, depending on their density relative to that of the flow. We also study how these aggregates are affected when a simple chemical reaction mechanism is introduced through a Eulerian scheme. We find that inertia can be responsible for maintaining a stationary concentration pattern even under nonfavorable reactive conditions or destroying it under favorable ones.

Excitability transitions and wave dynamics under spatiotemporal structured noise

We present an analytic and numerical study of the effects of external fluctuations in active media. Our analytical methodology transforms the initial stochastic partial differential equations into an effective set of deterministic reaction-diffusion equations. As a result we are able to explain and make quantitative predictions on the systematic and constructive effects of the noise, for example, target patterns created out of noise and traveling or spiral waves sustained by noise. Our study includes the case of realistic noises with temporal and spatial structures.