Les émissions liées à la combustion de la biomasse et des carburants fossiles constituent des sources significatives de particules dans l'atmosphère. Les particules de suie émises par les turboréacteurs aéronautiques dans la basse stratosphère et la haute troposphère, zones sensibles de l'atmosphère, ont un impact sur le forçage radiatif global de par la formation de<br />traînées de condensation qui peuvent évoluer en nuage de type cirrus. Par ailleurs, ces particules de taille submicrométrique, dont l'impact sur la qualité de l'air est avéré, sont également source de questionnement quant à leur(s) effet(s) sur la santé publique.<br />La présente étude a pour objet la détermination expérimentale des caractéristiques physiques et chimiques des particules de suie émises par des turboréacteurs d'avions commerciaux. Les prélèvements sont effectués sur un banc d'essai de la société SNECMA<br />(groupe SAFRAN) par un dispositif mis au point au cours de ce travail et localisé derrière des turboréacteurs CFM56-5C et CFM56-5B qui équipent respectivement des avions long-courriers et moyen-courriers.<br />Diverses techniques expérimentales telles que la microscopie électronique en transmission (MET), la microscopie électronique à balayage (MEB), la spectrométrie dispersive d'énergie de rayons X (XREDS), et la spectroscopie infrarouge à transformée de Fourier (IRTF) sont mises en oeuvre pour étudier la morphologie, la structure, la distribution en taille, la composition chimique élémentaire des particules primaires de suie, et la nature des groupes fonctionnels présents à leur surface. La dimension fractale ( Df ) des agrégats de suie et leur diamètre de giration ( d g ) sont également établis.<br />Les résultats expérimentaux montrent que les suies émises par les deux turboréacteurs<br />étudiés présentent :<br />- des valeurs similaires de la dimension fractale des agrégats<br />- des distributions en taille log-normales pour les particules primaires et les agrégats<br />- des compositions chimiques élémentaires et des groupes fonctionnels très similaires<br />- des diamètres de giration des agrégats qui diffèrent<br />Par ailleurs, nous montrons que la phase de préparation des échantillons, préalable aux observations en MEB, a une influence importante sur la valeur déterminée de la moyenne géométrique de la distribution en taille des particules primaires. Par conséquent, nous préconisons de réaliser l'étude des caractéristiques physiques des particules primaires de suie aéronautique par MET plutôt que par MEB.

Si l’Inde est un pays majoritairement rural, sa population urbaine est plus grande que la population totale des États-Unis. Ainsi, à la fin du 20ème siècle, le phénomène urbain indien est déjà bien développé, d’autant plus qu’il possède une histoire plurimillénaire. La présence de deux perturbations ayant modifié les logiques urbaines au cours du 20ème siècle, nous a amené à nous questionner sur la résilience du système urbain indien. Afin d’en présenter les principaux mécanismes, nous avons retracé l’évolution du système urbain indien au cours du 20ème siècle. Pour ce faire, nous avons utilisé une approche géographique par la population, qu’une approche qualitative et davantage historique est venue compléter. Nous nous sommes ainsi intéressés à l’évolution de la hiérarchie urbaine et à celle des rythmes de croissance des villes. Après avoir mis en lumière l’existence de mécanismes de rééquilibrage de la trame urbaine au cours du 20ème siècle, nous nous sommes penchés sur la stabilité organisationnelle et structurelle du système urbain. Cela nous a permis de mettre en évidence l’existence de classes dynamiques au sein de la hiérarchie urbaine. Leur évolution met en exergue le dynamisme du système urbain durant la seconde moitié du siècle dernier et nous présente les étapes de sa résilience. / If India is a predominantly rural country, its urban population is larger than the total population of the United States. Thus, at the end of the 20th century, the Indian urban phenomenon is already well developed, especially as it has a plurimillenary history. The presence of two disturbances have changed the urban logic in the 20th century and led us to question the resilience of the Indian urban system.To present the main mechanisms, we traced the evolution of the Indian urban system across the 20th century. To do this, we used a geographical approach by the population and a more historical one in complement. We emphasized the evolution of the urban hierarchy and that of the cities growth rates.After showing the existence of balancing mechanisms of the urban framework in the 20th century, we focused on the organizational and structural stability of the urban system. This allowed us to highlight the existence of dynamic classes in the urban hierarchy. Their evolution highlights the dynamism of the urban system in the second half of the last century and shows us the stages of its resilience.

