In vivo study comparing an X-shaped polymethylmethacrylate and a cylindrical collagen implant for deep sclerectomy.

Background:  Study in vivo characteristics of a polymethylmethacrylate (PMMA) implant compared to the standard cylindrical collagen implant for deep sclerectomy (DS). Design:  Six-month comparative study. Samples:  Twenty eyes of ten rabbits. Methods:  Eyes were randomized to have DS with PMMA implant in one eye and collagen implant in the opposite eye. The growth of the new subconjunctival drainage vessels was assessed by combined fluorescein and indocyanin green anterior segment angiography; intrascleral and subconjunctival blebs were imaged by ultrasound biomicroscopy (UBM). At six months, outflow facility (C) was measured by anterior chamber perfusion and portions of one side of the DS were compared to portions on the 180° opposite side and native sclera on histology. Results:  The mean IOP preoperatively and at one, four, twelve, and twenty-four weeks was comparable in both groups (P > 0.1). UBM showed a statistically insignificant quicker regression of the subconjunctival bleb as well as a durable intrascleral lake in the PMMA group (P > 0.05). New drainage vessels were initially observed one month after surgery; they were more numerous in the PMMA group on angiographic and histological findings at 6 months (P < 0.05). The mean C increased significantly after surgery compared to preoperative values (P < 0.05) and no difference was observed between the implants (0.24 ± 0.06 µl/min/mmHg [PMMA] and 0.23 ± 0.07 µl/min/mmHg [collagen implant]) (P = 0.39). Conclusions:  Deep sclerectomy performed with PMMA or collagen implants showed similar IOP lowering effects, outflow facility increase, and degree of inflammatory reaction.

In vivo study comparing an X-shaped polymethylmethacrylate and a cylindrical collagen implant for deep sclerectomy

Il s'agit de comparer in vivo la sécurité et l'efficacité d'un implant en polyméthylméthacrylate (PMMA) avec un implant standard en collagène dans la sclérectomie profonde (SP) sur une durée de six mois. La population étudiée comprend vingt lapins, chaque lapin étant randomisé pour une SP avec implant en PMMA dans un oeil et implant de collagène dans l'autre oeil. Plusieurs éléments ont été pris en compte dans la comparaison : - la mesure de la pression intraoculaire - l'évolution de l'espace de drainage intrascléral et de la bulle de filtration sous-conjonctivale, suivie par ultrasonographic biomicroscopique (UBM) - la croissance de nouveaux vaisseaux de drainage sous-conjonctivaux, croissance quantifiée par angiographie du segment antérieur à la fluorescéine combinée au vert d'indocyanine - la facilité à l'écoulement de l'humeur aqueuse (C), mesurée à six mois par cannulation-perfusion de la chambre antérieur - la sclère au site de SP, histologiquement comparée à la sclère native opposée à 180°, également à six mois La pression intraoculaire moyenne préopératoire à une, quatre, douze et 24 semaines postopératoires est comparable dans les deux groupes (P>0.1). L'UBM montre une régression légèrement plus rapide (statistiquement non significative) de la bulle de filtration sous-conjonctivale et la persistance d'un espace de drainage intrascléral dans le groupe PMMA (P>0.05). De nouveaux vaisseaux de drainage sont observés à un mois de la chirurgie ; à six mois, ces vaisseaux sont plus nombreux dans le groupe PMMA, tant sur l'analyse angiographique que sur l'analyse histologique (P>0.05). La facilité moyenne à l'écoulement de l'humeur aqueuse est significativement plus élevées à six mois dans les deux groupes par rapport aux valeurs préopératoires (P>0.05), sans qu'il n'y ait de différence entre les deux implants (0.24 ± 0.06 μΙ/min/mmHg [PMMA] et 0.23 ± 0.07 μΙ/min/mmHg [implant en collagène]) (Ρ = 0.39). Cette étude a pu démontrer que la sclérectomie profonde avec implant en collagène ou en PMMA donne des résultats similaires en terme de diminution de l'IOP et d'augmentation de la facilité à l'écoulement de l'humeur aqueuse, sans différence sur le plan des réactions inflammatoires post-intervention.

Ultrastructure of the membrane attack complex of complement: detection of the tetramolecular C9-polymerizing complex C5b-8.

