Detection and Characterization of Preferential Flow Paths in the Downstream Area of a Hazardous Landfill

We have used surface-based electrical resistivity tomography to detect and characterize preferential hydraulic pathways in the immediate downstream area of an abandoned, hazardous landfill. The landfill occupies the void left by a former gravel pit and its base is close to the groundwater table and lacking an engineered barrier. As such, this site is remarkably typical of many small- to medium-sized waste deposits throughout the densely populated and heavily industrialized foreland on both sides of the Alpine arc. Outflows of pollutants lastingly contaminated local drinking water supplies and necessitated a partial remediation in the form of a synthetic cover barrier, which is meant to prevent meteoric water from percolating through the waste before reaching the groundwater table. Any future additional isolation of the landfill in the form of lateral barriers thus requires adequate knowledge of potential preferential hydraulic pathways for outflowing contaminants. Our results, inferred from a suite of tomographically inverted surfaced-based electrical resistivity profiles oriented roughly perpendicular to the local hydraulic gradient, indicate that potential contaminant outflows would predominantly occur along an unexploited lateral extension of the original gravel deposit. This finds its expression as a distinct and laterally continuous high-resistivity anomaly in the resistivity tomograms. This interpretation is ground-truthed through a litholog from a nearby well. Since the probed glacio-fluvial deposits are largely devoid of mineralogical clay, the geometry of hydraulic and electrical pathways across the pore space of a given lithological unit can be assumed to be identical, which allows for an order-of-magnitude estimation of the overall permeability structure. These estimates indicate that the permeability of the imaged extension of the gravel body is at least two to three orders-of-magnitude higher than that of its finer-grained embedding matrix. This corroborates the preeminent role of the high-resistivity anomaly as a potential preferential flow path.

Contribution of (sub)domains of Staphylococcus aureus fibronectin-binding protein to the proinflammatory and procoagulant response of human vascular endothelial cells.

The Staphylococcus aureus fibronectin (Fn) -binding protein A (FnBPA) is involved in bacterium-endothelium interactions which is one of the crucial events leading to infective endocarditis (IE). We previously showed that the sole expression of S. aureus FnBPA was sufficient to confer to non-invasive Lactococcus lactis bacteria the capacity to invade human endothelial cells (ECs) and to launch the typical endothelial proinflammatory and procoagulant responses that characterize IE. In the present study we further questioned whether these bacterium-EC interactions could be reproduced by single or combined FnBPA sub-domains (A, B, C or D) using a large library of truncated FnBPA constructs expressed in L. lactis. Significant invasion of cultured ECs was found for L. lactis expressing the FnBPA subdomains CD (aa 604-877) or A4(+16) (aa 432-559). Moreover, this correlates with the capacity of these fragments to elicit in vitro a marked increase in EC surface expression of both ICAM-1 and VCAM-1 and secretion of the CXCL8 chemokine and finally to induce a tissue factor-dependent endothelial coagulation response. We thus conclude that (sub)domains of the staphylococcal FnBPA molecule that express Fn-binding modules, alone or in combination, are sufficient to evoke an endothelial proinflammatory as well as a procoagulant response and thus account for IE severity.

Morphological changes in convection-diffusion-limited deposition

The effect of hydrodynamic flow upon diffusion-limited deposition on a line is investigated using a Monte Carlo model. The growth process is governed by the convection and diffusion field. The convective diffusion field is simulated by the biased-random walker resulting from a superimposed drift that represents the convective flow. The development of distinct morphologies is found with varying direction and strength of drift. By introducing a horizontal drift parallel to the deposition plate, the diffusion-limited deposit changes into a single needle inclined to the plate. The width of the needle decreases with increasing strength of drift. The angle between the needle and the plate is about 45° at high flow rate. In the presence of an inclined drift to the plate, the convection-diffusion-limited deposit leads to the formation of a characteristic columnar morphology. In the limiting case where the convection dominates, the deposition process is equivalent to ballistic deposition onto an inclined surface.

The effect of lattice defects on Hele-Shaw flow over an etched lattice