Parmi les différentes technologies pour produire un effluent conforme aux normes de rejet en vigueur, les Procédés par Boues Activées sont largement utilisés en tant que traitement biologique du fait de leur faible coût et de leur facilité d’implémentation. Ces procédés réalisent l’abattement des polluants en deux étapes : (1) la biodégradation par des microorganismes dans un bassin d’aération et (2) la séparation de la biomasse par décantation. Les paramètres opératoires peuvent agir sur l’efficacité du procédé, notamment en modifiant les propriétés des boues . L’objectif de ce travail est d’apporter un nouvel éclairage sur cette relation triangulaire en mettant en perspectives les propriétés de taille et de forme des flocs de boues et les performances du procédé de traitement. La première étape pour parvenir à cet objectif a consisté à développer une méthodologie de caractérisation des flocs de boues par analyse d’images à l’aide du microscope automatisé Morphologi G3 de Malvern Instruments. Le protocole de mesure proposé repose sur le choix de la cellule d’analyse, du facteur de dilution, du grossissement, du seuil de niveau de gris mis en œuvre pour la binarisation des images et du nombre total de particules analysées. Des boues provenant de plusieurs stations d’épuration ont été étudiées afin d’évaluer la robustesse de la méthodologie. Des différences en termes de propriétés de taille et forme - Diamètre équivalent (Deq), Circularité (C), Convexité (Cx), Solidité (S) et Rapport d’Aspect (RA) - entre les différentes boues sont clairement apparues. Le protocole de caractérisation de la morphologie étant établi et validé, deux campagnes expérimentales ont été mises en œuvre en utilisant un système à boues activées à l’échelle pilote alimenté par un effluent synthétique. Les descripteurs morphologiques ont été mesurés sur des échantillons d’environ 100 000 flocs permettant d’atteindre une statistique fiable aussi bien sur la taille que sur la forme des agrégats. La performance du procédé a été évaluée en suivant des indicateurs physicochimiques classiquement associés au traitement des eaux usées : Matières en Suspension (MES), Indice de Boue (IB), taux d’abattement de la Demande Chimique d’Oxygène (DCO) et de l’Azote (N). Cette stratégie a permis d’identifier des relations entre les performances d’épuration, la morphologie des agrégats et les paramètres opératoires. Pendant la première campagne, trois valeurs différentes de l’âge de boue ont été imposées : 15, 20 et 30 jours. Lors d’un changement d’âge de boue, la taille et la circularité des flocs biologiques augmentent temporairement avant de se stabiliser à nouveau à leurs valeurs initiales. Ainsi, en régime permanent, les distributions en taille et en circularité semblent indépendantes de l’âge de boue. Lors de la deuxième campagne, trois paramètres opératoires ont été modifiés : la charge organique, le taux de recirculation et le rapport DCO:N. Une charge organique élevée a provoqué une forte augmentation de la taille des flocs et une diminution de leur solidité. Ce résultat peut être expliqué par la croissance de bactéries filamenteuses qui impacte fortement la morphologie des flocs et la performance du procédé. L’augmentation du taux de recirculation et la diminution du rapport DCO:N ont, en outre, induit une légère réduction de la taille des flocs et une évolution peu significative des valeurs de circularité et convexité. L’analyse des distributions en 3 Dimensions couplant les données de taille et de forme s’est révélée un outil prometteur, permettant de proposer des mécanismes responsables de l’évolution des propriétés morphologiques des boues. Enfin, une autre approche a été développée sur la base de deux dimensions fractales reliant l’aire et le périmètre des flocs (Dpf) ou l’aire et la longueur maximale (D2). De brusques variations de Dpf et D2 ont été corrélées avec les phénomènes liés au développement de bactéries filamenteuses.

