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.

Hybrid Multiscale Finite Volume method for two-phase flow in porous media

We present a novel hybrid (or multiphysics) algorithm, which couples pore-scale and Darcy descriptions of two-phase flow in porous media. The flow at the pore-scale is described by the Navier?Stokes equations, and the Volume of Fluid (VOF) method is used to model the evolution of the fluid?fluid interface. An extension of the Multiscale Finite Volume (MsFV) method is employed to construct the Darcy-scale problem. First, a set of local interpolators for pressure and velocity is constructed by solving the Navier?Stokes equations; then, a coarse mass-conservation problem is constructed by averaging the pore-scale velocity over the cells of a coarse grid, which act as control volumes; finally, a conservative pore-scale velocity field is reconstructed and used to advect the fluid?fluid interface. The method relies on the localization assumptions used to compute the interpolators (which are quite straightforward extensions of the standard MsFV) and on the postulate that the coarse-scale fluxes are proportional to the coarse-pressure differences. By numerical simulations of two-phase problems, we demonstrate that these assumptions provide hybrid solutions that are in good agreement with reference pore-scale solutions and are able to model the transition from stable to unstable flow regimes. Our hybrid method can naturally take advantage of several adaptive strategies and allows considering pore-scale fluxes only in some regions, while Darcy fluxes are used in the rest of the domain. Moreover, since the method relies on the assumption that the relationship between coarse-scale fluxes and pressure differences is local, it can be used as a numerical tool to investigate the limits of validity of Darcy's law and to understand the link between pore-scale quantities and their corresponding Darcy-scale variables.

A two-phase composite in simple shear: Effective mechanical anisotropy development and localization potential

We present a combined shape and mechanical anisotropy evolution model for a two-phase inclusion-bearing rock subject to large deformation. A single elliptical inclusion embedded in a homogeneous but anisotropic matrix is used to represent a simplified shape evolution enforced on all inclusions. The mechanical anisotropy develops due to the alignment of elongated inclusions. The effective anisotropy is quantified using the differential effective medium (DEM) approach. The model can be run for any deformation path and an arbitrary viscosity ratio between the inclusion and host phase. We focus on the case of simple shear and weak inclusions. The shape evolution of the representative inclusion is largely insensitive to the anisotropy development and to parameter variations in the studied range. An initial hardening stage is observed up to a shear strain of gamma = 1 irrespective of the inclusion fraction. The hardening is followed by a softening stage related to the developing anisotropy and its progressive rotation toward the shear direction. The traction needed to maintain a constant shear rate exhibits a fivefold drop at gamma = 5 in the limiting case of an inviscid inclusion. Numerical simulations show that our analytical model provides a good approximation to the actual evolution of a two-phase inclusion-host composite. However, the inclusions develop complex sigmoidal shapes resulting in the formation of an S-C fabric. We attribute the observed drop in the effective normal viscosity to this structural development. We study the localization potential in a rock column bearing varying fraction of inclusions. In the inviscid inclusion case, a strain jump from gamma = 3 to gamma = 100 is observed for a change of the inclusion fraction from 20% to 33%.

A framework for hybrid simulations of two-phase flow in porous media

Rough a global coarse problem. Although these techniques are usually employed for problems in which the fine-scale processes are described by Darcy's law, they can also be applied to pore-scale simulations and used as a mathematical framework for hybrid methods that couples a Darcy and pore scales. In this work, we consider a pore-scale description of fine-scale processes. The Navier-Stokes equations are numerically solved in the pore geometry to compute the velocity field and obtain generalized permeabilities. In the case of two-phase flow, the dynamics of the phase interface is described by the volume of fluid method with the continuum surface force model. The MsFV method is employed to construct an algorithm that couples a Darcy macro-scale description with a pore-scale description at the fine scale. The hybrid simulations results presented are in good agreement with the fine-scale reference solutions. As the reconstruction of the fine-scale details can be done adaptively, the presented method offers a flexible framework for hybrid modeling.

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.

Yeast vacuoles fragment in an asymmetrical two-phase process with distinct protein requirements.

