Estimating traveltimes and groundwater flow patterns using 3D time-lapse crosshole ERT imaging of electrical resistivity fluctuations induced by infiltrating river water

The infiltration of river water into aquifers is of high relevance to drinking-water production and is a key driver of biogeochemical processes in the hyporheic and riparian zone, but the distribution and quantification of the infiltrating water are difficult to determine using conventional hydrological methods (e.g., borehole logging and tracer tests). By time-lapse inverting crosshole ERT (electrical resistivity tomography) monitoring data, we imaged groundwater flow patterns driven by river water infiltrating a perialpine gravel aquifer in northeastern Switzerland. This was possible because the electrical resistivity of the infiltrating water changed during rainfall-runoff events. Our time-lapse resistivity models indicated rather complex flow patterns as a result of spatially heterogeneous bank filtration and aquifer heterogeneity. The upper part of the aquifer was most affected by the river infiltrate, and the highest groundwater velocities and possible preferential flow occurred at shallow to intermediate depths. Time series of the reconstructed resistivity models matched groundwater electrical resistivity data recorded on borehole loggers in the upper and middle parts of the aquifer, whereas the resistivity models displayed smaller variations and delayed responses with respect to the logging data. in the lower part. This study demonstrated that crosshole ERT monitoring of natural electrical resistivity variations of river infiltrate could be used to image and quantify 3D bank filtration and aquifer dynamics at a high spatial resolution.

Fluid dynamical prediction of changed v1-flow at energies available at the CERN Large Hadron Collider

Substantial collective flow is observed in collisions between lead nuclei at Large Hadron Collider (LHC) as evidenced by the azimuthal correlations in the transverse momentum distributions of the produced particles. Our calculations indicate that the global v1-flow, which at RHIC peaked at negative rapidities (named third flow component or antiflow), now at LHC is going to turn toward forward rapidities (to the same side and direction as the projectile residue). Potentially this can provide a sensitive barometer to estimate the pressure and transport properties of the quark-gluon plasma. Our calculations also take into account the initial state center-of-mass rapidity fluctuations, and demonstrate that these are crucial for v1 simulations. In order to better study the transverse momentum flow dependence we suggest a new"symmetrized" vS1(pt) function, and we also propose a new method to disentangle global v1 flow from the contribution generated by the random fluctuations in the initial state. This will enhance the possibilities of studying the collective Global v1 flow both at the STAR Beam Energy Scan program and at LHC.

Minor salivary gland tumors: a clinicopathological study of 18 cases

Introduction: Minor salivary gland tumors (MSGTs) are infrequent, representing 10-15% of all salivary neoplasms. Despite this low frequency, MSGTs conform a heterogeneous group of neoplasms characterized by a broad range of histological types. Patients and method: We identified cases of MSGT in a retrospective study of the biopsies made in the period 1997-2007 in the Service of Oral Surgery (Dental Clinic of the University of Barcelona, Spain). The data collected comprised patient age and sex, the clinical characteristics and location of the tumor, the duration of the lesion, its size, the treatment provided, and the histopathological findings. Results: Of the 18 cases of MSGT studied, 12 corresponded to women (66.7%) and 6 to men (33.3%). The great majority (94.4%) were benign tumors. The preferential location was the posterior third of the hard palate (33.2%), followed by the soft palate (16.7%) and the mucosa of the upper lip (16.7%). The histopathological diagnoses of our MSGTs comprised 10 pleomorphic adenomas (55.3%), 2 cystadenomas (11.1%), 1 myoepithelioma (5.6%), 1 sialadenoma papilliferum (5.6%), 1 basal cell adenoma (5.6%), 1 Warthin"s tumor (5.6%), 1 canalicular adenoma (5.6%), and 1 low-grade polymorphic adenocarcinoma (5.6%). Discussion and conclusions: Coinciding with our own results, the literature describes a high recurrence rate for MSGTs (5-30%) when surgical removal is incomplete. Six percent of all benign minor salivary gland tumors are considered to relapse, versus 65% of all malignant lesions. Periodic clinical controls are required, since the possibility of malignant transformation must be taken into account

Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