Plusieurs études sur la toxicité des particules dans l'air ambiant montrent que la phase polaire, qui contient entre autres des Hydrocarbures Aromatiques Polycycliques oxygénés et nitrés (OHAP et NHAP), est plus toxique que la fraction aromatique contenant les HAP classiques. Dans ce contexte, l'objectif de ce travail était de caractériser chimiquement les NHAP et OHAP dans l'air ambiant de différents milieux (trafic, urbain, péri-urbain, rural et altitude). Les prélèvements ont été réalisés dans le cadre du programme de recherche POllution des Vallées Alpines (POVA) au cours de l'hiver 2002-2003 et de l'été 2003 et dans la région de Marseille l'été 2004. Les deux phases, gazeuse et particulaire, ont été étudiées et une caractérisation des particules par tranche granulométrique a également été effectuée. Les résultats obtenus sur les OHAP et les NHAP ont permis d'alimenter la base de données très limitée sur ce type de composés et ce pour des atmosphères variées. Le développement d'une méthode d'analyse simultanée par GC/NICI-MS des NHAP et OHAP sur des matrices complexes, facile à mettre en œuvre et utilisable en routine, nous a permis de montrer que les niveaux de concentrations des OHAP dans l'air ambiant sont du même ordre de grandeur que ceux des HAP alors que les concentrations des NHAP sont de 1 à 2 ordres de grandeur inférieures. L'estimation du risque cancérigène à partir de facteurs d'équivalence toxique montre que le risque attribué aux NHAP peut atteindre 20% du risque total. L'étude de rapports spécifiques de NHAP a montré que sur les sites proches des sources de pollution, les NHAP ont une origine primaire forte alors que leur formation par réaction en phase gazeuse est évidente en été sur les sites ruraux éloignés des sources directes de pollution et qu'elle a lieu principalement de jour, initiée par OH. Il apparaît que les émissions des véhicules essence sont une source non négligeable de NHAP et de OHAP dans l'atmosphère. Le OHAP 9,10-anthraquinone présente une double origine (primaire diesel et formée par photochimie) sans pour autant pouvoir distinguer le lieu de la formation (sur le support du prélèvement ou dans l'atmosphère). La partition entre les deux phases des HAP, OHAP et NHAP, est extrêmement dépendante de la masse moléculaire du composé considéré et de la température. L'étude de cette partition semble pouvoir apporter des informations quant à l'origine primaire ou secondaire des dérivés des HAP. Enfin, les HAP, OHAP et NHAP sont essentiellement (>80%) associés à la fraction fine de l'aérosol (Dp<1,3 µm) et ce quelque soit la saison considérée (hiver et été). Les caractéristiques chimiques des différentes classes de composés (polarité…) peuvent néanmoins jouer un rôle quant à leur distribution en taille.

Les objectifs principaux des travaux de recherche réalisés sont de caractériser, quantifier et corréler les changements induits lors du recyclage sur les fibres, la pâte et le papier. Pour ce faire, des techniques de caractérisation spécifique ont été utilisées telles que la chromatographie inverse d'exclusion stérique (ISEC), l'analyse mécanique dynamique (DMA), la microscopie électronique à balayage environnementale (ESEM), la microscopie à force atomique (AFM) et la tomographie à rayons X. Le racornissement des fibres à elle seule ne peut pas expliquer totalement la perte de résistance des fibres. Lors du recyclage la largeur des fibres, l'épaisseur des parois fibreuses, la courbure, le nombre de coudes et d'irrégularité diminuent. Les points faibles à l'intérieur de la paroi augmentent dans les premiers et les derniers cycles de recyclage. Les fibres deviennent plus dures et cassantes à l'état sec. Les forces capillaires et la friction de surface augmentent à l'état humide. La surface de liaison entre fibres dans le réseau fibreux diminue initialement à cause de la perte de la flexibilité des fibres à l'état humide et des éléments fins alors que l'augmentation qui suit peut être reliée au collapse du lumen. Comme la solidité des fibres de diminue pas, la baisse des caractéristiques mécaniques du papier pourrait être attribuée à la dégradation de la qualité des liaisons et plus particulièrement à une délamination partielle de la couche P/S1.

