MORM--a Petri net based model for assessing OH&S risks in industrial processes: modeling qualitative aspects.

Because of the increase in workplace automation and the diversification of industrial processes, workplaces have become more and more complex. The classical approaches used to address workplace hazard concerns, such as checklists or sequence models, are, therefore, of limited use in such complex systems. Moreover, because of the multifaceted nature of workplaces, the use of single-oriented methods, such as AEA (man oriented), FMEA (system oriented), or HAZOP (process oriented), is not satisfactory. The use of a dynamic modeling approach in order to allow multiple-oriented analyses may constitute an alternative to overcome this limitation. The qualitative modeling aspects of the MORM (man-machine occupational risk modeling) model are discussed in this article. The model, realized on an object-oriented Petri net tool (CO-OPN), has been developed to simulate and analyze industrial processes in an OH&S perspective. The industrial process is modeled as a set of interconnected subnets (state spaces), which describe its constitutive machines. Process-related factors are introduced, in an explicit way, through machine interconnections and flow properties. While man-machine interactions are modeled as triggering events for the state spaces of the machines, the CREAM cognitive behavior model is used in order to establish the relevant triggering events. In the CO-OPN formalism, the model is expressed as a set of interconnected CO-OPN objects defined over data types expressing the measure attached to the flow of entities transiting through the machines. Constraints on the measures assigned to these entities are used to determine the state changes in each machine. Interconnecting machines implies the composition of such flow and consequently the interconnection of the measure constraints. This is reflected by the construction of constraint enrichment hierarchies, which can be used for simulation and analysis optimization in a clear mathematical framework. The use of Petri nets to perform multiple-oriented analysis opens perspectives in the field of industrial risk management. It may significantly reduce the duration of the assessment process. But, most of all, it opens perspectives in the field of risk comparisons and integrated risk management. Moreover, because of the generic nature of the model and tool used, the same concepts and patterns may be used to model a wide range of systems and application fields.

Modelling hydrodynamics in the Rio Parana, Argentina : an evaluation and inter-comparison of reduced-complexity and physics based models applied to a large sand-bed river

Depth-averaged velocities and unit discharges within a 30 km reach of one of the world's largest rivers, the Rio Parana, Argentina, were simulated using three hydrodynamic models with different process representations: a reduced complexity (RC) model that neglects most of the physics governing fluid flow, a two-dimensional model based on the shallow water equations, and a three-dimensional model based on the Reynolds-averaged Navier-Stokes equations. Row characteristics simulated using all three models were compared with data obtained by acoustic Doppler current profiler surveys at four cross sections within the study reach. This analysis demonstrates that, surprisingly, the performance of the RC model is generally equal to, and in some instances better than, that of the physics based models in terms of the statistical agreement between simulated and measured flow properties. In addition, in contrast to previous applications of RC models, the present study demonstrates that the RC model can successfully predict measured flow velocities. The strong performance of the RC model reflects, in part, the simplicity of the depth-averaged mean flow patterns within the study reach and the dominant role of channel-scale topographic features in controlling the flow dynamics. Moreover, the very low water surface slopes that typify large sand-bed rivers enable flow depths to be estimated reliably in the RC model using a simple fixed-lid planar water surface approximation. This approach overcomes a major problem encountered in the application of RC models in environments characterised by shallow flows and steep bed gradients. The RC model is four orders of magnitude faster than the physics based models when performing steady-state hydrodynamic calculations. However, the iterative nature of the RC model calculations implies a reduction in computational efficiency relative to some other RC models. A further implication of this is that, if used to simulate channel morphodynamics, the present RC model may offer only a marginal advantage in terms of computational efficiency over approaches based on the shallow water equations. These observations illustrate the trade off between model realism and efficiency that is a key consideration in RC modelling. Moreover, this outcome highlights a need to rethink the use of RC morphodynamic models in fluvial geomorphology and to move away from existing grid-based approaches, such as the popular cellular automata (CA) models, that remain essentially reductionist in nature. In the case of the world's largest sand-bed rivers, this might be achieved by implementing the RC model outlined here as one element within a hierarchical modelling framework that would enable computationally efficient simulation of the morphodynamics of large rivers over millennial time scales. (C) 2012 Elsevier B.V. All rights reserved.