We examine the patterns formed by injecting nitrogen gas into the center of a horizontal, radial Hele-Shaw cell filled with paraffin oil. We use smooth plates and etched plates with lattices having different amounts of defects (010 %). In all cases, a quantitative measure of the pattern ramification shows a regular trend with injection rate and cell gap, such that the dimensionless perimeter scales with the dimensionless time. By adding defects to the lattice, we observe increased branching in the pattern morphologies. However, even in this case, the scaling behavior persists. Only the prefactor of the scaling function shows a dependence on the defect density. For different lattice defect densities, we examine the nature of the different morphology phases.

Patterns of gravity induced aggregate migration during casting of fluid concretes

In this paper, aggregate migration patterns during fluid concrete castings are studied through experiments, dimensionless approach and numerical modeling. The experimental results obtained on two beams show that gravity induced migration is primarily affecting the coarsest aggregates resulting in a decrease of coarse aggregates volume fraction with the horizontal distance from the pouring point and in a puzzling vertical multi-layer structure. The origin of this multi layer structure is discussed and analyzed with the help of numerical simulations of free surface flow. Our results suggest that it finds its origin in the non Newtonian nature of fresh concrete and that increasing casting rate shall decrease the magnitude of gravity induced particle migration. (C) 2012 Elsevier Ltd. All rights reserved.

Characterization of surface functional groups present on laboratory-generated and ambient aerosol particles by means of heterogeneous titration reactions

A Knudsen flow reactor has been used to quantify surface functional groups on aerosols collected in the field. This technique is based on a heterogeneous titration reaction between a probe gas and a specific functional group on the particle surface. In the first part of this work, the reactivity of different probe gases on laboratory-generated aerosols (limonene SOA, Pb(NO3)2, Cd(NO3)2) and diesel reference soot (SRM 2975) has been studied. Five probe gases have been selected for the quantitative determination of important functional groups: N(CH3)3 (for the titration of acidic sites), NH2OH (for carbonyl functions), CF3COOH and HCl (for basic sites of different strength), and O3 (for oxidizable groups). The second part describes a field campaign that has been undertaken in several bus depots in Switzerland, where ambient fine and ultrafine particles were collected on suitable filters and quantitatively investigated using the Knudsen flow reactor. Results point to important differences in the surface reactivity of ambient particles, depending on the sampling site and season. The particle surface appears to be multi-functional, with the simultaneous presence of antagonistic functional groups which do not undergo internal chemical reactions, such as acid-base neutralization. Results also indicate that the surface of ambient particles was characterized by a high density of carbonyl functions (reactivity towards NH2OH probe in the range 0.26-6 formal molecular monolayers) and a low density of acidic sites (reactivity towards N(CH3)3 probe in the range 0.01-0.20 formal molecular monolayer). Kinetic parameters point to fast redox reactions (uptake coefficient ?0>10-3 for O3 probe) and slow acid-base reactions (?0<10-4 for N(CH3)3 probe) on the particle surface. [Authors]

Foot skin blood flow following infrainguinal revascularization for critical lower limb ischemia.

INTRODUCTION: The aim of this study was to assess the blood flow in the feet before and after lower limb revascularization using laser Doppler imaging (LDI). METHODS: Ten patients with critical lower limb ischemia were prospectively enrolled from June to October 2004. All patients underwent successful unilateral surgical interventions including above-knee bypass, distal bypass and endarterectomy. Skin blood flow (SBF) over the plantar surface of both forefeet and heels was measured by LDI 24h before and 10 days after revascularization, expressed in perfusion units (PU), and reported as mean+/-SD. RESULTS: Measurements in the forefoot and heel were similar. Before revascularization mean SBF was significantly lower in the ischemic foot (130+/-71 PU) compared to the contralateral foot (212+/-68 PU), p<0.05. After revascularization a significant increase of the SBF in the forefoot (from 135+/-67 to 202+/-86 PU, p=0.001) and hindfoot (from 148+/-58 to 203+/-83, p=0.001) was observed on the treatment side. However, a large decrease of the SBF was seen in forefoot and hindfoot on the untreated side (from 250+/-123 PU to 176+/-83 and from 208+/-116 to 133+/-40, p=0.001, respectively). CONCLUSION: This study confirms the benefits of revascularization in patients with nonhealing foot lesions due to critical limb ischemia. A significant increase of the SBF was observed on the treatment side. However, an unexpected decrease was observed on the untreated side.