L’objectif de ce travail était d’étudier l’évolution de la morphologie des agrégats soumis à des conditions hydrodynamiques turbulentes. Pour cela, des expériences de floculation de particules de latex ont été réalisées dans un réacteur de Taylor-Couette en mettant en œuvre un séquençage des conditions hydrodynamiques, consistant en six étapes effectuées alternativement à faible et fort taux de cisaillement. Pour l’ensemble des expériences, un suivi en ligne de la taille des flocs a été réalisé par granulométrie laser. Par ailleurs, des paramètres caractéristiques de la forme des agrégats, ont été déterminés par analyse d’images sur des échantillons de suspension prélevés à la fin des différentes étapes du séquençage hydrodynamique. Dans le cas où le mécanisme de floculation est la neutralisation de charges (par ajout de NaCl ou de Al2(SO4)3), les résultats ont permis de confirmer la forte influence de l’hydrodynamique sur la taille des agrégats mais ont aussi révélé l’existence d’un lien entre la forme des agrégats et leur histoire. Contrairement à la taille, la morphologie des flocs, caractérisée par leur circularité ou leur convexité est peu sensible à l’hydrodynamique. L’étude des distributions surfaciques de ces paramètres de forme ainsi que la représentation des résultats expérimentaux sous forme bidimensionnelle (reliant taille et forme) ont en effet montré qu’après un certain nombre de cycles rupture-refloculation, la circularité comme la convexité évoluent vers une distribution unique, indépendante de l’hydrodynamique. Dans le cas où le mécanisme de floculation est le pontage (le polymère utilisé étant le PolyDADMAC), pour des taux de cisaillement faible, les agrégats sont plus grands, mais plus circulaires et plus résistants à la rupture que ceux créés par neutralisation de charge. Une analyse complémentaire est ensuite proposée permettant d’analyser l’effet des conditions hydrodynamiques sur la dimension fractale des flocs déduite des mesures par diffraction laser ou de l’analyse d’images. A des fins de comparaison, une nouvelle expression de la dimension fractale issue de l’analyse morphologique Df* a été introduite. Enfin, la modélisation de la floculation par résolution de l’équation de bilan de population basée sur la méthode de la quadrature des moments (QMOM) a été abordée ; les paramètres variables du modèle étant les valeurs expérimentales du taux de cisaillement moyen et la dimension fractale. Les résultats préliminaires montrent qu’il est possible de reproduire correctement l’évolution des six premiers moments de la distribution en taille des agrégats pour une phase d’agrégation du séquençage hydrodynamique. / The aim of this work was to investigate the evolution of the morphology of aggregates under turbulent hydrodynamic conditions. For this purpose, latex particles flocculation experiments were conducted in a Taylor-Couette reactor by carrying out a sequencing of the hydrodynamic conditions, consisting of six steps performed alternately at low and high shear rate. For all the experiments, an on-line monitoring of the size of the flocs was carried out by laser diffraction. Moreover, characteristic parameters of the shape of the aggregates were determined by image analysis on suspension samples taken at the end of the various stages of the hydrodynamic sequencing. When the flocculation mechanism is the charge neutralization (by addition of NaCl or Al2(SO4)3), the results have confirmed the strong influence of hydrodynamics on aggregate size but also revealed the existence of a relationship between the shape of the aggregates and their history. Contrary to size, the morphology of the flocs, characterized either by their circularity or convexity, is not sensitive to hydrodynamics. The study of the surface distributions of these shape parameters as well as the representation of the experimental results in two-dimensional form (linking size and shape) have shown that, after a number of breakage-reflocculation cycles, both circularity and convexity evolve towards a single distribution which is independent of hydrodynamics. When the bridging flocculation is the main mechanism (the polymer used being PolyDADMAC), at low shear rates, the aggregates are larger, but more circular and more resistant to breakage than those issued from charge neutralization. A complementary analysis is then proposed to analyze the effect of hydrodynamics on the fractal dimension measured either by laser diffraction or by image analysis. In order to compare both results, a new expression of the fractal dimension derived from the morphological analysis Df * was introduced. Finally, the modeling of the flocculation by solving the population balance equation using the quadrature method of moments (QMOM) is presented; the only variable parameters being the experimental values of the average shear rate and the fractal dimension. The preliminary results show that it is possible to reproduce correctly the evolution of the first six moments of the aggregate size distribution for an aggregation phase of the hydrodynamic sequencing.