Properties of debris flow deposits and source materials compared: implications for debris flow characterization

In the Alps, debris flow deposits generally contain < 5% clay-size particles, and the role of the surface-charged < 2 mu m particles is often neglected, although these particles may have a significant impact on the rheological properties of the interstitial fluid. The objective of this study was to compare debris flow deposits and parent materials from two neighbouring catchments of the Swiss Alps, with special emphasis on the colloidal constituents. The catchments are small in area (4 km(2)), 2.5 km long, similar in morphology, but different in geology. The average slopes are 35-40%. The catchments were monitored for debris flow events and mapped for surface aspect and erosion activity. Debris flow deposits and parent materials were sampled, the clay and silt fractions extracted and the bulk density, < 2 mm fraction bulk density, particle size distribution, chemical composition, cation exchange capacity (CEC) and mineralogy analysed. The results show that the deposits are similar to the parent screes in terms of chemical composition, but differ in terms of: (i) particle size distribution; and (ii) mineralogy, reactivity and density of the < 2 mm fraction. In this fraction, compared with the parent materials the deposits show dense materials enriched in coarse monocrystalline particles, of which the smallest and more reactive particles were leached. The results suggest that deposit samples should not be considered as representative of source or flow materials, particularly with respect to their physical properties.

Ambient vertical flow in long-screen wells: a case study in the Fontainebleau Sands Aquifer (France)

A tritium (H-3) profile was constructed in a long-screened well (LSW) of the Fontainebleau Sands Aquifer (France), and the data were combined with temperature logs to gain insight into the potential effects of the ambient vertical flow (AVF) of water through the well on the natural aquifer stratification. AVF is commonly taken into account in wells located in fracture aquifers or intercepting two different aquifers with distinct hydraulic heads. However, due to the vertical hydraulic gradient of the flow lines intercepted by wells, AVF of groundwater is a common process within any type of aquifer. The detection of 3H in the deeper parts of the studied well ( approximate depth 50m), where H-3-free groundwater is expected, indicates that shallow young water is being transported downwards through the well itself. The temperature logs show a nearly zero gradient with depth, far below the mean geothermal gradient in sedimentary basins. The results show that the age distribution of groundwater samples might be biased in relation to the age distribution in the surroundings of the well. The use of environmental tracers to investigate aquifer properties, particularly in LSWs, is then limited by the effects of the AVF of water that naturally occurs through the well.

Anastomotic longitudinal stress due to modification of arterial longitudinal properties after anastomosis.

BACKGROUND: In our hands, in vivo segmental vessel length changes up to 5% because of blood pressure: increasing in arterial pressure is associated to decrease in segmental vessel length. METHODS AND MATERIAL: Using two piezoelectric crystals sutured on vessel wall and a high fidelity pressure probe, we recorded artery length variations as function of blood pressure, before and after an end-to-end anastomosis on four pigs carotid arteries. RESULTS: Mean arterial pressure before anastomosis = 73 mmHg (+/- 12); mean arterial pressure after anastomosis = 91 mmHg (+/- 14); mean crystals displacement before anastomosis during systole = -0.21 mm; mean crystals displacement after anastomosis during systole = +0.24 mm; mean distance between crystals before anastomosis = 12.3 mm (+/- 0.8) and after anastomosis = 11.2 mm (+/- 0.5). CONCLUSIONS: In the acute phase following an end-to-end anastomosis, an increase in blood pressure causes increasing in vessel length, with an exponential correlation. The anastomosis is constantly subjected to a longitudinal traction whose magnitude depends on blood pressure.

Magnetic and in vitro heating properties of implants formed in situ from injectable formulations and containing superparamagnetic iron oxide nanoparticles (SPIONs) embedded in silica microparticles for magnetically induced local hyperthermia