The ultrastructure of the membrane attack complex (MAC) of complement had been described as representing a hollow cylinder of defined dimensions that is composed of the proteins C5b, C6, C7, C8, and C9. After the characteristic cylindrical structure was identified as polymerized C9 [poly(C9)], the question arose as to the ultrastructural identity and topology of the C9-polymerizing complex C5b-8. An electron microscopic analysis of isolated MAC revealed an asymmetry of individual complexes with respect to their length. Whereas the length of one boundary (+/- SEM) was always 16 +/- 1 nm, the length of the other varied between 16 and 32 nm. In contrast, poly(C9), formed spontaneously from isolated C9, had a uniform tubule length (+/- SEM) of 16 +/- 1 nm. On examination of MAC-phospholipid vesicle complexes, an elongated structure was detected that was closely associated with the poly(C9) tubule and that extended 16-18 nm beyond the torus of the tubule and 28-30 nm above the membrane surface. The width of this structure varied depending on its two-dimensional projection in the electron microscope. By using biotinyl C5b-6 in the formation of the MAC and avidin-coated colloidal gold particles for the ultrastructural analysis, this heretofore unrecognized subunit of the MAC could be identified as the tetramolecular C5b-8 complex. Identification also was achieved by using anti-C5 Fab-coated colloidal gold particles. A similar elongated structure of 25 nm length (above the surface of the membrane) was observed on single C5b-8-vesicle complexes. It is concluded that the C5b-8 complex, which catalyzes poly(C9) formation, constitutes a structure of discrete morphology that remains as such identifiable in the fully assembled MAC, in which it is closely associated with the poly(C9) tubule.

Standardisation of 18F by a coincidence method using full solid angle detectors.

A solution of (18)F was standardised with a 4pibeta-4pigamma coincidence counting system in which the beta detector is a one-inch diameter cylindrical UPS89 plastic scintillator, positioned at the bottom of a well-type 5''x5'' NaI(Tl) gamma-ray detector. Almost full detection efficiency-which was varied downwards electronically-was achieved in the beta-channel. Aliquots of this (18)F solution were also measured using 4pigamma NaI(Tl) integral counting and Monte Carlo calculated efficiencies as well as the CIEMAT-NIST method. Secondary measurements of the same solution were also performed with an IG11 ionisation chamber whose equivalent activity is traceable to the Système International de Référence through the contribution IRA-METAS made to it in 2001; IRA's degree of equivalence was found to be close to the key comparison reference value (KCRV). The (18)F activity predicted by this coincidence system agrees closely with the ionisation chamber measurement and is compatible within one standard deviation of the other primary measurements. This work demonstrates that our new coincidence system can standardise short-lived radionuclides used in nuclear medicine.

Real-time recording of circadian liver gene expression in freely moving mice reveals the phase-setting behavior of hepatocyte clocks.

The mammalian circadian timing system consists of a master pacemaker in the suprachiasmatic nucleus (SCN) in the hypothalamus, which is thought to set the phase of slave oscillators in virtually all body cells. However, due to the lack of appropriate in vivo recording technologies, it has been difficult to study how the SCN synchronizes oscillators in peripheral tissues. Here we describe the real-time recording of bioluminescence emitted by hepatocytes expressing circadian luciferase reporter genes in freely moving mice. The technology employs a device dubbed RT-Biolumicorder, which consists of a cylindrical cage with reflecting conical walls that channel photons toward a photomultiplier tube. The monitoring of circadian liver gene expression revealed that hepatocyte oscillators of SCN-lesioned mice synchronized more rapidly to feeding cycles than hepatocyte clocks of intact mice. Hence, the SCN uses signaling pathways that counteract those of feeding rhythms when their phase is in conflict with its own phase.

Analysis of secondary growth in the Arabidopsis shoot reveals a positive role of jasmonate signalling in cambium formation.

After primary growth, most dicotyledonous plants undergo secondary growth. Secondary growth involves an increase in the diameter of shoots and roots through formation of secondary vascular tissue. A hallmark of secondary growth initiation in shoots of dicotyledonous plants is the initiation of meristematic activity between primary vascular bundles, i.e. in the interfascicular regions. This results in establishment of a cylindrical meristem, namely the vascular cambium. Surprisingly, despite its major implications for plant growth and the accumulation of biomass, the molecular regulation of secondary growth is only poorly understood. Here, we combine histological, molecular and genetic approaches to characterize interfascicular cambium initiation in the Arabidopsis thaliana inflorescence shoot. Using genome-wide transcriptional profiling, we show that stress-related and touch-inducible genes are up-regulated in stem regions where secondary growth takes place. Furthermore, we show that the products of COI1, MYC2, JAZ7 and the touch-inducible gene JAZ10, which are components of the JA signalling pathway, are cambium regulators. The positive effect of JA application on cambium activity confirmed a stimulatory role of JA in secondary growth, and suggests that JA signalling triggers cell divisions in this particular context.