Development of Elispot and CFSE flow cytometric assays for detecting beryllium sensitivity

Background: Beryllium sensitization (BeS) is caused by exposure to beryllium in the workplace and may progress to chronic beryllium disease (CBD). This granulomatous lung disorder mimicks sarcoidosis clinically, but is characterized by beryllium specific CD4+ T-cells immune response. BeS is classically detected by beryllium lymphocyte proliferation test (BeLPT), but this assay requires radioactivity and is not very sensitive. In the context of a study aiming to evaluate if CBD patients are misdiagnosed as sarcoidosis patients in Switzerland, we developed EliSpot and CFSE beryllium flow cytometric test. Methods: 23 patients considered as having sarcoidosis (n = 21), CBD (n = 1) and possible CBD (n = 1) were enrolled. Elispot was performed using plate covered with gamma-IFN mAb. Cells were added to wells and incubated overnight at 37 °C with medium (neg ctrl), SEB (pos ctrl) or BeSO4 at 1, 10 and 100 microM. Anti-IFN-gamma biotinylated mAb were added and spots were visualized using streptavidinhorseradish peroxidase and AEC substrate reagent. Results were reported as spot forming unit (SFU). For Beryllium specific CFSE flow cytometry analysis, CFSE labelled cells were cultured in the presence of SEB and 1, 10 or 100 microM BeSO4. Unstimulated CFSE labeled cells were defined as controls. The cells were incubated for 6 days at 37 °C and 5% CO2. Surface labelling of T-lymphocytes and vivid as control of cells viability was performed at the time of harvest. Results: Using EliSpot technology, we were able to detect a BeS in 1/23 enrolled patients with a mean of 780 SFU (cut off value at 50 SFU). This positive result was confirmed using different concentration of BeSO4. Among the 23 patients tested, 22 showed negative results with EliSpot. Using CFSE flow cytometry, 1/7 tested patients showed a positive result with a beryllium specific CD4+ count around 30% versus 45% for SEB stimulation as positif control and 0.6 % for negative control. This patient was the one with a positive EliSpot assay. Conclusions: The preliminary data demonstrated the feasibility of Elispot and CFSE flow cytometry to detect BeS. The patient with a beryllium specific positive EliSpot and CFSE flow cytometry result had been exposed to beryllium at her workplace 20 years ago and is still regularly controlled for her pulmonary status. A positive BeLPT had already been described in 2001 in France for this patient. Further validation of these techniques are in progress.

LiveCount Assay: concomitant measurement of cytolytic activity and phenotypic characterisation of CD8(+) T-cells by flow cytometry.

Tumour immunologists strive to develop efficient tumour vaccination and adoptive transfer therapies that enlarge the pool of tumour-specific and -reactive effector T-cells in vivo. To assess the efficiency of the various strategies, ex vivo assays are needed for the longitudinal monitoring of the patient's specific immune responses providing both quantitative and qualitative data. In particular, since tumour cell cytolysis is the end goal of tumour immunotherapy, routine immune monitoring protocols need to include a read-out for the cytolytic efficiency of Ag-specific cells. We propose to combine current immune monitoring techniques in a highly sensitive and reproducible multi-parametric flow cytometry based cytotoxicity assay that has been optimised to require low numbers of Ag-specific T-cells. The possibility of re-analysing those T-cells that have undergone lytic activity is illustrated by the concomitant detection of CD107a upregulation on the surface of degranulated T-cells. To date, the LiveCount Assay provides the only possibility of assessing the ex vivo cytolytic activity of low-frequency Ag-specific cytotoxic T-lymphocytes from patient material.

Is it appropriate to index glomerular filtration rate for body surface area