The biological and therapeutic responses to hyperthermia, when it is envisaged as an anti-tumor treatment modality, are complex and variable. Heat delivery plays a critical role and is counteracted by more or less efficient body cooling, which is largely mediated by blood flow. In the case of magnetically mediated modality, the delivery of the magnetic particles, most often superparamagnetic iron oxide nanoparticles (SPIONs), is also critically involved. We focus here on the magnetic characterization of two injectable formulations able to gel in situ and entrap silica microparticles embedding SPIONs. These formulations have previously shown suitable syringeability and intratumoral distribution in vivo. The first formulation is based on alginate, and the second on a poly(ethylene-co-vinyl alcohol) (EVAL). Here we investigated the magnetic properties and heating capacities in an alternating magnetic field (141 kHz, 12 mT) for implants with increasing concentrations of magnetic microparticles. We found that the magnetic properties of the magnetic microparticles were preserved using the formulation and in the wet implant at 37 degrees C, as in vivo. Using two orthogonal methods, a common SLP (20 Wg(-1)) was found after weighting by magnetic microparticle fraction, suggesting that both formulations are able to properly carry the magnetic microparticles in situ while preserving their magnetic properties and heating capacities. (C) 2010 Elsevier B.V. All rights reserved.

Comparative Study of Chemoembolization, Loadable Beads: In vitro Drug Release and Physical Properties of DC Bead and Hepasphere, Loaded with Doxorubicin and Irinotecan

Purpose To characterize in vitro the loadability, physical properties, and release of irinotecan and doxorubicin from two commercially available embolization microspheres. Materials and Methods DC Bead (500-700 μm) and Hepasphere (400-600 μm) microspheres were loaded with either doxorubicin or irinotecan solutions. Drug amount was quantified with spectrophotometry, bead elasticity was measured under compression, and bead size and loading homogeneity were assessed with microscopy image analysis. Drug release was measured over 1-week periods in saline by using a pharmacopeia flow-through method. Results Almost complete drug loading was obtained for both microsphere types and drugs. Doxorubicin-loaded DC Beads maintained their spherical shape throughout the release. In contrast, Hepaspheres showed less homogeneous doxorubicin loading and, after release, some fractured microspheres. Incomplete doxorubicin release was observed in saline over 1 week (27% ± 2 for DC beads and 18% ± 7 for Hepaspheres; P = .013). About 75% of this amount was released within 2.2 hours for both beads. For irinotecan, complete release was obtained for both types of beads, in a sustained manner over 2-3 hours for DC Beads, and in a significantly faster manner as a 7-minute burst for Hepaspheres. Conclusions The two drug-eluting microspheres could be efficiently loaded with both drugs. Incomplete doxorubicin release was attributed to strong drug-bead ionic interactions. Weaker interactions were observed with irinotecan, which led to faster drug release.

Formaldehyde-fixation of platelets for flow cytometric measurement of phosphatidylserine exposure is feasible.

Strong platelet activation results in a redistribution of negatively charged phospholipids from the cytosolic to the outer leaflet of the cellular membrane. Annexin V has a high affinity to negatively charged phospholipids and can be used to identify procoagulant platelets. Formaldehyde fixation can cause factitious Annexin V binding. Our aim was to evaluate a method for fixing platelets avoiding additional Annexin V binding. We induced expression of negatively charged phospholipids on the surface of a fraction of platelets by combined activation with convulxin and thrombin in the presence of Annexin V-fluorescein isothiocyanate and calcium. Aliquots of resting and activated platelets were fixed with a low concentration, calcium-free formaldehyde solution. Both native platelets and fixed platelets were analyzed by flow cytometry immediately and after a 24-h storage at 4°C. We observed that the percentage of Annexin V positive resting platelets ranged from 1.5 to 9.3% for the native samples and from 0.4 to 12.8% for the fixed samples (P=0.706, paired t-test). The amount of Annexin V positive convulxin/thrombin activated platelets varied from 12.9 to 35.4% without fixation and from 15.3 to 36.3% after formalin fixation (P=0.450). After a 24-h storage at 4°C, Annexin V positive platelets significantly increased both in the resting and in the convulxin/thrombin activated samples of native platelets (both P<0.001), while results for formalin fixed platelets did not differ from baseline values (P=0.318 for resting fixed platelets; P=0.673 for activated fixed platelets). We conclude that platelet fixation with a low concentration, calcium-free formaldehyde solution does not alter the proportion of Annexin V positive platelets. This method can be used to investigate properties of procoagulant platelets by multicolor flow-cytometric analysis requiring fixation steps.

