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.

Hyperventilation in anticipatory music performance anxiety

Objectives and Methods: Self-report studies have shown an association between music performance anxiety (MPA) and hyperventilation complaints. However, hyperventilation was never assessed physiologically in MPA. This study investigated the self-reported affective experience, self-reported physiological symptoms, and cardiorespiratory variables including partial pressure of end-tidal CO(2) (Petco(2)), which is an indicator for hyperventilation, in 67 music students before a private and a public performance. The response coherence between these response domains was also investigated.ResultsFrom the private to the public session, the intensity of all self-report variables increased (all p values < .001). As predicted, the higher the musician's usual MPA level, the larger were these increases (p values < .10). With the exception of Petco(2), the main cardiorespiratory variables also increased from the private to the public session (p values < .05). These increases were not modulated by the usual MPA level (p values > .10). Petco(2) showed a unique response pattern reflected by an MPA-by-session interaction (p < .01): it increased from the private to the public session for musicians with low MPA levels and decreased for musicians with high MPA levels. Self-reported physiological symptoms were related to the self-reported affective experience (p values < .05) rather than to physiological measures (p values > .17).ConclusionsThese findings show for the first time how respiration is stimulated before a public performance in music students with different MPA levels. The hypothesis of a hyperventilation tendency in high-performance-anxious musicians is supported. The response coherence between physiological symptoms and physiological activation is weak.

A multi-perspective approach to performance anxiety in music students

AbstractPerforming publicly has become increasingly important in a variety of professions. This condition is associated with performance anxiety in almost all performers. Whereas some performers successfully cope with this anxiety, for others it represents a major problem and even threatens their career. Musicians and especially music students were shown to be particularly affected by performance anxiety.Therefore, the goal of this PhD thesis was to gain a better understanding of performance anxiety in university music students. More precisely, the first part of this thesis aimed at increasing knowledge on the occurrence, the experience, and the management of performance anxiety (Article 1). The second part aimed at investigating the hypothesis that there is an underlying hyperventilation problem in musicians with a high level of anxiety before a performance. This hypothesis was addressed in two ways: firstly, by investigating the association between the negative affective dimension of music performance anxiety (MPA) and self-perceived physiological symptoms that are known to co-occur with hyperventilation (Article 2) and secondly, by analyzing this association on the physiological level before a private (audience-free) and a public performance (Article 3). Article 4 places some key variables of Article 3 in a larger context by jointly analyzing the phases before, during, and after performing.The main results of the self-report data show (a) that stage fright is experienced as a problem by one-third of the surveyed students, (b) that the students express a considerable need for more help to better cope with it, and (c) that there is a positive association between negative feelings of MPA and the self-reported hyperventilation complaints before performing. This latter finding was confirmed on the physiological level in a tendency of particularly high performance-anxious musicians to hyperventilate. Furthermore, the psycho-physiological activation increased from a private to a public performance, and was higher during the performances than before or after them. The physiological activation was mainly independent of the MPA score. Finally, there was a low response coherence between the actual physiological activation and the self-reports on the instantaneous anxiety, tension, and perceived physiological activation.Given the high proportion of music students who consider stage fright as a problem and given the need for more help to better cope with it, a better understanding of this phenomenon and its inclusion in the educational process is fundamental to prevent future occupational problems. On the physiological level, breathing exercises might be a good means to decrease - but also to increase - the arousal associated with a public performance in order to meet an optimal level of arousal needed for a good performance.

High cefepime plasma concentrations and neurological toxicity in febrile neutropenic patients with mild impairment of renal function.

High-dose cefepime therapy is recommended for febrile neutropenia. Safety issues have been raised in a recent meta-analysis reporting an increased risk of mortality during cefepime therapy. Cefepime-related neurological toxicity has been associated with overdosing due to severe renal dysfunction. This study aimed to investigate the association between cefepime plasma concentrations and neurological toxicity in febrile neutropenic patients. Cefepime trough concentrations (by high-performance liquid chromatography) were retrospectively analyzed for 30 adult febrile neutropenic patients receiving the recommended high-dose regimen (6 g/day for a glomerular filtration rate [GFR] of >50 ml/min). The dose adjustment to renal function was evaluated by the ratio of the cefepime daily dose per 100 ml/min of glomerular filtration. The association between cefepime plasma concentrations and neurological toxicity was assessed on the basis of consistent neurological symptoms and/or signs (by NCI criteria). The median cefepime concentration was 8.7 mg/liter (range, 2.1 to 38 mg/liter) at a median of 4 days (range, 2 to 15 days) after the start of therapy. Neurological toxicity (altered mental status, hallucinations, or myoclonia) was attributed to cefepime in 6/30 (20%) patients (median GFR, 45 ml/min; range, 41 to 65 ml/min) receiving a median dose of 13.2 g/day per 100 ml/min GFR (range, 9.2 to 14.3 g/day per 100 ml/min GFR). Cefepime discontinuation resulted in complete neurological recovery for five patients and improvement for one patient. A multivariate logistic regression model confirmed high cefepime concentrations as an independent predictor of neurological toxicity, with a 50% probability threshold at ≥22 mg/liter (P = 0.05). High cefepime plasma concentrations are associated with neurological toxicity in febrile neutropenic patients with mild renal dysfunction. Careful adherence to normalized dosing per 100 ml/min GFR is crucial. Monitoring of plasma concentrations may contribute to preventing neurological toxicity of high-dose therapy for this life-threatening condition.

