Foam pore size is a critical interface parameter of suction-based wound healing devices.

BACKGROUND: Suction-based wound healing devices with open-pore foam interfaces are widely used to treat complex tissue defects. The impact of changes in physicochemical parameters of the wound interfaces has not been investigated. METHODS: Full-thickness wounds in diabetic mice were treated with occlusive dressing or a suction device with a polyurethane foam interface varying in mean pore size diameter. Wound surface deformation on day 2 was measured on fixed tissues. Histologic cross-sections were analyzed for granulation tissue thickness (hematoxylin and eosin), myofibroblast density (α-smooth muscle actin), blood vessel density (platelet endothelial cell adhesion molecule-1), and cell proliferation (Ki67) on day 7. RESULTS: Polyurethane foam-induced wound surface deformation increased with polyurethane foam pore diameter: 15 percent (small pore size), 60 percent (medium pore size), and 150 percent (large pore size). The extent of wound strain correlated with granulation tissue thickness that increased 1.7-fold in small pore size foam-treated wounds, 2.5-fold in medium pore size foam-treated wounds, and 4.9-fold in large pore size foam-treated wounds (p < 0.05) compared with wounds treated with an occlusive dressing. All polyurethane foams increased the number of myofibroblasts over occlusive dressing, with maximal presence in large pore size foam-treated wounds compared with all other groups (p < 0.05). CONCLUSIONS: The pore size of the interface material of suction devices has a significant impact on the wound healing response. Larger pores increased wound surface strain, tissue growth, and transformation of contractile cells. Modification of the pore size is a powerful approach for meeting biological needs of specific wounds.

A single point mutation in the pore region of the epithelial Na+ channel changes ion selectivity by modifying molecular sieving.

The epithelial Na+ channel (ENaC) belongs to a new class of channel proteins called the ENaC/DEG superfamily involved in epithelial Na+ transport, mechanotransduction, and neurotransmission. The role of ENaC in Na+ homeostasis and in the control of blood pressure has been demonstrated recently by the identification of mutations in ENaC beta and gamma subunits causing hypertension. The function of ENaC in Na+ reabsorption depends critically on its ability to discriminate between Na+ and other ions like K+ or Ca2+. ENaC is virtually impermeant to K+ ions, and the molecular basis for its high ionic selectivity is largely unknown. We have identified a conserved Ser residue in the second transmembrane domain of the ENaC alpha subunit (alphaS589), which when mutated allows larger ions such as K+, Rb+, Cs+, and divalent cations to pass through the channel. The relative ion permeability of each of the alphaS589 mutants is related inversely to the ionic radius of the permeant ion, indicating that alphaS589 mutations increase the molecular cutoff of the channel by modifying the pore geometry at the selectivity filter. Proper geometry of the pore is required to tightly accommodate Na+ and Li+ ions and to exclude larger cations. We provide evidence that ENaC discriminates between cations mainly on the basis of their size and the energy of dehydration.

Simulation of interface dynamics across pore discontinuities

We study the dynamics of a water-oil meniscus moving from a smaller to a larger pore. The process is characterised by an abrupt change in the configuration, yielding a sudden energy release. A theoretic study for static conditions provides analytical solutions of the surface energy content of the system. Although the configuration after the sudden energy release is energetically more convenient, an energy barrier must be overcome before the process can happen spontaneously. The energy barrier depends on the system geometry and on the flow parameters. The analytical results are compared to numerical simulations that solve the full Navier-Stokes equation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. First, the numerical simulations of a quasi-static process are validated by comparison with the analytical solutions for a static meniscus, then numerical simulations with varying injection velocity are used to investigate dynamic effects on the configuration change. During the sudden energy jump the system exhibits an oscillatory behaviour. Extension to more complex geometries might elucidate the mechanisms leading to a dynamic capillary pressure and to bifurcations in final distributions of fluid phases in porous

Is it the grain size or the characteristic pore size that controls the induced polarization relaxation time of clean sands and sandstones?