Purpose: Obesity is an established independent risk factor for chronic kidney disease. Thus, measurement of glomerular filtration rate (GFR) is important in this population. Traditionally, GFR has been indexed for body surface area (BSA), but this indexation may not be appropriate in obese individuals. Therefore, the objective of the study was to compare absolute GFR with GFR indexed for BSA and with GFR indexed for height. Methods and materials: The study was conducted in 66 families from the Seychelles islands that included several members with hypertension. GFR and effective renal plasma flow (ERPF) were measured using inulin and PAH clearances, respectively. Antihypertensive treatment, if used, was withheld 2 weeks before conducting the clearances. Participants with diabetes mellitus were excluded from the analysis. BSA was calculated using the Dubois formula. We assessed trend across BMI categories using a non parametric test. Results: Participants included 174 women and 127 men. The prevalence of hypertension was 61%, of which 68% were treated. The table shows that absolute GFR, GFR indexed for height, ERPF, filtration fraction were significantly higher across BMI categories. When GFR was indexed for BSA, the association between GFR and BMI categories was lost. Conclusion: Indexing GFR for BSA in overweight and obese individuals leads to a substantial underestimation of GFR. Filtration fraction, which does not depend on BSA, is higher in obese individuals, which suggests glomerular hyperfiltration. Indexing GFR for BSA therefore would mask the underlying glomerular hyperfiltration. As the number of nephrons does not increase with weight gain, absolute GFR represents a better marker of single nephron GFR and is more appropriate.

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.

Direct numerical simulations of interface dynamics to link capillary pressure and total surface energy

The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work we perform direct numerical simulations of drainage by solving Navier-Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters.

Fate of the herbicide glyphosate and its metabolite AMPA in soils and their transfer to surface waters: A multi-scale approach in the Lavaux vineyards, western Switzerland

The use of herbicides in agriculture may lead to environmental problems, such as surface water pollution, with a potential risk for aquatic organisms. The herbicide glyphosate is the most used active ingredient in the world and in Switzerland. In the Lavaux vineyards it is nearly the only molecule applied. This work aimed at studying its fate in soils and its transfer to surface waters, using a multi-scale approach: from molecular (10-9 m) and microscopic scales (10-6 m), to macroscopic (m) and landscape ones (103 m). First of all, an analytical method was developed for the trace level quantification of this widely used herbicide and its main by-product, aminomethylphosphonic acid (AMPA). Due to their polar nature, their derivatization with 9-fluorenylmethyl chloroformate (FMOC-Cl) was done prior to their concentration and purification by solid phase extraction. They were then analyzed by ultra performance liquid chromatography coupled with tandem mass spectrometry (UPLC-MS/MS). The method was tested in different aqueous matrices with spiking tests and validated for the matrix effect correction in relevant environmental samples. Calibration curves established between 10 and 1000ng/l showed r2 values above 0.989, mean recoveries varied between 86 and 133% and limits of detection and quantification of the method were as low as 5 and 10ng/l respectively. At the parcel scale, two parcels of the Lavaux vineyard area, located near the Lutrive River at 6km to the east of Lausanne, were monitored to assess to which extent glyphosate and AMPA were retained in the soil or exported to surface waters. They were equipped at their bottom with porous ceramic cups and runoff collectors, which allowed retrieving water samples for the growing seasons 2010 and 2011. Results revealed that the mobility of glyphosate and AMPA in the unsaturated zone was likely driven by the precipitation regime and the soil characteristics, such as slope, porosity structure and layer permeability discrepancy. Elevated glyphosate and AMPA concentrations were measured at 60 and 80 cm depth at parcel bottoms, suggesting their infiltration in the upper parts of the parcels and the presence of preferential flow in the studied parcels. Indeed, the succession of rainy days induced the gradual saturation of the soil porosity, leading to rapid infiltration through macropores, as well as surface runoff formation. Furthermore, the presence of more impervious weathered marls at 100 cm depth induced throughflows, the importance of which for the lateral transport of the herbicide molecules was determined by the slope steepness. Important rainfall events (>10 mm/day) were clearly exporting molecules from the soil top layer, as indicated by important concentrations in runoff samples. A mass balance showed that total loss (10-20%) mainly occurred through surface runoff (96%) and, to a minor extent, by throughflows in soils (4%), with subsequent exfiltration to surface waters. Observations made in the Lutrive River revealed interesting details of glyphosate and AMPA dynamics in urbanized landscapes, such as the Lavaux vineyards. Indeed, besides their physical and chemical properties, herbicide dynamics at the catchment level strongly depend on application rates, precipitation regime, land use and also on the presence of drains or constructed channels. Elevated concentrations, up to 4970 ng/l, observed just after the application, confirmed the diffuse export of these compounds from the vineyard area by surface runoff during main rain events. From April to September 2011, a total load of 7.1 kg was calculated, with 85% coming from vineyards and minor urban sources and 15% from arable crops. Small vineyard surfaces could generate high concentrations of herbicides and contribute considerably to the total load calculated at the outlet, due to their steep slopes (~10%). The extrapolated total amount transferred yearly from the Lavaux vineyards to the Lake of Geneva was of 190kg. At the molecular scale, the possible involvement of dissolved organic matter (DOM) in glyphosate and copper transport was studied using UV/Vis fluorescence spectroscopy. Combined with parallel factor (PARAFAC) analysis, this technique allowed characterizing DOM of soil and surface water samples from the studied vineyard area. Glyphosate concentrations were linked to the fulvic-like spectroscopic signature of DOM in soil water samples, as well as to copper, suggesting the formation of ternary complexes. In surface water samples, its concentrations were also correlated to copper ones, but not in a significant way to the fulvic-like signature. Quenching experiments with standards confirmed field tendencies in the laboratory, with a stronger decrease in fluorescence intensity for fulvic-like fluorophore than for more aromatic ones. Lastly, based on maximum concentrations measured in the river, an environmental risk for these compounds was assessed, using laboratory tests and ecotoxicity data from the literature. In our case and with the methodology applied, the risk towards aquatic species was found negligible (RF<1).