cAMP-dependent chloride secretion mediates tubule enlargement and cyst formation by cultured mammalian collecting duct cells.

Polycystic kidney diseases result from disruption of the genetically defined program that controls the size and geometry of renal tubules. Cysts which frequently arise from the collecting duct (CD) result from cell proliferation and fluid secretion. From mCCD(cl1) cells, a differentiated mouse CD cell line, we isolated a clonal subpopulation (mCCD-N21) that retains morphogenetic capacity. When grown in three-dimensional gels, mCCD-N21 cells formed highly organized tubular structures consisting of a palisade of polarized epithelial cells surrounding a cylindrical lumen. Subsequent addition of cAMP-elevating agents (forskolin or cholera toxin) or of membrane-permeable cAMP analogs (CPT-cAMP) resulted in rapid and progressive dilatation of existing tubules, leading to the formation of cystlike structures. When grown on filters, mCCD-N21 cells exhibited a high transepithelial resistance as well as aldosterone- and/or vasopressin-induced amiloride-sensitive and -insensitive current. The latter was in part inhibited by Na(+)-K(+)-2Cl(-) cotransporter (bumetanide) and chloride channel (NPPB) inhibitors. Real-time PCR analysis confirmed the expression of NKCC1, the ubiquitous Na(+)-K(+)-2Cl(-) cotransporter and cystic fibrosis transmembrane regulator (CFTR) in mCCD-N21 cells. Tubule enlargement and cyst formation were prevented by inhibitors of Na(+)-K(+)-2Cl(-) cotransporters (bumetanide or ethacrynic acid) or CFTR (NPPB or CFTR inhibitor-172). These results further support the notion that cAMP signaling plays a key role in renal cyst formation, at least in part by promoting chloride-driven fluid secretion. This new in vitro model of tubule-to-cyst conversion affords a unique opportunity for investigating the molecular mechanisms that govern the architecture of epithelial tubes, as well as for dissecting the pathophysiological processes underlying cystic kidney diseases.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

A Miocene mud volcano and its plumbing system: A chaotic complex revisited (Monferrato, NW Italy)

Chaotic deposits are frequently reported in the geological literature and are commonly interpreted as olistostromes or tectonic melanges. A chaotic complex in the Cenozoic succession of Monferrato (NW Italy) consists of interbedded mud breccia and burrowed silty clays that are pierced by sheared mud breccias and embed carbonate-cemented blocks. These may be represented by microcrystalline limestones or strongly cemented matrix-supported breccias locally containing remains of chemosymbiotic organisms (lucinid bivalves). Moreover, cylindrical concretions, up to 15 cm in diameter and 1 m long, occur in the chaotic complex and crosscut bedding planes at high angles. The cement of all these lithified portions is mainly dolomite characterized by low delta(13)C values (from -10.3 to -23parts per thousand PDB) and delta(18)O values up to + 7parts per thousand PDB. The delta(13)C values testify to precipitation of carbonates induced by microbial oxidation of methane, whereas the markedly positive delta(18)C signature, ubiquitous in the cylindrical concretions, is the evidence for the presence and destabilization of gas hydrates. The studied section provides a well-exposed example of the geological record of the birth, life, and death of a mud volcano. Unsheared, soft mud breccias represent mud flows along the flanks of the volcano, whereas sheared mud breccias are the result of the injection of unconsolidated overpressured fine-grained sediments, both taking place during ``eruptive'' phases. They were followed by more quiet stages of hemipelagic sedimentation, burrowing, and CH4 seeping. The cylindrical concretions represent the first described ancient example of the chimneys observed in present-day mud-volcano settings. They are the remnants of a cold-seep plumbing network that crosscut the mud volcano edifice. The chimneys were the pathway for the expulsion toward the sea floor of gas- and sediment-charged fluids likely originated from destabilization of methane gas hydrates. The association of mud breccias and methane-derived carbonates may not be due to mass gravity flows but can be primary and, therefore, is a diagnostic criterion for recognizing chaotic deposits due to mud volcano activity in the geological record.