There is a wide range of evidence to suggest that permeability can be constrained through of induced polarization measurements. For clean sands and sandstones, current mechanistic models of induced polarization predict a relationship between the low-frequency time constant inferred from induced polarization measurements and the grain diameter. A number of observations do, however, disagree with this and indicate that the observed relaxation behavior is rather governed by the so-called dynamic pore radius L. To test this hypothesis, we have developed a set of new scaling relationships, which allow the relaxation time to be computed from the pore size and the permeability to be computed from both the Cole-Cole time constant and the formation factor. Moreover, these new scaling relationships can be also used to predict the dependence of the Cole-Cole time constant as a function of the water saturation under unsaturated conditions. Comparative tests of the proposed new relationships with regard to various published experimental results for saturated clean sands and sandstones as well as for partially saturated clean sandstones, do indeed confirm that the dynamic pore radius L is a much more reliable indicator of the observed relaxation behavior than grain-size-based models.

Crystallization and preliminary X-ray analysis of monalysin, a novel β-pore-forming toxin from the entomopathogen Pseudomonas entomophila.

Monalysin was recently described as a novel pore-forming toxin (PFT) secreted by the Drosophila pathogen Pseudomonas entomophila. Recombinant monalysin is multimeric in solution, whereas PFTs are supposed to be monomeric until target membrane association. Monalysin crystals were obtained by the hanging-drop vapour-diffusion method using PEG 8000 as precipitant. Preliminary X-ray diffraction analysis revealed that monalysin crystals belonged to the monoclinic space group C2, with unit-cell parameters a = 162.4, b = 146.2, c = 144.4 Å, β = 122.8°, and diffracted to 2.85 Å resolution using synchrotron radiation. Patterson self-rotation analysis and Matthews coefficient calculation indicate that the asymmetric unit contains nine copies of monalysin. Heavy-atom derivative data were collected and a Ta6Br14 cluster derivative data set confirmed the presence of ninefold noncrystallographic symmetry.

Transition from hemifusion to pore opening is rate limiting for vacuole membrane fusion.

Fusion pore opening and expansion are considered the most energy-demanding steps in viral fusion. Whether this also applies to soluble N-ethyl-maleimide sensitive fusion protein attachment protein receptor (SNARE)- and Rab-dependent fusion events has been unknown. We have addressed the problem by characterizing the effects of lysophosphatidylcholine (LPC) and other late-stage inhibitors on lipid mixing and pore opening during vacuole fusion. LPC inhibits fusion by inducing positive curvature in the bilayer and changing its biophysical properties. The LPC block reversibly prevented formation of the hemifusion intermediate that allows lipid, but not content, mixing. Transition from hemifusion to pore opening was sensitive to guanosine-5'-(gamma-thio)triphosphate. It required the vacuolar adenosine triphosphatase V0 sector and coincided with its transformation. Pore opening was rate limiting for the reaction. As with viral fusion, opening the fusion pore may be the most energy-demanding step for intracellular, SNARE-dependent fusion reactions, suggesting that fundamental aspects of lipid mixing and pore opening are related for both systems.

Impact of controlled changes in grain size and pore space characteristics on the hydraulic conductivity and spectral induced polarization response of "proxies" of saturated alluvial sediments

Understanding the influence of pore space characteristics on the hydraulic conductivity and spectral induced polarization (SIP) response is critical for establishing relationships between the electrical and hydrological properties of surficial unconsolidated sedimentary deposits, which host the bulk of the world's readily accessible groundwater resources. Here, we present the results of laboratory SIP measurements on industrial-grade, saturated quartz samples with granulometric characteristics ranging from fine sand to fine gravel, which can be regarded as proxies for widespread alluvial deposits. We altered the pore space characteristics by changing (i) the grain size spectra, (ii) the degree of compaction, and (iii) the level of sorting. We then examined how these changes affect the SIP response, the hydraulic conductivity, and the specific surface area of the considered samples. In general, the results indicate a clear connection between the SIP response and the granulometric as well as pore space characteristics. In particular, we observe a systematic correlation between the hydraulic conductivity and the relaxation time of the Cole-Cole model describing the observed SIP effect for the entire range of considered grain sizes. The results do, however, also indicate that the detailed nature of these relations depends strongly on variations in the pore space characteristics, such as, for example, the degree of compaction. The results of this study underline the complexity of the origin of the SIP signal as well as the difficulty to relate it to a single structural factor of a studied sample, and hence raise some fundamental questions with regard to the practical use of SIP measurements as site- and/or sample-independent predictors of the hydraulic conductivity.