Assessment of diesel exhaust particulate exposure and surface characteristics in association with levels of oxidative stress biomarkers

Exposure to PM10 and PM2.5 (particulate matter with aerodynamic diameter smaller than 10 μm and 2.5 μm, respectively) is associated with a range of adverse health effects, including cancer, pulmonary and cardiovascular diseases. Surface characteristics (chemical reactivity, surface area) are considered of prime importance to understand the mechanisms which lead to harmful effects. A hypothetical mechanism to explain these adverse effects is the ability of components (organics, metal ions) adsorbed on these particles to generate Reactive Oxygen Species (ROS), and thereby to cause oxidative stress in biological systems (Donaldson et al., 2003). ROS can attack almost any cellular structure, like DNA or cellular membrane, leading to the formation of a wide variety of degradation products which can be used as a biomarker of oxidative stress. The aim of the present research project is to test whether there is a correlation between the exposure to Diesel Exhaust Particulate (DEP) and the oxidative stress status. For that purpose, a survey has been conducted in real occupational situations where workers were exposed to DEP (bus depots). Different exposure variables have been considered: - particulate number, size distribution and surface area (SMPS); - particulate mass - PM2.5 and PM4 (gravimetry); - elemental and organic carbon (coulometry); - total adsorbed heavy metals - iron, copper, manganese (atomic adsorption); - surface functional groups present on aerosols (Knudsen flow reactor). Several biomarkers of oxidative stress (8-hydroxy-2'-deoxyguanosine and several aldehydes) have been determined either in urine or serum of volunteers. Results obtained during the sampling campaign in several bus depots indicated that the occupational exposure to particulates in these places was rather low (40-50 μg/m3 for PM4). Bimodal size distributions were generally observed (5 μm and <1 μm). Surface characteristics of PM4 varied strongly, depending on the bus depot. They were usually characterized by high carbonyl and low acidic sites content. Among the different biomarkers which have been analyzed within the framework of this study, mean urinary levels of 8-hydroxy-2'-deoxyguanosine increased significantly (p<0.05) during two consecutive days of exposure for non-smoker workers. On the other hand, no statistically significant differences were observed for serum levels of hexanal, nonanal and 4- hydroxy-nonenal (p>0.05). Biomarkers levels will be compared to exposure variables to gain a better understanding of the relation between the particulate characteristics and the formation of ROS by-products. This project is financed by the Swiss State Secretariat for Education and Research. It is conducted within the framework of the COST Action 633 "Particulate Matter - Properties Related to Health Effects".

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.