A pseudo-spectral method for simulating of poro-elastic seismic wave propagation in complex borehole environments

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in cylindrical coordinates. An important application of this method is the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the center, the borehole casing and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction, Fourier expansions in the other directions, and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method based on the method of characteristics is used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces. The viability and accuracy of the proposed method has been tested and verified in 2D polar coordinates through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently benchmarked solution for 2D Cartesian coordinates. The proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is handled adequately.

Geology of the NW Indian Himalaya

This review paper deals with the geology of the NW Indian Himalaya situated in the states of Jammu and Kashmir, Himachal Pradesh and Garhwal. The models and mechanisms discussed, concerning the tectonic and metamorphic history of the Himalayan range, are based on a new compilation of a geological map and cross sections, as well as on paleomagnetic, stratigraphic, petrologic, structural, metamorphic, thermobarometric and radiometric data. The protolith of the Himalayan range, the North Indian flexural passive margin of the Neo-Tethys ocean, consists of a Lower Proterozoic basement, intruded by 1.8-1.9 Ga bimodal magmatites, overlain by a horizontally stratified sequence of Upper Proterozoic to Paleocene sediments, intruded by 470-500 Ma old Ordovician mainly peraluminous s-type granites, Carboniferous tholeiitic to alkaline basalts and intruded and overlain by Permian tholeiitic continental flood basalts. No elements of the Archaen crystalline basement of the South Indian shield have been identified in the Himalayan range. Deformation of the Himalayan accretionary wedge resulted from the continental collision of India and Asia beginning some 65-55 Ma ago, after the NE-directed underthrusting of the Neo-Tethys oceanic crust below Asia and the formation of the Andean-type 103-50 (-41) Ma old Ladakh batholith to the north of the Indus Suture. Cylindrical in geometry, the Himalayan range consists, from NE to SW, from older to younger tectonic elements, of the following zones: 1) The 25 km wide Ladakh batholith and the Asian mantle wedge form the backstop of the growing Himalayan accretionary wedge. 2) The Indus Suture zone is composed of obducted slices of the oceanic crust, island arcs, like the Dras arc, overlain by Late Cretaceous fore arc basin sediments and the mainly Paleocene to Early Eocene and Miocene epi-sutural intra-continental Indus molasse. 3) The Late Paleocene to Eocene North Himalayan nappe stack, up to 40 km thick prior to erosion, consists of Upper Proterozoic to Paleocene rocks, with the eclogitic and coesite bearing Tso Morari gneiss nappe at its base. It includes a branch of the Central Himalayan detachment, the 22-18 Ma old Zanskar Shear zone that is intruded and dated by the 22 Ma Gumburanjun leucogranite; it reactivates the frontal thrusts of the SW-verging North Himalayan nappes. 4) The late Eocene-Miocene SW-directed High Himalayan or ``Crystalline'' nappe comprises Upper Proterozoic to Mesozoic sediments and Ordovician granites, identical to those of the North Himalayan nappes. The Main Central thrust at its base was created in a zone of Eocene to Early Oligocene anatexis by ductile detachment of the subducted Indian crust, below the pre-existing 25-35 km thick NE-directed Shikar Beh and SW-directed North Himalayan nappe stacks. 5) The late Miocene Lesser Himalayan thrust with the Main Boundary Thrust at its base consists of early Proterozoic to Cambrian rocks intruded by 1.8-1.9 Ga bimodal magmatites. The Subhimalaya is a thrust wedge of Himalayan fore deep basin sediments, composed of the Early Eocene marine Subathu marls and sandstones as well as the up to 8'000 m-thick Miocene to recent Ganga molasse, a coarsening upwards sequence of shales, sandstones and conglomerates. The active frontal thrust is covered by the sediments of the Indus-Ganga plains.

Analysis of bioavailable arsenic in rice with whole cell living bioreporter bacteria.

A multiwell plate bioassay was developed using genetically modified bacteria (bioreporter cells) to detect inorganic arsenic extracted from rice. The bacterial cells expressed luciferase upon exposure to arsenite, the activity of which was detected by measurement of cellular bioluminescence. The bioreporter cells detected arsenic in all rice varieties tested, with averages of 0.02-0.15 microg of arsenite equivalent per gram of dry weight and a method detection limit of 6 ng of arsenite per gram of dry rice. This amounted to between approximately 20 and 90% of the total As content reported by chemical methods for the same sample and suggested that a major proportion of arsenic in rice is in the inorganic form. Calibrations of the bioassay with pure inorganic and organic arsenic forms showed that the bacterial cells react to arsenite with highest affinity, followed by arsenate (with 25% response relative to an equivalent arsenite concentration) and trimethylarsine oxide (at 10% relative response). A method for biocompatible arsenic extraction was elaborated, which most optimally consisted of (i) grinding rice to powder, (ii) mixing with an aqueous solution containing pancreatic enzymes, (iii) mechanical shearing, (iv) extraction in mild acid conditions and moderate heat, and (v) centrifugation and pH neutralization. Detection of mainly inorganic arsenic by the bacterial cells may have important advantages for toxicity assessment of rice consumption and would form a good complement to total chemical arsenic determination.