Cette thèse a pour objectif de fournir aux utilisateurs de SMPS (Scanning Mobility Particle Sizer) une méthodologie pour calculer les incertitudes associées à l'estimation de la granulométrie en nombre des aérosols. Le résultat de mesure est le comptage des particules de l'aérosol en fonction du temps. Estimer la granulométrie en nombre de l'aérosol à partir des mesures CNC revient à considérer un problème inverse sous incertitudes.Une revue des modèles existants est présentée dans le premier chapitre. Le modèle physique retenu fait consensus dans le domaine d'application.Dans le deuxième chapitre, un critère pour l'estimation de la granulométrie en nombre qui couple les techniques de régularisation et de la décomposition sur une base d'ondelettes est décrit.La nouveauté des travaux présentés réside dans l'estimation de granulométries présentant à la fois des variations lentes et des variations rapides. L'approche multi-échelle que nous proposons pour la définition du nouveau critère de régularisation permet d'ajuster les poids de la régularisation sur chaque échelle du signal. La méthode est alors comparée avec la régularisation classique. Les résultats montrent que les estimations proposées par la nouvelle méthode sont meilleures (au sens du MSE) que les estimations classiques.Le dernier chapitre de cette thèse traite de la propagation de l'incertitude à travers le modèle d'inversiondes données. C'est une première dans le domaine d'application car aucune incertitude n'est associée actuellement au résultat de mesure. Contrairement à l'approche classique qui utilise un modèle fixe pour l'inversion en faisant porter l'incertitude sur les entrées, nous proposons d'utiliser un modèle d'inversion aléatoire (tirage Monte-Carlo) afin d'intégrer les erreurs de modèle. Une estimation moyenne de la granulométrie en nombre de l'aérosol et une incertitude associée sous forme d'une région de confiance à 95 % est finalement présentée sur quelques mesures réelles.

Cette étude porte sur le broyage en voie humide (nanobroyage) de produits organiques à usage pharmaceutique dans un broyeur à billes agité. L'objectif des travaux réalisés était d'analyser la faisabilité de l'opération de nanobroyage et d'améliorer la compréhension des processus de réduction de taille dans le cas de produits organiques cristallins. Dans un premier temps, des expériences préliminaires ont été menées pour le choix d'agents mouillants et dispersants adéquats pour les matériaux choisis (ciclopirox, ibuprofène) et pour la mise au point du mode de fonctionnement du procédé de broyage. Nous avons également testé différentes techniques comme la spectroscopie acoustique et la diffusion dynamique de la lumière pour caractériser la distribution de taille des particules broyées, en se basant sur une étude expérimentale conduite sur la calcite, utilisée dans cette étude comme matériau de référence. La seconde partie a été consacrée au broyage par microbilles appliqué aux différents matériaux choisis. Nous avons étudié l'influence des paramètres opératoires tels que le débit de suspension, la vitesse de rotation de l'agitateur, la taille des billes de broyage et la concentration en solide sur l'efficacité énergétique du procédé et la qualité des produits broyés. Les critères de qualité pris en compte dans cette étude sont la distribution de taille des particules, la stabilité et le comportement rhéologique des suspensions broyées ainsi que les propriétés structurales du produit. L'effet du broyage sur les changements microstructuraux des différents produits a enfin été analysé. Une attention particulière a également été portée sur l'effet de la température sur le procédé de réduction de taille et les propriétés de l'ibuprofène broyé. / This work focuses on wet grinding (nanogrinding) of organic materials for pharmaceuticals using a stirred media mill. The aim of the work was to examine the feasibility of nanogrinding process and to improve the understanding of size reduction processes applied on crystalline organic products. Firstly, preliminary investigations were conducted to make a choice of appropriate wetting agents and dispersants for the selected materials (ciclopirox, ibuprofen), and for the development of the operation mode of the grinding process. We also tested different techniques such as acoustic attenuation spectroscopy and dynamic light scattering to characterize the ground particle size distribution, based on an experimental study conducted on calcite, used in this study as reference material. The second part was devoted to the milling process by grinding beads applied to the selected materials. The influence of operating parameters such as suspension flow rate, stirrer tip speed, grinding media diameter and solid mass concentration on grinding process efficiency and ground products quality were investigated. The quality criteria considered in this study are particle size distribution, stability and rheological behaviour of the products, as well as their structural properties. The effect of grinding process on the microstructural change of different materials was finally analysed. Particular attention was also focused on the effect of temperature on size reduction process and ground ibuprofen properties.