Numerical analysis of wave-induced fluid flow effects on seismic data: Application to monitoring of CO2 storage at the Sleipner Field

In this work we analyze how patchy distributions of CO2 and brine within sand reservoirs may lead to significant attenuation and velocity dispersion effects, which in turn may have a profound impact on surface seismic data. The ultimate goal of this paper is to contribute to the understanding of these processes within the framework of the seismic monitoring of CO2 sequestration, a key strategy to mitigate global warming. We first carry out a Monte Carlo analysis to study the statistical behavior of attenuation and velocity dispersion of compressional waves traveling through rocks with properties similar to those at the Utsira Sand, Sleipner field, containing quasi-fractal patchy distributions of CO2 and brine. These results show that the mean patch size and CO2 saturation play key roles in the observed wave-induced fluid flow effects. The latter can be remarkably important when CO2 concentrations are low and mean patch sizes are relatively large. To analyze these effects on the corresponding surface seismic data, we perform numerical simulations of wave propagation considering reservoir models and CO2 accumulation patterns similar to the CO2 injection site in the Sleipner field. These numerical experiments suggest that wave-induced fluid flow effects may produce changes in the reservoir's seismic response, modifying significantly the main seismic attributes usually employed in the characterization of these environments. Consequently, the determination of the nature of the fluid distributions as well as the proper modeling of the seismic data constitute important aspects that should not be ignored in the seismic monitoring of CO2 sequestration problems.

Measurements of seismic attenuation and transient fluid pressure in partially saturated Berea sandstone: Evidence of fluid flow on the mesoscopic scale

A novel laboratory technique is proposed to investigate wave-induced fluid flow on the mesoscopic scale as a mechanism for seismic attenuation in partially saturated rocks. This technique combines measurements of seismic attenuation in the frequency range from 1 to 100?Hz with measurements of transient fluid pressure as a response of a step stress applied on top of the sample. We used a Berea sandstone sample partially saturated with water. The laboratory results suggest that wave-induced fluid flow on the mesoscopic scale is dominant in partially saturated samples. A 3-D numerical model representing the sample was used to verify the experimental results. Biot's equations of consolidation were solved with the finite-element method. Wave-induced fluid flow on the mesoscopic scale was the only attenuation mechanism accounted for in the numerical solution. The numerically calculated transient fluid pressure reproduced the laboratory data. Moreover, the numerically calculated attenuation, superposed to the frequency-independent matrix anelasticity, reproduced the attenuation measured in the laboratory in the partially saturated sample. This experimental?numerical fit demonstrates that wave-induced fluid flow on the mesoscopic scale and matrix anelasticity are the dominant mechanisms for seismic attenuation in partially saturated Berea sandstone.

Estimating vadose zone hydraulic properties using ground penetrating radar: The impact of prior information

A number of geophysical methods, such as ground-penetrating radar (GPR), have the potential to provide valuable information on hydrological properties in the unsaturated zone. In particular, the stochastic inversion of such data within a coupled geophysical-hydrological framework may allow for the effective estimation of vadose zone hydraulic parameters and their corresponding uncertainties. A critical issue in stochastic inversion is choosing prior parameter probability distributions from which potential model configurations are drawn and tested against observed data. A well chosen prior should reflect as honestly as possible the initial state of knowledge regarding the parameters and be neither overly specific nor too conservative. In a Bayesian context, combining the prior with available data yields a posterior state of knowledge about the parameters, which can then be used statistically for predictions and risk assessment. Here we investigate the influence of prior information regarding the van Genuchten-Mualem (VGM) parameters, which describe vadose zone hydraulic properties, on the stochastic inversion of crosshole GPR data collected under steady state, natural-loading conditions. We do this using a Bayesian Markov chain Monte Carlo (MCMC) inversion approach, considering first noninformative uniform prior distributions and then more informative priors derived from soil property databases. For the informative priors, we further explore the effect of including information regarding parameter correlation. Analysis of both synthetic and field data indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when we combine these data with a realistic, informative prior.