Yeast vacuoles fragment and fuse in response to environmental conditions, such as changes in osmotic conditions or nutrient availability. Here we analyze osmotically induced vacuole fragmentation by time-lapse microscopy. Small fragmentation products originate directly from the large central vacuole. This happens by asymmetrical scission rather than by consecutive equal divisions. Fragmentation occurs in two distinct phases. Initially, vacuoles shrink and generate deep invaginations that leave behind tubular structures in their vicinity. Already this invagination requires the dynamin-like GTPase Vps1p and the vacuolar proton gradient. Invaginations are stabilized by phosphatidylinositol 3-phosphate (PI(3)P) produced by the phosphoinositide 3-kinase complex II. Subsequently, vesicles pinch off from the tips of the tubular structures in a polarized manner, directly generating fragmentation products of the final size. This phase depends on the production of phosphatidylinositol-3,5-bisphosphate and the Fab1 complex. It is accelerated by the PI(3)P- and phosphatidylinositol 3,5-bisphosphate-binding protein Atg18p. Thus vacuoles fragment in two steps with distinct protein and lipid requirements.

p21 as a transcriptional co-repressor of S-phase and mitotic control genes.

It has been previously described that p21 functions not only as a CDK inhibitor but also as a transcriptional co-repressor in some systems. To investigate the roles of p21 in transcriptional control, we studied the gene expression changes in two human cell systems. Using a human leukemia cell line (K562) with inducible p21 expression and human primary keratinocytes with adenoviral-mediated p21 expression, we carried out microarray-based gene expression profiling. We found that p21 rapidly and strongly repressed the mRNA levels of a number of genes involved in cell cycle and mitosis. One of the most strongly down-regulated genes was CCNE2 (cyclin E2 gene). Mutational analysis in K562 cells showed that the N-terminal region of p21 is required for repression of gene expression of CCNE2 and other genes. Chromatin immunoprecipitation assays indicated that p21 was bound to human CCNE2 and other p21-repressed genes gene in the vicinity of the transcription start site. Moreover, p21 repressed human CCNE2 promoter-luciferase constructs in K562 cells. Bioinformatic analysis revealed that the CDE motif is present in most of the promoters of the p21-regulated genes. Altogether, the results suggest that p21 exerts a repressive effect on a relevant number of genes controlling S phase and mitosis. Thus, p21 activity as inhibitor of cell cycle progression would be mediated not only by the inhibition of CDKs but also by the transcriptional down-regulation of key genes.

Vaniprevir with pegylated interferon alpha-2a and ribavirin in treatment-naïve patients with chronic hepatitis C: a randomized phase II study.

Vaniprevir (MK-7009) is a macrocyclic hepatitis C virus (HCV) nonstructural protein 3/4A protease inhibitor. The aim of the present phase II study was to examine virologic response rates with vaniprevir in combination with pegylated interferon alpha-2a (Peg-IFN-α-2a) plus ribavirin (RBV). In this double-blind, placebo-controlled, dose-ranging study, treatment-naïve patients with HCV genotype 1 infection (n = 94) were randomized to receive open-label Peg-IFN-α-2a (180 μg/week) and RBV (1,000-1,200 mg/day) in combination with blinded placebo or vaniprevir (300 mg twice-daily [BID], 600 mg BID, 600 mg once-daily [QD], or 800 mg QD) for 28 days, then open-label Peg-IFN-α-2a and RBV for an additional 44 weeks. The primary efficacy endpoint was rapid viral response (RVR), defined as undetectable plasma HCV RNA at week 4. Across all doses, vaniprevir was associated with a rapid two-phase decline in viral load, with HCV RNA levels approximately 3 log(10) IU/mL lower in vaniprevir-treated patients, compared to placebo recipients. Rates of RVR were significantly higher in each of the vaniprevir dose groups, compared to the control regimen (68.8%-83.3% versus 5.6%; P < 0.001 for all comparisons). There were numerically higher, but not statistically significant, early and sustained virologic response rates with vaniprevir, as compared to placebo. Resistance profile was predictable, with variants at R155 and D168 detected in a small number of patients. No relationship between interleukin-28B genotype and treatment outcomes was demonstrated in this study. The incidence of adverse events was generally comparable between vaniprevir and placebo recipients; however, vomiting appeared to be more common at higher vaniprevir doses. CONCLUSION: Vaniprevir is a potent HCV protease inhibitor with a predictable resistance profile and favorable safety profile that is suitable for QD or BID administration.