Deux approches pour caractériser les milieux poreux en terme de distribution de taille de pores (DTP) sont développées au sein de l'équipe ECPS. Ce travail a pour but de confirmer expérimentalement leurs validités. A l'instar des autres méthodes utilisant l'intrusion du mercure, l'adsorption isotherme ou la thermoporosimétrie, la première méthode consiste à utiliser un fluide à seuil d'écoulement. En effet, l'utilisation de l'écoulement d'un fluide à seuil de type Herschel-Bulkley, au travers d'un poreux, en fonction du gradient de pression permet (en utilisant les solutions analytique et numérique du problème inverse) de déterminer la fonction de distribution de la taille de pores. La seconde méthode utilise l'admittance complexe d'un milieu poreux, mesurée à partir de la réponse en débit à une sollicitation harmonique du gradient de pression. Comme la fréquence de la sollicitation est reliée aux rayons des pores par le biais de la profondeur de pénétration hydrodynamique, l'admittance permet de retrouver la distribution de taille de pores par la résolution numérique du problème inverse associé. Ces deux techniques sont basées sur le modèle de faisceaux de capillaires parallèles employé dans la plupart des autres études qui traitent du même problème. Nos expériences s'appuient sur des milieux poreux calibrés. L'application de ces techniques aux milieux poreux réels se fait actuellement en collaboration avec le TREFLE de Bordeaux. Les résultats expérimentaux obtenus affirment clairement la validité et l'applicabilité de ces deux méthodes pour la caractérisation de la DTP. Il est désormais envisageable de les transférer pour un usage industriel. / Two approaches to characterize porous media in terms of pore size distribution (PSD) are developed within our ECPS team. The aim of this study is to experimentally confirm their validity. Like the other methods using mercury intrusion, adsorption isotherm or thermoporometry, the first method consists in the use of a fluid flow threshold. Indeed, the use of flow of a yield-stress fluid like those of Herschel-Bulkley's, through a porous media, versus the pressure gradient, permits (using the analytical and numerical solutions of the inverse problem) to determine the distribution function of its PSD. The second method uses the complex admittance of a porous medium, measured from the flow rate response to a harmonic pressure gradient. As the frequency of the sollicitation is related to the pore radius through the hydrodynamic penetration depth, the admittance allows to determine the PSD using numerical solution of the associated inverse problem. Both techniques are based on the parallel capillaries bundle model, as used in most of similar studies. Our experiments use calibrated porous media. The application of these techniques to real porous media is currently lead in collaboration with the TREFLE (Bordeaux). Our experimental results clearly confirm the validity and the applicability of these methods for the characterization of the PSD. It is now possible to transfer them for an industrial use.

Les méthodes expérimentales utilisées actuellement pour déterminer la distribution de taille des pores (DTP) dans les milieux poreux présentent des inconvénients, tels que par exemple, la toxicité des fluides employés (porosimétrie à mercure). La base théoriques d'une nouvelle méthode pour obtenir la DTP a été proposée dans la littérature. Celle-ci est fondée sur l'injection de fluides à seuil, caractérisés par une contrainte de cisaillement en deçà de laquelle ils ne s'écoulent pas. L'idée principale de ces travaux théoriques est que l'écoulement de fluides à seuil à travers un milieu poreux permet d'obtenir sa DTP à partir de la mesure des débits correspondant à différents gradients de pression Q(∇P). L'objectif du travail proposé ici est de présenter une nouvelle méthode d'exploitation des données expérimentales Q(∇P) permettant d'obtenir de façon simple, robuste et reproductible les DTPs des milieux poreux analysés. La démarche consiste à évaluer la contribution au débit total des nouveaux pores qui s'incorporent à l'écoulement entre deux valeurs de ∇P. Ces nouveaux pores sont caractérisés par un rayon représentatif qui est fonction de la contrainte seuil du fluide et de ∇P. L'importance de leur contribution au débit total par rapport à celle d'un seul pore donne le nombre de pores dans l'échantillon ayant ce rayon représentatif. Cette méthode est d'abord testée et validée avec des expériences générées numériquement. Ensuite, elle est utilisée pour exploiter des données provenant d'expériences de laboratoire réalisées avec de différents milieux poreux. Les résultats obtenus en termes de DTPs sont comparés avec ceux fournis par d'autres techniques: porosimétrie à mercure et microtomographie. / Current experimental methods used to determine pore size distributions (PSD)of porous media present several drawbacks such as toxicity of the employed fluids (e.g., mercury porosimetry). The theoretical basis of a new method to obtain the PSD by injecting yield stress fluids through porous media and measuring the flow rate Q at several pressure gradients ∇P was proposed in the literature. On the basis of these theoretical considerations,an intuitive approach to obtain PSD from Q(∇P) is presented in this work. It relies on considering the extra increment of Q when ∇P is increased, as a consequence of the pores of smaller radius newly incorporated to the flow. This procedure is first tested and validated on numerically generated experiments. Then, it is applied to exploit data coming from laboratory experiments and the obtained PSDs are compared to those deduced by mercury porosimetry and micro tomography.