In vitro characterization of immune-related properties of human fetal bone cells for potential tissue engineering applications.

We describe herein some immunological properties of human fetal bone cells recently tested for bone tissue-engineering applications. Adult mesenchymal stem cells (MSCs) and osteoblasts were included in the study for comparison. Surface markers involved in bone metabolism and immune recognition were analyzed using flow cytometry before and after differentiation or treatment with cytokines. Immunomodulatory properties were studied on activated peripheral blood mononuclear cells (PBMCs). The immuno-profile of fetal bone cells was further investigated at the gene expression level. Fetal bone cells and adult MSCs were positive for Stro-1, alkaline phosphatase, CD10, CD44, CD54, and beta2-microglobulin, but human leukocyte antigen (HLA)-I and CD80 were less present than on adult osteoblasts. All cells were negative for HLA-II. Treatment with recombinant human interferon gamma increased the presence of HLA-I in adult cells much more than in fetal cells. In the presence of activated PBMCs, fetal cells had antiproliferative effects, although with patterns not always comparable with those of adult MSCs and osteoblasts. Because of the immunological profile, and with their more-differentiated phenotype than of stem cells, fetal bone cells present an interesting potential for allogeneic cell source in tissue-engineering applications.

Effects of a nonpressor dose of neuropeptide Y on cardiac output, regional blood flow distribution and plasma renin, vasopressin and catecholamine levels.

Neuropeptide Y (NPY) is a peptide with vasoconstrictor properties known to be present in the central nervous system as well as in sympathetic nerve endings and the adrenal medulla. The purposes of this study were to investigate in normotensive conscious rats the effects of nonpressor doses of NPY on cardiac output and regional blood flow distribution (using radiolabeled microspheres) as well as on plasma renin activity, plasma catecholamine and vasopressin levels. NPY (0.1 microgram/min) infused i.v. for 30 min modified neither blood pressure nor heart rate. Cardiac index was at comparable levels in NPY- as in vehicle-treated rats (17.7 +/- 1.6, n = 8, vs. 21.3 +/- 0.9 ml/min/100 g, n = 8, mean +/- S.E.M.). There was no significant difference in regional blood flow distribution between the two groups of rats, except for the large intestine (0.42 +/- 0.06 vs. 0.71 +/- 0.1 ml/min/g in NPY- and vehicle-treated rats, respectively, P less than .05). Basal plasma renin activity and catecholamine levels were not modified by NPY whereas plasma vasopressin levels were lower (P less than .05) in rats given NPY (0.76 +/- 0.3 pg/ml, n = 8) than in those having received the vehicle (2.2 +/- 0.4 pg/ml).(ABSTRACT TRUNCATED AT 250 WORDS)

Statistical properties of population differentiation estimators under stepwise mutation in a finite island model.

Microsatellite loci mutate at an extremely high rate and are generally thought to evolve through a stepwise mutation model. Several differentiation statistics taking into account the particular mutation scheme of the microsatellite have been proposed. The most commonly used is R(ST) which is independent of the mutation rate under a generalized stepwise mutation model. F(ST) and R(ST) are commonly reported in the literature, but often differ widely. Here we compare their statistical performances using individual-based simulations of a finite island model. The simulations were run under different levels of gene flow, mutation rates, population number and sizes. In addition to the per locus statistical properties, we compare two ways of combining R(ST) over loci. Our simulations show that even under a strict stepwise mutation model, no statistic is best overall. All estimators suffer to different extents from large bias and variance. While R(ST) better reflects population differentiation in populations characterized by very low gene-exchange, F(ST) gives better estimates in cases of high levels of gene flow. The number of loci sampled (12, 24, or 96) has only a minor effect on the relative performance of the estimators under study. For all estimators there is a striking effect of the number of samples, with the differentiation estimates showing very odd distributions for two samples.

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.