Towards the validation of "in silico" models and physicochemical filters to identify and characterize new chemical entities

Summary: Lipophilicity plays an important role in the determination and the comprehension of the pharmacokinetic behavior of drugs. It is usually expressed by the partition coefficient (log P) in the n-octanol/water system. The use of an additional solvent system (1,2-dichlorethane/water) is necessary to obtain complementary information, as the log Poct values alone are not sufficient to explain ail biological properties. The aim of this thesis is to develop tools allowing to predict lipophilicity of new drugs and to analyze the information yielded by those log P values. Part I presents the development of theoretical models used to predict lipophilicity. Chapter 2 shows the necessity to extend the existing solvatochromic analyses in order to predict correctly the lipophilicity of new and complex neutral compounds. In Chapter 3, solvatochromic analyses are used to develop a model for the prediction of the lipophilicity of ions. A global model was obtained allowing to estimate the lipophilicity of neutral, anionic and cationic solutes. Part II presents the detailed study of two physicochemical filters. Chapter 4 shows that the Discovery RP Amide C16 stationary phase allows to estimate lipophilicity of the neutral form of basic and acidic solutes, except of lipophilic acidic solutes. Those solutes present additional interactions with this particular stationary phase. In Chapter 5, 4 different IANI stationary phases are investigated. For neutral solutes, linear data are obtained whatever the IANI column used. For the ionized solutes, their retention is due to a balance of electrostatic and hydrophobie interactions. Thus no discrimination is observed between different series of solutes bearing the same charge, from one column to an other. Part III presents two examples illustrating the information obtained thanks to Structure-Properties Relationships (SPR). Comparing graphically lipophilicity values obtained in two different solvent systems allows to reveal the presence of intramolecular effects .such as internai H-bond (Chapter 6). SPR is used to study the partitioning of ionizable groups encountered in Medicinal Chemistry (Chapter7). Résumé La lipophilie joue un .rôle important dans la détermination et la compréhension du comportement pharmacocinétique des médicaments. Elle est généralement exprimée par le coefficient de partage (log P) d'un composé dans le système de solvants n-octanol/eau. L'utilisation d'un deuxième système de solvants (1,2-dichloroéthane/eau) s'est avérée nécessaire afin d'obtenir des informations complémentaires, les valeurs de log Poct seules n'étant pas suffisantes pour expliquer toutes les propriétés biologiques. Le but de cette thèse est de développer des outils permettant de prédire la lipophilie de nouveaux candidats médicaments et d'analyser l'information fournie par les valeurs de log P. La Partie I présente le développement de modèles théoriques utilisés pour prédire la lipophilie. Le chapitre 2 montre la nécessité de mettre à jour les analyses solvatochromiques existantes mais inadaptées à la prédiction de la lipophilie de nouveaux composés neutres. Dans le chapitre 3, la même méthodologie des analyses solvatochromiques est utilisée pour développer un modèle permettant de prédire la lipophilie des ions. Le modèle global obtenu permet la prédiction de la lipophilie de composés neutres, anioniques et cationiques. La Partie II présente l'étude approfondie de deux filtres physicochimiques. Le Chapitre 4 montre que la phase stationnaire Discovery RP Amide C16 permet la détermination de la lipophilie de la forme neutre de composés basiques et acides, à l'exception des acides très lipophiles. Ces derniers présentent des interactions supplémentaires avec cette phase stationnaire. Dans le Chapitre 5, 4 phases stationnaires IAM sont étudiées. Pour les composés neutres étudiés, des valeurs de rétention linéaires sont obtenues, quelque que soit la colonne IAM utilisée. Pour les composés ionisables, leur rétention est due à une balance entre des interactions électrostatiques et hydrophobes. Donc aucune discrimination n'est observée entre les différentes séries de composés portant la même charge d'une colonne à l'autre. La Partie III présente deux exemples illustrant les informations obtenues par l'utilisation des relations structures-propriétés. Comparer graphiquement la lipophilie mesurée dans deux différents systèmes de solvants permet de mettre en évidence la présence d'effets intramoléculaires tels que les liaisons hydrogène intramoléculaires (Chapitre 6). Cette approche des relations structures-propriétés est aussi appliquée à l'étude du partage de fonctions ionisables rencontrées en Chimie Thérapeutique (Chapitre 7) Résumé large public Pour exercer son effet thérapeutique, un médicament doit atteindre son site d'action en quantité suffisante. La quantité effective de médicament atteignant le site d'action dépend du nombre d'interactions entre le médicament et de nombreux constituants de l'organisme comme, par exemple, les enzymes du métabolisme ou les membranes biologiques. Le passage du médicament à travers ces membranes, appelé perméation, est un paramètre important à optimiser pour développer des médicaments plus puissants. La lipophilie joue un rôle clé dans la compréhension de la perméation passive des médicaments. La lipophilie est généralement exprimée par le coefficient de partage (log P) dans le système de solvants (non miscibles) n-octanol/eau. Les valeurs de log Poct seules se sont avérées insuffisantes pour expliquer la perméation à travers toutes les différentes membranes biologiques du corps humain. L'utilisation d'un système de solvants additionnel (le système 1,2-dichloroéthane/eau) a permis d'obtenir les informations complémentaires indispensables à une bonne compréhension du processus de perméation. Un grand nombre d'outils expérimentaux et théoriques sont à disposition pour étudier la lipophilie. Ce travail de thèse se focalise principalement sur le développement ou l'amélioration de certains de ces outils pour permettre leur application à un champ plus large de composés. Voici une brève description de deux de ces outils: 1)La factorisation de la lipophilie en fonction de certaines propriétés structurelles (telle que le volume) propres aux composés permet de développer des modèles théoriques utilisables pour la prédiction de la lipophilie de nouveaux composés ou médicaments. Cette approche est appliquée à l'analyse de la lipophilie de composés neutres ainsi qu'à la lipophilie de composés chargés. 2)La chromatographie liquide à haute pression sur phase inverse (RP-HPLC) est une méthode couramment utilisée pour la détermination expérimentale des valeurs de log Poct.