Diffusion in stationary flow from mesoscopic nonequilibrium thermodynamics

We analyze the diffusion of a Brownian particle in a fluid under stationary flow. By using the scheme of nonequilibrium thermodynamics in phase space, we obtain the Fokker-Planck equation that is compared with others derived from the kinetic theory and projector operator techniques. This equation exhibits violation of the fluctuation-dissipation theorem. By implementing the hydrodynamic regime described by the first moments of the nonequilibrium distribution, we find relaxation equations for the diffusion current and pressure tensor, allowing us to arrive at a complete description of the system in the inertial and diffusion regimes. The simplicity and generality of the method we propose makes it applicable to more complex situations, often encountered in problems of soft-condensed matter, in which not only one but more degrees of freedom are coupled to a nonequilibrium bath.

Impaired chemical coupling of cerebral blood flow is compatible with intact neurological outcome in neonates with perinatal risk factors.

Early detection of pathophysiological factors associated with permanent brain damage is a major issue in neonatal medicine. The aim of our study was to evaluate the significance of the CO2 reactivity of cerebral blood flow (CBF) in neonates with perinatal risk factors. Fourteen ventilated neonates with perinatal risk factors (pathological cardiotocogramm, low cord pH, postpartal encephalopathy) were enrolled into this prospective study. The study was performed 18-123 h after birth. CBF was measured using the noninvasive intravenous 133Xe method. Two measurements were taken with a minimal PaCO2-difference of 5 mm Hg. From the two CBF values the CO2 reactivity was calculated. Outcome was evaluated 1 year after birth. The CBF values at a lower PaCO2 ranged from 6.6 to 115. 2 ml/100 g brain issue/min (median = 18.2) and at a higher PaCO2 level from 7.1 to 125.7 ml/100 g brain tissue/min (median = 18.75). The calculated CO2 reactivity ranged from -9.6 to 6.6% (median 1.1%) change in CBF/mm Hg change in PaCO2. CO2 reactivity correlated with lowest pH (r2 = 0.35, p = 0.02). Two infants died, one of neonatal sepsis, the other of heart failure. Neurological outcome at the age of 1 year was normal in 11 patients, 1 had severe cerebral palsy. From the 12 surviving patients the patient with severe neurological deficit showed the highest CBF values (125.7 ml/100 g/min). Impaired chemical coupling of cerebral blood flow is compatible with intact neurological outcome in neonates with perinatal risk factors. CO2 reactivity in these newborns correlates with the lowest pH and may reflect the severity of perinatal asphyxia.

A hybrid algorithm combining metaheuristic with Monte Carlo simulation for solving the Stochastic Flow Shop problem

In this paper, a hybrid simulation-based algorithm is proposed for the StochasticFlow Shop Problem. The main idea of the methodology is to transform the stochastic problem into a deterministic problem and then apply simulation to the latter. In order to achieve this goal, we rely on Monte Carlo Simulation and an adapted version of a deterministic heuristic. This approach aims to provide flexibility and simplicity due to the fact that it is not constrained by any previous assumption and relies in well-tested heuristics.

A hybrid algorithm combining metaheuristic with Monte Carlo simulation for solving the Stochastic Flow Shop problem

In this paper, a hybrid simulation-based algorithm is proposed for the StochasticFlow Shop Problem. The main idea of the methodology is to transform the stochastic problem into a deterministic problem and then apply simulation to the latter. In order to achieve this goal, we rely on Monte Carlo Simulation and an adapted version of a deterministic heuristic. This approach aims to provide flexibility and simplicity due to the fact that it is not constrained by any previous assumption and relies in well-tested heuristics.

Trust region versus line search for computing the optical flow