Performance of blood pressure-to-height ratio at a single screening visit for the identification of hypertension in children.

BACKGROUND: The diagnosis of hypertension in children is difficult because of the multiple sex-, age-, and height-specific thresholds to define elevated blood pressure (BP). Blood pressure-to-height ratio (BPHR) has been proposed to facilitate the identification of elevated BP in children. OBJECTIVE: We assessed the performance of BPHR at a single screening visit to identify children with hypertension that is sustained elevated BP. METHOD: In a school-based study conducted in Switzerland, BP was measured at up to three visits in 5207 children. Children had hypertension if BP was elevated at the three visits. Sensitivity, specificity, negative predictive value (NPV), and positive predictive value (PPV) for the identification of hypertension were assessed for different thresholds of BPHR. The ability of BPHR at a single screening visit to discriminate children with and without hypertension was evaluated with receiver operating characteristic (ROC) curve analyses. RESULTS: The prevalence of systolic/diastolic hypertension was 2.2%. Systolic BPHR had a better performance to identify hypertension compared with diastolic BPHR (area under the ROC curve: 0.95 vs. 0.84). The highest performance was obtained with a systolic BPHR threshold set at 0.80 mmHg/cm (sensitivity: 98%; specificity: 85%; PPV: 12%; and NPV: 100%) and a diastolic BPHR threshold set at 0.45 mmHg/cm (sensitivity: 79%; specificity: 70%; PPV: 5%; and NPV: 99%). The PPV was higher among tall or overweight children. CONCLUSION: BPHR at a single screening visit had a high performance to identify hypertension in children, although the low prevalence of hypertension led to a low PPV.

SlimCodeML: An Optimized Version of CodeML for the Branch-Site Model

CodeML (part of the PAML package) im- plements a maximum likelihood-based approach to de- tect positive selection on a specific branch of a given phylogenetic tree. While CodeML is widely used, it is very compute-intensive. We present SlimCodeML, an optimized version of CodeML for the branch-site model. Our performance analysis shows that SlimCodeML substantially outperforms CodeML (up to 9.38 times faster), especially for large-scale genomic analyses.

Three essays on major trends in a slow clockspeed industry : the case of industrial automation