X-ray mosaic nanotomography of large microorganisms.

Full-field X-ray microscopy is a valuable tool for 3D observation of biological systems. In the soft X-ray domain organelles can be visualized in individual cells while hard X-ray microscopes excel in imaging of larger complex biological tissue. The field of view of these instruments is typically 10(3) times the spatial resolution. We exploit the assets of the hard X-ray sub-micrometer imaging and extend the standard approach by widening the effective field of view to match the size of the sample. We show that global tomography of biological systems exceeding several times the field of view is feasible also at the nanoscale with moderate radiation dose. We address the performance issues and limitations of the TOMCAT full-field microscope and more generally for Zernike phase contrast imaging. Two biologically relevant systems were investigated. The first being the largest known bacteria (Thiomargarita namibiensis), the second is a small myriapod species (Pauropoda sp.). Both examples illustrate the capacity of the unique, structured condenser based broad-band full-field microscope to access the 3D structural details of biological systems at the nanoscale while avoiding complicated sample preparation, or even keeping the sample environment close to the natural state.

A pea plasma membrane protein exhibiting blue light-induced phosphorylation retains photosensitivity following triton solubilization.

Phosphorylation of a polypeptide of approximately 120 kD in pea (Pisum sativum L.) plasma membranes in response to blue light has been shown to be involved in phototropic curvature, but the relationship of this protein to the kinase and photoreceptor acting upon it is uncertain. Using two-phase aqueous partitioning to isolate right-side-out plasma membrane vesicles, we have obtained evidence suggesting that the photoreceptor, kinase, and substrate are localized to the plasma membrane fraction. Latent phosphorylation accessible through Triton X-100 or freeze/thaw treatments of purified plasma membrane vesicles indicates that at least the kinase moiety is present on the internal face of the plasma membrane. Effects of solubilization of vesicles on fluence-response characteristics and on phosphorylation levels provide evidence that the receptor, kinase, and protein substrate are present together in individual mixed detergent micelles, either as a stable complex or as domains of a single polypeptide. In vivo blue-light irradiation results in a small but significant decrease in mobility of the 120-kD phosphorylated protein on sodium dodecylsulfate gel electrophoresis. This mobility shift is evident on Coomassie-stained gels and on western blots probed with polyclonal antibodies raised against the 120-kD protein. Among the plasma membrane proteins bound to the reactive nucleotide analog fluorosulfonylbenzoyladenine (FSBA), a distinct protein band at 120 kD can be detected on blots probed with anti-FSBA antibodies. This band exhibits an in vivo light-dependent mobility shift identical to that observed for the protein band and antibodies specific for the 120-kD protein, implying that the 120-kD protein has an integral nucleotide binding site and consistent with the possibility that the substrate protein is also a kinase.