The birth of the Rheic Ocean - Early Palaeozoic subsidence patterns and subsequent tectonic plate scenarios

New plate-tectonic reconstructions of the Gondwana margin suggest that the location of Gondwana-derived terranes should not only be guided by the models, but should also consider the possible detrital input from some Asian blocks (Hunia), supposed to have been located along the Cambrian Gondwana margin, and accreted in the Silurian to the North-Chinese block. Consequently, the Gondwana margin has to be subdivided into a more western domain, where the future Avalonian blocks will be separated from Gondwana by the opening Rheic Ocean, whereas in its eastern continuation, hosting the future basement areas of Central Europe, different periods of crustal extension should be distinguished. Instead of applying a rather cylindrical model, it is supposed that crustal extension follows a much more complex pattern, where local back-arcs or intra-continental rifts are involved. Guided by the age data of magmatic rocks and the pattern of subsidence curves, the following extensional events can be distinguished: During the early to middle Cambrian, a back-arc setting guided the evolution at the Gondwana margin. Contemporaneous intra-continental rift basins developed at other places related to a general post-PanAfrican extensional phase affecting Africa Upper Cambrian formation of oceanic crust is manifested in the Chamrousse area, and may have lateral cryptic relics preserved in other places. This is regarded as the oceanisation of some marginal basins in a context of back-arc rifting. These basins were closed in a mid-Ordovician tectonic phase, related to the subduction of buoyant material (mid-ocean ridge?) Since the Early Ordovician, a new phase of extension is observed, accompanied by a large-scale volcanic activity, erosion of the rift shoulders generated detritus (Armorican Quartzite) and the rift basins collected detrital zircons from a wide hinterland. This phase heralded the opening of Palaeotethys, but it failed due to the Silurian collision (Eo-Variscan phase) of an intra-oceanic arc with the Gondwana margin. During this time period, at the eastern wing of the Gondwana margin begins the drift of the future Hunia microcontinents, through the opening of an eastern prolongation of the already existing Rheic Ocean. The passive margin of the remaining Gondwana was composed of the Galatian superterranes, constituents of the future Variscan basement areas. Remaining under the influence of crustal extension, they will start their drift to Laurussia since the earliest Devonian during the opening of the Palaeotethys Ocean. (C) 2008 Elsevier B.V. All rights reserved.

Update on the non-prewhitening model observer in computed tomography for the assessment of the adaptive statistical and model-based iterative reconstruction algorithms.

The state of the art to describe image quality in medical imaging is to assess the performance of an observer conducting a task of clinical interest. This can be done by using a model observer leading to a figure of merit such as the signal-to-noise ratio (SNR). Using the non-prewhitening (NPW) model observer, we objectively characterised the evolution of its figure of merit in various acquisition conditions. The NPW model observer usually requires the use of the modulation transfer function (MTF) as well as noise power spectra. However, although the computation of the MTF poses no problem when dealing with the traditional filtered back-projection (FBP) algorithm, this is not the case when using iterative reconstruction (IR) algorithms, such as adaptive statistical iterative reconstruction (ASIR) or model-based iterative reconstruction (MBIR). Given that the target transfer function (TTF) had already shown it could accurately express the system resolution even with non-linear algorithms, we decided to tune the NPW model observer, replacing the standard MTF by the TTF. It was estimated using a custom-made phantom containing cylindrical inserts surrounded by water. The contrast differences between the inserts and water were plotted for each acquisition condition. Then, mathematical transformations were performed leading to the TTF. As expected, the first results showed a dependency of the image contrast and noise levels on the TTF for both ASIR and MBIR. Moreover, FBP also proved to be dependent of the contrast and noise when using the lung kernel. Those results were then introduced in the NPW model observer. We observed an enhancement of SNR every time we switched from FBP to ASIR to MBIR. IR algorithms greatly improve image quality, especially in low-dose conditions. Based on our results, the use of MBIR could lead to further dose reduction in several clinical applications.

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.