The motivation for this research initiated from the abrupt rise and fall of minicomputers which were initially used both for industrial automation and business applications due to their significantly lower cost than their predecessors, the mainframes. Later industrial automation developed its own vertically integrated hardware and software to address the application needs of uninterrupted operations, real-time control and resilience to harsh environmental conditions. This has led to the creation of an independent industry, namely industrial automation used in PLC, DCS, SCADA and robot control systems. This industry employs today over 200'000 people in a profitable slow clockspeed context in contrast to the two mainstream computing industries of information technology (IT) focused on business applications and telecommunications focused on communications networks and hand-held devices. Already in 1990s it was foreseen that IT and communication would merge into one Information and communication industry (ICT). The fundamental question of the thesis is: Could industrial automation leverage a common technology platform with the newly formed ICT industry? Computer systems dominated by complex instruction set computers (CISC) were challenged during 1990s with higher performance reduced instruction set computers (RISC). RISC started to evolve parallel to the constant advancement of Moore's law. These developments created the high performance and low energy consumption System-on-Chip architecture (SoC). Unlike to the CISC processors RISC processor architecture is a separate industry from the RISC chip manufacturing industry. It also has several hardware independent software platforms consisting of integrated operating system, development environment, user interface and application market which enables customers to have more choices due to hardware independent real time capable software applications. An architecture disruption merged and the smartphone and tablet market were formed with new rules and new key players in the ICT industry. Today there are more RISC computer systems running Linux (or other Unix variants) than any other computer system. The astonishing rise of SoC based technologies and related software platforms in smartphones created in unit terms the largest installed base ever seen in the history of computers and is now being further extended by tablets. An underlying additional element of this transition is the increasing role of open source technologies both in software and hardware. This has driven the microprocessor based personal computer industry with few dominating closed operating system platforms into a steep decline. A significant factor in this process has been the separation of processor architecture and processor chip production and operating systems and application development platforms merger into integrated software platforms with proprietary application markets. Furthermore the pay-by-click marketing has changed the way applications development is compensated: Three essays on major trends in a slow clockspeed industry: The case of industrial automation 2014 freeware, ad based or licensed - all at a lower price and used by a wider customer base than ever before. Moreover, the concept of software maintenance contract is very remote in the app world. However, as a slow clockspeed industry, industrial automation has remained intact during the disruptions based on SoC and related software platforms in the ICT industries. Industrial automation incumbents continue to supply systems based on vertically integrated systems consisting of proprietary software and proprietary mainly microprocessor based hardware. They enjoy admirable profitability levels on a very narrow customer base due to strong technology-enabled customer lock-in and customers' high risk leverage as their production is dependent on fault-free operation of the industrial automation systems. When will this balance of power be disrupted? The thesis suggests how industrial automation could join the mainstream ICT industry and create an information, communication and automation (ICAT) industry. Lately the Internet of Things (loT) and weightless networks, a new standard leveraging frequency channels earlier occupied by TV broadcasting, have gradually started to change the rigid world of Machine to Machine (M2M) interaction. It is foreseeable that enough momentum will be created that the industrial automation market will in due course face an architecture disruption empowered by these new trends. This thesis examines the current state of industrial automation subject to the competition between the incumbents firstly through a research on cost competitiveness efforts in captive outsourcing of engineering, research and development and secondly researching process re- engineering in the case of complex system global software support. Thirdly we investigate the industry actors', namely customers, incumbents and newcomers, views on the future direction of industrial automation and conclude with our assessments of the possible routes industrial automation could advance taking into account the looming rise of the Internet of Things (loT) and weightless networks. Industrial automation is an industry dominated by a handful of global players each of them focusing on maintaining their own proprietary solutions. The rise of de facto standards like IBM PC, Unix and Linux and SoC leveraged by IBM, Compaq, Dell, HP, ARM, Apple, Google, Samsung and others have created new markets of personal computers, smartphone and tablets and will eventually also impact industrial automation through game changing commoditization and related control point and business model changes. This trend will inevitably continue, but the transition to a commoditized industrial automation will not happen in the near future.

Photodynamic therapy with mTHPC and polyethylene glycol-derived mTHPC: a comparative study on human tumour xenografts.

The photosensitizing properties of m-tetrahydroxyphenylchlorin (mTHPC) and polyethylene glycol-derivatized mTHPC (pegylated mTHPC) were compared in nude mice bearing human malignant mesothelioma, squamous cell carcinoma and adenocarcinoma xenografts. Laser light (20 J/cm2) at 652 nm was delivered to the tumour (surface irradiance) and to an equal-sized area of the hind leg of the animals after i.p. administration of 0.1 mg/kg body weight mTHPC and an equimolar dose of pegylated mTHPC, respectively. The extent of tumour necrosis and normal tissue injury was assessed by histology. Both mTHPC and pegylated mTHPC catalyse photosensitized necrosis in mesothelioma xenografts at drug-light intervals of 1-4 days. The onset of action of pegylated mTHPC seemed slower but significantly exceeds that of mTHPC by days 3 and 4 with the greatest difference being noted at day 4. Pegylated mTHPC also induced significantly larger photonecrosis than mTHPC in squamous cell xenografts but not in adenocarcinoma at day 4, where mTHPC showed greatest activity. The degree of necrosis induced by pegylated mTHPC was the same for all three xenografts. mTHPC led to necrosis of skin and underlying muscle at a drug-light interval of 1 day but minor histological changes only at drug-light intervals from 2-4 days. In contrast, pegylated mTHPC did not result in histologically detectable changes in normal tissues under the same treatment conditions at any drug-light interval assessed. In this study, pegylated mTHPC had advantages as a photosensitizer compared to mTHPC. Tissue concentrations of mTHPC and pegylated mTHPC were measured by high-performance liquid chromatography in non-irradiated animals 4 days after administration. There was no significant difference in tumour uptake between the two sensitizers in mesothelioma, adenocarcinoma and squamous cell carcinoma xenografts. Tissue concentration measurements were of limited use for predicting photosensitization in this model.

Photodynamic therapy enhances liposomal doxorubicin distribution in tumors during isolated perfusion of rodent lungs.