Structural and metamorphic evolution of the northern Himachal Himalaya, NW India - (Spiti-eastern Lahul-Parvati valley traverse)

The Himalayan orogen is the result of the collision between the Indian and Asian continents that began 55-50 Ma ago, causing intracontinental thrusting and nappe formation. Detailed mapping as well as structural and microfabric analyses on a traverse from the Tethyan Himalaya southwestward through the High Himalayan Crystalline and the Main Central Thrust zone (MCT zone) to the Lesser Himalayan Sequence in the Spiti-eastern Lahul-Parvati valley area reveal eight main phases of deformation, a series of late stage phases and five stages of metamorphic crystallization. This sequence of events is integrated into a reconstruction of the tectonometamorphic evolution of the Himalayan orogen in northern Himachal Pradesh. The oldest phase D-1 is preserved as relies in the High Himalayan Crystalline. Its deformational conditions are poorly known, but the metamorphic evolution is well documented by a prograde metamorphism reaching peak conditions within the upper amphibolite facies. This indicates that D-1 was an important tectonometamorphic event including considerable crustal thickening. The structural, metamorphic and sedimentary record suggest that D-1 most probably represents an early stage of continental collision. The first event clearly attributed to the collision between India and Asia is documented by two converging nappe systems, the NE-verging Shikar Beh Nappe and the SW-verging north Himalayan nappes. The D-2 Shikar Beh Nappe is characterized by isoclinal folding and top-to-the NE shearing, representing the main deformation in the High Himalayan Crystalline. D-2 also caused the main metamorphism in the High Himalayan Crystalline that was of a Barrovian-type, reaching upper amphibolite facies peak conditions. The Shikar Beh Nappe is interpreted to have formed within the Indian crust SW of the subduction zone. Simultaneously with NE-directed nappe formation, incipient subduction of India below Asia caused stacking of the SW-verging north Himalayan Nappes, that were thrust from the northern edge of the subducted continent toward the front of the Shikar Beh Nappe. As a result, the SW-verging folds of the D-3 Main Fold Zone formed in the Tethyan Himalaya below the front of the north Himalayan nappes. D-3 represents the main deformation in the Tethyan Himalaya, associated with a greenschist facies metamorphism. Folding within the Main Fold Zone subsequently propagated toward SW into the High Himalayan Crystalline, where it overprinted the preexisting D-2 structures. After subduction at the base of the north Himalayan nappes, the subduction zone stepped to the base of the High Himalayan Crystalline, where D-3 folds were crosscut by SW-directed D-4 thrusting. During D-4, the Crystalline Nappe, comprising the Main Fold Zone and relies of the Shikar Beh Nappe was thrust toward SW over the Lesser Himalayan Sequence along the 4 to 5 kms thick Main Central Thrust zone. Thrusting was related to a retrograde greenschist facies overprint at the base of the Crystalline Nappe and to pro-grade greenschist facies conditions in the Lesser Himalayan Sequence. Simultaneously with thrusting at the base of the Crystalline Nappe, higher crustal levels were affected by NE-directed D-5 normal extensional shearing and by dextral strike-slip motion, indicating that the high-grade metamorphic Crystalline Nappe was extruded between the low-grade metamorphic Lesser Himalayan Sequence at the base and the north Himalayan nappes at the top. The upper boundary of the Crystalline Nappe is not clearly delimited and passes gradually into the low-grade rocks at the front of the north Himalayan nappes. Extrusion of the Crystalline Nappe was followed by the phase D-6, characterized by large-scale, upright to steeply inclined, NE-verging folds and by another series of normal and extensional structures D-7+D-8 that may be related to ongoing extrusion of the Crystalline Nappe. The late stage evolution is represented by the phases D-A and D-B that indicate shortening parallel to the axis of the mountain chain and by D-C that is interpreted to account for the formation of large-scale domes with NNW-SSE-trending axes, an example of which is exposed in the Larji-Kullu-Rampur tectonic window.

Two randomized trials of canakinumab in systemic juvenile idiopathic arthritis.