Geochemical evolution of active porphyry copper tailings impoundments : from alkaline deposition towards acidification

J. Smuda: Geochemical evolution of active porphyry copper tailings impoundments Thesis abstract Mine waste is the largest volume of materials handled in the world. The oxidation of sulfidic mine waste may result in the release of acid mine drainage (AMD) rich in heavy metals and arsenic to the environment, one of the major problems the mining industry is facing today. To control and reduce this environmental impact, it is crucial to identify the main geochemical and hydrological processes influencing contaminant liberation, transport, and retention. This thesis presents the results of a geochemical, mineralogical and stable isotope study (δ2H, δ18O, δ34S) from two active porphyry copper tailings impoundments in Mediterranean (Carén tailings impoundment, El Teniente mine, Central Chile) and hyper-arid climate (Talabre tailings impoundment, Chuquicamata, Northern Chile) from the deposition in alkaline environment (pH 10.5) towards acidification after several years of exposure. The major hydrological results were the identification of vertical contaminant and water transport in the uppermost, not water-saturated zone, triggered by capillary rise due to evaporation, and infiltration downwards due to new tailings deposition, and of horizontal transport in the groundwater zone. At the surface of the sedimented tailings, evaporation of pore water led to the precipitation of Na-Ca-Mg sulfates (e.g., gypsum, tenorite), in hyper-arid climate also halite. At the Carén tailings impoundment, renewed deposition in a 4-week interval inhibited a pH decrease below neutral values and the formation of an efflorescent salt crust. At the Talabre tailings impoundment, deposition breaks of several years resulted in the formation of acidic oxidation zones in the timeframe of less than 4 years. This process enabled the transport of liberated Cu, Zn, and Fe via capillary rise to the surface, where these metals precipitated as heavy-metal sulfates (e.g., devilline, krohnkite) and chlorides (eriochalcite, atacamite). Renewed depositing may dissolve efflorescent salts and transport liberated elements towards the groundwater zone. This zone was found to be highly dynamic due to infiltration and mixing with water from different sources, like groundwater, catchment water, and infiltration from superficial waters. There, Cu was found to be partially mobile due to complexation with Cl (in Cl-rich groundwater, Talabre) and dissolved organic matter (in zones with infiltration of catchment water rich in dissolved organic matter, Carén). A laboratory study on the isotopic fractionation of sulfur and oxygen of sulfate in different minerals groups (water-soluble sulfates, low- and high-crystalline Fe(III) oxyhydroxides) contributed to the use of stable isotopes as tracer of geochemical and transport processes for environmental studies. The results highlight that a detailed geochemical, stable isotope and mineralogical study permits the identification of contamination processes and pathways already during the deposition of mine tailings. This knowledge allows the early planning of adequate actions to reduce and control the environmental impact during tailings deposition and after the closing of the impoundment. J. Smuda: Geochemical evolution of active porphyry copper tailings impoundments Résumé de these Les déchets miniers constituent les plus grands volumes de matériel gérés dans le monde. L'oxydation des déchets miniers sulfuriques peut conduire à la libération de drainages miniers acides (DMA) riches en métaux et arsenic dans l'environnement, ce qui est l'un des principaux problèmes de l'industrie minière aujourd'hui. Pour contrôler et réduire ces impacts sur l'environnement, il est crucial d'identifier les principaux processus géochimiques et hydrologiques influençant la libération, le transport et la rétention des contaminants. Cette thèse présente les résultats d'une étude géochimique, minéralogique et des isotopes stables (δ2H, δ18O, δ34S) sur des déchets miniers de 2 sites de dépôt actifs en climat méditerranéen (Dépôt de déchets de Carén, mine de El Teniente, Centre du Chili) et en climat hyper-aride (Dépôt de déchets de Talabre, mine de Chuquicamata, Nord du Chili). L'objectif était d'étudier l'évolution des déchets de la déposition en milieu alcalin (pH = 10.5) vers l'acidification après plusieurs années d'exposition. Le principal résultat hydrologique a été l'identification de 2 types de transport : un transport vertical de l'eau et des contaminants dans la zone non saturée en surface, induit par la montée capillaire due à l'évaporation et par l'infiltration subséquente de la déposition de sédiments frais ; et un transport horizontal dans la zone des eaux souterraines. À la surface des déchets, l'évaporation de l'eau interstitielle conduit à la précipitation de sulfates de Na-Ca-Mg (ex. gypse, ténorite) et halite en climat hyper-aride. Dans le site de Carén, une nouvelle déposition de déchets frais à 4 semaines intervalle a empêché la baise du pH en deçà des valeurs neutres et la formation d'une croûte de sels efflorescentes en surface. Dans le site de Talabre, les fentes de dessiccation des dépôts ont entraîné la formation d'une zone d'oxydation à pH acide en moins de 4 ans. Ce processus a permis la libération et le transport par capillarité de Cu, Zn, Fe vers la surface, où ces éléments précipitent sous forme de sulfates de métaux lourds (ex., dévilline, krohnkite) de chlorures (ex. ériochalcite, atacamite). Une nouvelle déposition de sédiments frais pourrait dissoudre ces sels et les transporter vers la zone des eaux souterraines. Cette dernière zone était très dynamique en raison du mélange d'eaux provenant de différentes sources, comme les eaux souterraines, l'eau de captage et l'infiltration des eaux superficielles. Egalement dans cette zone, le cuivre était partiellement mobile à cause de la formation de complexe avec le chlore (dans les zone riche en Cl, Talabre) et avec la matière organique dissoute (dans les zones où s'infiltre l'eau de captage riche en matière organique, Carén). Une étude en laboratoire sur le fractionnement des isotopes stables de sulfure et d'oxygène des sulfates dans différents groupes de minéraux (sulfates hydrosolubles, sulfures de oxy-hydroxyde de Fe(III) faiblement ou fortement cristallins) a permis d'apporter une contribution à leur utilisation comme traceurs dans l'étude des processus géochimiques et de transport lors d'études environnementales. Les résultats montrent qu'une étude détaillée de la géochimie, des isotopes stables et de la minéralogie permet d'identifier les processus et les voies de contamination déjà pendant la période de dépôt des déchets miniers. Cette connaissance permet de planifier, dès le début de l'exploitation, des mesures adéquates pour réduire et contrôler l'impact sur l'environnement pendant la période de dépôts de déchets miniers et après la fermeture du site.

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).

Spectral induced polarization: Laboratory measurements on sandy quartz samples of varied pore characteristics

Understanding the influence of pore space characteristics on the hydraulic conductivity and spectral induced polarization (SIP) response is critical for establishing relationships between the electrical and hydrological properties of surficial sedimentary deposits. Here, we present the results of laboratory SIP measurements on saturated quartz samples with granulometric characteristics ranging from fine sand to fine gravel. We alter the pore characteristics using three principal methods: (i) variation of the grain sizes, (ii) changing the degree of compaction, and (iii) changing the level of sorting. We then examine how these changes affect both the SIP response and the hydraulic conductivity. In general, the results indicate a clear connection between the applied changes in pore characteristics and the SIP response. In particular, we observe a systematic correlation between the hydraulic conductivity and the relaxation time of the Cole-Cole model describing the observed SIP effect for the whole range of considered grain sizes.

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.