BACKGROUND: Photodynamic therapy (PDT) at low drug-light conditions can enhance the transport of intravenously injected macromolecular therapeutics through the tumor vasculature. Here we determined the impact of PDT on the distribution of liposomal doxorubicin (Liporubicin™) administered by isolated lung perfusion (ILP) in sarcomas grown on rodent lungs. METHODS: A syngeneic methylcholanthrene-induced sarcoma cell line was implanted subpleurally in the left lung of Fischer rats. Treatment schemes consisted in ILP alone (400 μg of Liporubicin), low-dose (0.0625 mg/kg Visudyne®, 10 J/cm(2) and 35 mW/cm(2)) and high-dose left lung PDT (0.125 mg/kg Visudyne, 10 J/cm(2) and 35 mW/cm(2)) followed by ILP (400 μg of Liporubicin). The uptake and distribution of Liporubicin in tumor and lung tissues were determined by high-performance liquid chromatography and fluorescence microscopy in each group. RESULTS: Low-dose PDT significantly improved the distribution of Liporubicin in tumors compared to high-dose PDT (p < 0.05) and ILP alone (p < 0.05). However, both PDT pretreatments did not result in a higher overall drug uptake in tumors or a higher tumor-to-lung drug ratio compared to ILP alone. CONCLUSIONS: Intraoperative low-dose Visudyne-mediated PDT enhances liposomal doxorubicin distribution administered by ILP in sarcomas grown on rodent lungs which is predicted to improve tumor control by ILP.

The mouse Lyt-2/3 antigen complex--II. Structural analysis of the subunits.

Surface- or biosynthetically labeled Lyt-2/3 antigens were isolated from cell lysates by immunoprecipitation and affinity chromatography with a monoclonal antibody. Tryptic digests of the individual subunits of 37,000, 32,000 and 28,000 apparent mol. wts were analysed by reverse-phase high-performance liquid chromatography and by two-dimensional peptide mapping. The results indicate that the 37,000 and 32,000 mol. wt components are structurally very similar whereas the 28,000 mol. wt component appears as a different molecule.

Design and evaluation of a solid sampler for the monitoring of airborne 1,6-hexamethylene diisocyanate (HDI) and its prepolymers in 2-component spray painting

An active, solvent-free solid sampler was developed for the collection of 1,6-hexamethylene diisocyanate (HDI) aerosol and prepolymers. The sampler was made of a filter impregnated with 1-(2-methoxyphenyl)piperazine contained in a filter holder. Interferences with HDI were observed when a set of cellulose acetate filters and a polystyrene filter holder were used; a glass fiber filter and polypropylene filter cassette gave better results. The applicability of the sampling and analytical procedure was validated with a test chamber, constructed for the dynamic generation of HDI aerosol and prepolymers in commercial two-component spray paints (Desmodur(R) N75) used in car refinishing. The particle size distribution, temporal stability, and spatial uniformity of the simulated aerosol were established in order to test the sample. The monitoring of aerosol concentrations was conducted with the solid sampler paired to the reference impinger technique (impinger flasks contained 10 mL of 0.5 mg/mL 1-(2-methoxyphenyl)piperazine in toluene) under a controlled atmosphere in the test chamber. Analyses of derivatized HDI and prepolymers were carried out by using high-performance liquid chromatography and ultraviolet detection. The correlation between the solvent-free and the impinger techniques appeared fairly good (Y = 0.979X - 0.161; R = 0.978), when the tests were conducted in the range of 0.1 to 10 times the threshold limit value (TLV) for HDI monomer and up to 60-mu-g/m3 (3 U.K. TLVs) for total -N = C = O groups.

Plasma levels of trimipramine and metabolites in four patients: determination of the enantiomer concentrations of the hydroxy metabolites.

A two-step high-performance liquid chromatography method is described, using a CN column and an alpha 1-acid glycoprotein column, which allows the measurement of the enantiomers of the hydroxy metabolites of trimipramine in plasma of trimipramine-treated patients. Of the four patients analyzed, three showed approximately equimolar concentrations of the (D)- and (L)-enantiomers of the hydroxy metabolites (2-hydroxy-trimipramine and 2-hydroxy desmethyltrimipramine), and one was found to have roughly twice as much of the (L)-form and of the (D)-form of 2-hydroxy trimipramine and 2-hydroxy desmethyltrimipramine. From the data available on the pharmacological effects of the enantiomers of trimipramine, it is postulated that this interindividual variability in its pharmacokinetics is another factor that could contribute to the interindividual variability in its pharmacodynamics.