BACKGROUND: Interleukin-1 is pivotal in the pathogenesis of systemic juvenile idiopathic arthritis (JIA). We assessed the efficacy and safety of canakinumab, a selective, fully human, anti-interleukin-1β monoclonal antibody, in two trials. METHODS: In trial 1, we randomly assigned patients, 2 to 19 years of age, with systemic JIA and active systemic features (fever; ≥2 active joints; C-reactive protein, >30 mg per liter; and glucocorticoid dose, ≤1.0 mg per kilogram of body weight per day), in a double-blind fashion, to a single subcutaneous dose of canakinumab (4 mg per kilogram) or placebo. The primary outcome, termed adapted JIA ACR 30 response, was defined as improvement of 30% or more in at least three of the six core criteria for JIA, worsening of more than 30% in no more than one of the criteria, and resolution of fever. In trial 2, after 32 weeks of open-label treatment with canakinumab, patients who had a response and underwent glucocorticoid tapering were randomly assigned to continued treatment with canakinumab or to placebo. The primary outcome was time to flare of systemic JIA. RESULTS: At day 15 in trial 1, more patients in the canakinumab group had an adapted JIA ACR 30 response (36 of 43 [84%], vs. 4 of 41 [10%] in the placebo group; P<0.001). In trial 2, among the 100 patients (of 177 in the open-label phase) who underwent randomization in the withdrawal phase, the risk of flare was lower among patients who continued to receive canakinumab than among those who were switched to placebo (74% of patients in the canakinumab group had no flare, vs. 25% in the placebo group, according to Kaplan-Meier estimates; hazard ratio, 0.36; P=0.003). The average glucocorticoid dose was reduced from 0.34 to 0.05 mg per kilogram per day, and glucocorticoids were discontinued in 42 of 128 patients (33%). The macrophage activation syndrome occurred in 7 patients; infections were more frequent with canakinumab than with placebo. CONCLUSIONS: These two phase 3 studies show the efficacy of canakinumab in systemic JIA with active systemic features. (Funded by Novartis Pharma; ClinicalTrials.gov numbers, NCT00889863 and NCT00886769.).

Two randomized trials of canakinumab in systemic juvenile idiopathic arthritis.

BACKGROUND: Interleukin-1 is pivotal in the pathogenesis of systemic juvenile idiopathic arthritis (JIA). We assessed the efficacy and safety of canakinumab, a selective, fully human, anti-interleukin-1β monoclonal antibody, in two trials. METHODS: In trial 1, we randomly assigned patients, 2 to 19 years of age, with systemic JIA and active systemic features (fever; ≥2 active joints; C-reactive protein, >30 mg per liter; and glucocorticoid dose, ≤1.0 mg per kilogram of body weight per day), in a double-blind fashion, to a single subcutaneous dose of canakinumab (4 mg per kilogram) or placebo. The primary outcome, termed adapted JIA ACR 30 response, was defined as improvement of 30% or more in at least three of the six core criteria for JIA, worsening of more than 30% in no more than one of the criteria, and resolution of fever. In trial 2, after 32 weeks of open-label treatment with canakinumab, patients who had a response and underwent glucocorticoid tapering were randomly assigned to continued treatment with canakinumab or to placebo. The primary outcome was time to flare of systemic JIA. RESULTS: At day 15 in trial 1, more patients in the canakinumab group had an adapted JIA ACR 30 response (36 of 43 [84%], vs. 4 of 41 [10%] in the placebo group; P<0.001). In trial 2, among the 100 patients (of 177 in the open-label phase) who underwent randomization in the withdrawal phase, the risk of flare was lower among patients who continued to receive canakinumab than among those who were switched to placebo (74% of patients in the canakinumab group had no flare, vs. 25% in the placebo group, according to Kaplan-Meier estimates; hazard ratio, 0.36; P=0.003). The average glucocorticoid dose was reduced from 0.34 to 0.05 mg per kilogram per day, and glucocorticoids were discontinued in 42 of 128 patients (33%). The macrophage activation syndrome occurred in 7 patients; infections were more frequent with canakinumab than with placebo. CONCLUSIONS: These two phase 3 studies show the efficacy of canakinumab in systemic JIA with active systemic features. (Funded by Novartis Pharma; ClinicalTrials.gov numbers, NCT00889863 and NCT00886769.).