We consider the numerical treatment of the optical flow problem by evaluating the performance of the trust region method versus the line search method. To the best of our knowledge, the trust region method is studied here for the first time for variational optical flow computation. Four different optical flow models are used to test the performance of the proposed algorithm combining linear and nonlinear data terms with quadratic and TV regularization. We show that trust region often performs better than line search; especially in the presence of non-linearity and non-convexity in the model.

Estimating gene flow across hybrid zones: how reliable are microsatellites?

Due to their high polymorphism, microsatellites have become one of the most valued genetic markers in population biology. We review the first two published studies on hybrid zones of the common shrew based on microsatellites. Both reveal surprisingly high interracial gene flow. It can be shown that these are overestimates. Indeed, in classical population genetics models as F-statistics, mutation is neglected. This constitutes an acceptable assumption as long as migration is higher than mutation. However, in hybrid zones where genetic exchanges can be rare, neglecting mutation will lead to strong overestimates of migration when working with microsatellites which display mutation rates up to 10(-3). As there seems to be no straightforward way to correct for this bias, interracial gene flow estimates based on microsatellites should be taken with caution. This problem should however not conceal the enormous potential of microsatellites to unravel the genetics of hybrid zones.

Milk flow, teat morphology and subclinical mastitis prevalence in Gir cows

The aim of this work was to evaluate the association between milk flow, teat morphological measurements and subclinical mastitis prevalence in Gir cows. Eighty cows in the 2nd and 3rd lactations, with 90 to 200 days of lactation, were divided according to milk flow during milking into fast or slow groups. Teat morphometry was assessed by ultrasound scanning of the right anterior teat and external measurements. Milk samples were collected for somatic cells count (SCC) and microbiological culture. The effect of milk flow during milking was evaluated by analysis of variance of milk yield, SCC, morphometry and external measurements. The association of morphometry and external measurements of the teats with the SCC and microorganisms found in milk were analysed. Milk flow was significantly correlated to milk production. Gir cows with slower milk flow had longer teat canal and greater milk yield, in comparison to cows with fast milk flow. Teat-end to floor distance influenced SCC of Gir cows. Prevalence of subclinical mastitis and the type of mastitis-causing pathogens were not affected by milk flow during milking

A study of the effect of sartorius transposition on lymph flow after ilioinguinal node dissection.

Ilioinguinal dissection is associated with a high rate of lymphatic complications. Prolonged lymph flow causes greatest concern and preventive strategies are needed. A retrospective study of 28 consecutive patients undergoing groin dissection for melanoma metastases was performed to evaluate the influence of sartorius muscle transposition on lymph flow. Modification of the surgical technique with transposition of the sartorius muscle was not associated with reduced drainage time (P = 0.66). A 2-staged approach, with initial sentinel lymph node resection and lymph node dissection in a second operation, however, lead to shortened duration of the lymph flow (P = 0.01). Prolonged lymphorrhea was more frequent in older (P = 0.03), obese (P = 0.02) patients affected by diabetes mellitus (P = 0.03) and hypertension (P = 0.04).

Spatiotemporal adaptive multiscale multiphysics simulations of two-phase flow

We present a spatiotemporal adaptive multiscale algorithm, which is based on the Multiscale Finite Volume method. The algorithm offers a very efficient framework to deal with multiphysics problems and to couple regions with different spatial resolution. We employ the method to simulate two-phase flow through porous media. At the fine scale, we consider a pore-scale description of the flow based on the Volume Of Fluid method. In order to construct a global problem that describes the coarse-scale behavior, the equations are averaged numerically with respect to auxiliary control volumes, and a Darcy-like coarse-scale model is obtained. The space adaptivity is based on the idea that a fine-scale description is only required in the front region, whereas the resolution can be coarsened elsewhere. Temporal adaptivity relies on the fact that the fine-scale and the coarse-scale problems can be solved with different temporal resolution (longer time steps can be used at the coarse scale). By simulating drainage under unstable flow conditions, we show that the method is able to capture the coarse-scale behavior outside the front region and to reproduce complex fluid patterns in the front region.