Ce travail concerne la généralisation d'un modèle de dépôt de particules. La première partie présente l'approche de modélisation adoptée pour simuler le dépôt de particules colloïdales et traite du raccord entre le Modèle de Langevin Généralisé qui calcule le transport des particules dans le coeur de l'écoulement et le Modèle de Proche Paroi (MPP) destiné au voisinage de la paroi. Le raccord est assuré par l'équilibre du flux de particules traversant l'interface entre les modèles. La validation est assurée par le passage en particules fluides en vérifiant que le raccord est valable dans le cas limite du fluide. Ce travail a donné lieu à des développements théoriques et numériques sur les statistiques conditionnées d'une particule traversant une interface. Des calculs de dépôt de particules montrent que le MPP reproduit correctement les résultats expérimentaux pour des écoulements isothermes. Dans un écoulement non isotherme, les particules peuvent faire l'expérience du phénomène de thermophorèse (qui est le mécanisme de déplacement des particules induit par un gradient de température) qui peut considérablement influencer le dépôt de particules. Nous avons choisi de modéliser la thermophorèse dans le MPP. Nous étudions d'abord la thermophorèse en milieu gazeux car elle est très bien expliquée par la théorie cinétique des gaz. L'étude de la thermophorèse en milieu liquide est plus ardue car elle ne dispose pas encore de théorie unifié ni de travaux expérimentaux qui étudient ce phénomène. La dernière partie est consacrée à la prise en compte du phénomène de thermophorèse en liquide dans le MPP. Cette étude fournit de nouveaux résultats numériques de dépôt en milieu liquide / We present new results obtained with the generalisation to non isothermal flows of a Near-Wall particle deposition Model (NWM). This work is divided in two parts. The first part consists in coupling the Generalised Langevin Model used to compute the particles statistics in the core of the flow with the NWM that calculates particles motion in the vicinity of the wall. In order to achieve the coupling, we focus on the balance of fluxes of fluid particles crossing the interface between the models. This amounts to assessing if the coupling is valid in the fluid limit case, as non-physical accumulation of fluid particles can occur at the interface. It has lead to theoretical and numerical developments on the conditional statistics of particles crossing an interface. The particle deposition rate is then computed with inertial particles and the model reproduces satisfactorily experimental studies in isothermal flow. In the case of non isothermal flows, particles can experience thermophoresis (which is a mechanism that induces a particle flux in the presence of a temperature gradient) that may strongly affect the particle deposition rate. We thus choose to model thermophoresis in NWM to predict accurate particule deposition in non isothermal flows. We first investigate thermophoresis in gas as it is well explained. Thermophoresis in liquids is then studied. In contrast to the gases, the theory and experiment of thermophoresis in liquids are far from being well established. The last part of this work, dedicated to the modelisation of thermophoresis in liquids in the NWM has lead to new promising numerical results of particle deposition enhanced by thermophoresis

La pollution de l'air peut être transportée sur des distances de plusieurs milliers de kilomètres par la dynamique de grande échelle. Elle affecte ainsi le climat et la qualité de l'air loin des régions sources. Durant leur transport, les panaches de pollution voient leur composition altérée par le mélange avec l'air environnant. Les modèles globaux eulériens utilisés pour calculer les bilans d'export de polluants ne représentent pas bien le transport des panaches du fait de leur faible résolution, ce qui induit des erreurs dans le calcul de ces bilans et dans les vitesses de réaction du fait de la non-linéarité de la photochimie atmosphérique. Cette thèse a pour but de contribuer à améliorer la compréhension et la quantification du mélange subi par les panaches de pollution dans la troposphère libre lors du transport à longue distance, et sa représentation dans les modèles numériques. Elle se focalise sur l'Arctique, une région vulnérable au changement climatique, où la distribution des gaz à effet de serre et des aérosols est fortement influencée par le transport de polluants depuis les moyennes latitudes. Les travaux effectués s'organisent autour de trois études de cas de transport dans la troposphère arctique, dans le cadre de la campagne aéroportée POLARCAT durant l'été 2008. Les résultats obtenus permettent d'estimer que l'intensité du mélange dans cette région est plus faible qu'aux moyennes latitudes. Ils montrent également qu'une approche lagrangienne permet de reproduire la variabilité des concentrations dans les panaches âgés, de simuler le transport des structures de petite échelle, et d'évaluer la variabilité de l'évolution des concentrations en ozone dans les panaches induite par les différents scénarios de mélange. Les résultats de ces modèles lagrangiens permettent aussi de jauger les performances des modèles globaux eulériens en terme de transport des polluants en Arctique, et de définir la résolution spatiale nécessaire pour représenter ce transport de façon satisfaisante. Les trois études présentées concordent sur l'importance de combiner les approches lagrangienne et eulérienne pour simuler le transport à longue distance dans les modèles numériques.

Dans ce travail de thèse, réalisé à l'IMFT, nous nous sommes intéressés aux écoulements turbulents diphasiques gaz-solides et plus particulièrement au phénomène de couplage inverse qui correspond à la modulation de la turbulence par la phase dispersée. Ce mécanisme est crucial pour les écoulements à forts chargements massiques. Dans cette thèse, nous avons considéré une turbulence homogène isotrope stationnaire sans gravité dans laquelle des particules sont suivies individuellement d'une façon Lagrangienne. La turbulence du fluide porteur est obtenue par des simulations directes (DNS). Les particules sont sphériques, rigides et d'une taille inférieure aux plus petites échelles de la turbulence. Leur densité est bien plus grande que la densité du fluide. Dans ce cadre, la force la plus importante agissant sur les particules est celle de traînée. Les interactions inter-particules ainsi que la gravité ne sont pas prises en compte. Pour modéliser ce type d'écoulement, une approche stochastique est utilisée pour laquelle l'accélération du fluide est modélisée par une équation de Langevin. L'originalité de ce travail est la prise en compte de l'effet de la modulation de la turbulence par un terme additionnel. Nous avons proposé deux modèles : une force de couplage moyenne qui est définie à partir des vitesses moyennes des phases, et une force instantanée qui est définie à l'aide du formalisme mésoscopique Eulérien. La fermeture des modèles s’appuie sur la fonction d’autocorrélation Lagrangienne et l’équation de transport de l’énergie cinétique. Les modèles sont testés en terme de prédiction de la vitesse de dérive et des corrélations fluide-particule. Les résultats montrent que le modèle moyen, plus simple, prend en compte les effets principaux du couplage inverse. Cependant, le problème de fermeture pratique est reporté sur la modélisation de l’échelle intégrale Lagrangienne et l’énergie cinétique de la turbulence du fluide vue par les particules. / In this thesis, performed in IMFT, we are interested in the turbulent gas-solid flows and more specifically, in the phenomenon of turbulence modulation which is the modification of the structure of the turbulence due to the solid particles. This mechanism is crucial in flows with high particle mass-loadings. In this work, we considered a homogeneous isotropic turbulence without gravity kept stationary with stochastic type forcing. Discrete particles are tracked individually in Lagrangian manner. Turbulence of the carrier phase is obtained by using DNS. The particles are spherical, rigid and of a diameter smaller than the smallest scales of turbulence. Their density is very large in comparison to the density of the fluid. In this configuration the only force acting on the particles is the drag force. Volume fraction of particles is very small and inter-particle interactions are not considered. To model this type of flow, a stochastic approach is used where the fluid element accel- eration is modeled using stochastic Langevin equation. The originality in this work is an additional term in the stochastic equation which integrates the effect of the particles on the trajectory of fluid elements. To model this term, we proposed two types of modeling: a mean drag model which is defined using the mean velocities from the mean transport equations of the both phases and an instantaneous drag term which is written with the help of the Mesoscopic Eulerian Approach. The closure of the models is based on the Lagrangian auto- correlation function of the fluid velocity and on the transport equation of the fluid kinetic energies. The models are tested in terms of the fluid-particle correlations and fluid-particle turbulent drift velocity. The results show that the mean model, simple, takes into account the principal physical mechanism of turbulence modulation. However, practical closure problem is brought forward to the Lagrangian integral scale and the fluid kinetic energy of the fluid turbulence viewed by the particles.

Cette thèse s'inscrit dans un projet de modélisation numérique de la dispersion atmosphérique de polluants à travers le code de mécanique des fluides numérique Code_Saturne. L'objectif est de pouvoir simuler la dispersion atmosphérique de polluants en environnement complexe, c'est-à-dire autour de centrales, sites industriels ou en milieu urbain. Dans ce contexte, nous nous concentrons sur la modélisation de la dispersion des polluants à micro-échelle, c'est-à-dire pour des distances de l'ordre de quelques mètres à quelques kilomètres et correspondant à des échelles de temps de l'ordre de quelques dizaines de secondes à quelques dizaines de minutes : on parle de modélisation en champ proche. L’approche suivie dans ces travaux de recherche suit une formulation hybride eulérienne/lagrangienne, où les champs dynamiques moyens relatifs au fluide porteur (pression, vitesse, température, turbulence) sont calculés via une approche eulérienne et sont ensuite fournis au solveur lagrangien. Ce type de formulation est couramment utilisé dans la littérature atmosphérique pour son efficacité numérique. Le modèle lagrangien stochastique considéré dans nos travaux est le Simplified Langevin Model (SLM), développé par Pope (1985,2000). Ce modèle appartient aux méthodes communément appelées méthodes PDF (Probability Density Function), et, à notre connaissance, n'a pas été exploité auparavant dans le contexte de la dispersion atmosphérique. Premièrement, nous montrons que le SLM respecte le critère dit de mélange homogène (Thomson, 1987). Ce critère, essentiel pour juger de la bonne qualité d'un modèle lagrangien stochastique, correspond au fait que si des particules sont initialement uniformément réparties dans un fluide incompressible, alors elles doivent le rester. Nous vérifions le bon respect du critère de mélange homogène pour trois cas de turbulence inhomogène représentatifs d'une large gamme d'applications pratiques : une couche de mélange, un canal plan infini, ainsi qu'un cas de type atmosphérique mettant en jeu un obstacle au sein d'une couche limite neutre. Nous montrons que le bon respect du critère de mélange homogène réside simplement en la bonne introduction du terme de gradient de pression en tant que terme de dérive moyen dans le modèle de Langevin (Pope, 1987; Minier et al., 2014; Bahlali et al., 2018c). Nous discutons parallèlement de l'importance de la consistance entre champs eulériens et lagrangiens dans le cadre de telles formulations hybrides eulériennes/lagrangiennes. Ensuite, nous validons le modèle dans le cas d'un rejet de polluant ponctuel et continu, en conditions de vent uniforme et turbulence homogène. Dans ces conditions, nous disposons en effet d'une solution analytique nous permettant une vérification précise. Nous observons que dans ce cas, le modèle lagrangien discrimine bien les deux différents régimes de diffusion de champ proche et champ lointain, ce qui n'est pas le cas d'un modèle eulérien à viscosité turbulente (Bahlali et al., 2018b).Enfin, nous travaillons sur la validation du modèle sur plusieurs campagnes expérimentales en atmosphère réelle, en tenant compte de la stratification thermique de l'atmosphère et de la présence de bâtiments. Le premier programme expérimental considéré dans nos travaux concerne le site du SIRTA (Site Instrumental de Recherche par Télédétection Atmosphérique), dans la banlieue sud de Paris, et met en jeu une stratification stable de la couche limite atmosphérique. La seconde campagne étudiée est l'expérience MUST (Mock Urban Setting Test). Réalisée aux Etats-Unis, dans le désert de l'Utah, cette expérience a pour but de représenter une ville idéalisée, au travers d'un ensemble de lignées de conteneurs. Deux rejets ont été simulés et analysés, respectivement en conditions d'atmosphère neutre et stable (Bahlali et al., 2018a) / This Ph.D. thesis is part of a project that aims at modeling pollutant atmospheric dispersion with the Computational Fluid Dynamics code Code_Saturne. The objective is to simulate atmospheric dispersion of pollutants in a complex environment, that is to say around power plants, industrial sites or in urban areas. In this context, the focus is on modeling the dispersion at micro-scale, that is for distances of the order of a few meters to a few kilometers and corresponding to time scales of the order of a few tens of seconds to a few tens of minutes: this is also called the near field area. The approach followed in this thesis follows a hybrid Eulerian/Lagrangian formulation, where the mean dynamical fields relative to the carrier fluid (pressure, velocity, temperature, turbulence) are calculated through an Eulerian approach and are then provided to the Lagrangian solver. This type of formulation is commonly used in the atmospheric literature for its numerical efficiency. The Lagrangian stochastic model considered in our work is the Simplified Langevin Model (SLM), developed by Pope (1985,2000). This model belongs to the methods commonly referred to as PDF (Probability Density Function) methods, and, to our knowledge, has not been used before in the context of atmospheric dispersion. First, we show that the SLM meets the so-called well-mixed criterion (Thomson, 1987). This criterion, essential for any Lagrangian stochastic model to be regarded as acceptable, corresponds to the fact that if particles are initially uniformly distributed in an incompressible fluid, then they must remain so. We check the good respect of the well-mixed criterion for three cases of inhomogeneous turbulence representative of a wide range of practical applications: a mixing layer, an infinite plane channel, and an atmospheric-like case involving an obstacle within a neutral boundary layer. We show that the good respect of the well-mixed criterion lies simply in the good introduction of the pressure gradient term as the mean drift term in the Langevin model (Pope, 1987; Minier et al., 2014; Bahlali et al., 2018c). Also, we discuss the importance of consistency between Eulerian and Lagrangian fields in the framework of such Eulerian/Lagrangian hybrid formulations. Then, we validate the model in the case of continuous point source pollutant dispersion, under uniform wind and homogeneous turbulence. In these conditions, there is an analytical solution allowing a precise verification. We observe that in this case, the Lagrangian model discriminates well the two different near- and far-field diffusion regimes, which is not the case for an Eulerian model based on the eddy-viscosity hypothesis (Bahlali et al., 2018b).Finally, we work on the validation of the model on several experimental campaigns in real atmosphere, taking into account atmospheric thermal stratification and the presence of buildings. The first experimental program considered in our work has been conducted on the `SIRTA' site (Site Instrumental de Recherche par Télédétection Atmosphérique), in the southern suburb of Paris, and involves a stably stratified surface layer. The second campaign studied is the MUST (Mock Urban Setting Test) experiment. Conducted in the United States, in Utah's desert, this experiment aims at representing an idealized city, through several ranges of containers. Two cases are simulated and analyzed, respectively corresponding to neutral and stable atmospheric stratifications (Bahlali et al., 2018a)

L'objectif de cette thèse sur la modélisation des décanteurs automobiles est de mettre au point un outil de simulation numérique permettant de se substituer à un banc organe pour déterminer l'efficacité des décanteurs à l'étape de conception. Après avoir étudié la composition du blow-by sur banc moteur et cela pour différents points de fonctionnement du moteur, nous avons pu comparer la pertinence des différents outils de mesure granulométrique choisis (PCS, ELPI, TEOM) en fonction du spectre granulométrique émis par le moteur. Le PCS a été retenu comme l'outil le plus adéquat en fonction des tailles de gouttes sur lesquelles nous souhaitions avoir des informations. En parallèle, nous avons mené des investigations pour mettre au point une méthodologie de calcul permettant de simuler l'efficacité du décanteur. Le suivi des particules d'huile s'effectue grâce à une simulation lagrangienne. Les forces à prendre en compte ont été déterminées grâce à des expériences numériques simples permettant de verifier si les différentes expressions de ces forces trouvées dans la littérature étaient compatibles des conditions de calculs rencontrées dans les décanteurs. Une autre question soulevée par l'utilisation d'une modélisation lagrangienne dans la littérature, est celle de la prise en compte ou non de la dispersion turbulente. L'étude proposée dans cette thèse nous a permis de montrer que l'implémentation d'un modèle de dispersion turbulente rudimentaire ne permet pas d'obtenir des résultats satisfaisants. La faisabilité de la mise au point d'une méthodologie de calcul fiable de l'efficacité du décanteur repose finalement sur la modification du type de maillage utilisé ainsi que le traitement du calcul de l'écoulement à la paroi. Nous avons montré dans cette thèse que la mise en place d'un maillage raffiné aux parois reposant sur des prismes avec un Y+ =1, associé à une loi de parois à deux couches permettait d'obtenir des résultats satisfaisants en terme de prédiction de l'efficacité du décanteur.

Nous présentons la généralisation d'un modèle de dépôt de particules aux écoulements non isothermes. Ce travail comporte deux parties. La première partie présente l'approche de modélisation lagrangienne stochastique adoptée pour simuler le dépôt de particules colloïdales et traite du raccord entre le Modèle de Langevin Généralisé qui calcule le transport des particules dans le cœur de l'écoulement et le Modèle de Proche Paroi destiné au transport des particules au voisinage de la paroi. Le raccord est assuré par l'équilibre du flux de particules traversant l'interface entre les deux modèles. La validation du couplage s'effectue par le passage en particules fluides, il s'agit de vérifier que le raccord est valable dans le cas limite du fluide. En effet, des accumulations non physiques de particules fluides peuvent se produire autour de l'interface. Le travail de couplage a donné lieu à des développements théoriques et numériques sur les statistiques conditionnées d'une particule traversant une interface. Des calculs de dépôt de particules montrent que le modèle reproduit correctement les résultats expérimentaux pour des écoulements isothermes. Cependant, dans un écoulement non isotherme, les particules peuvent faire l'expérience du phénomène de thermophorèse (qui est le mécanisme de déplacement des particules induit par un gradient de température) qui peut considérablement influencer le dépôt de particules. Nous avons donc choisi d'intégrer la modélisation de la thermophorèse dans le Modèle de Proche Paroi afin d'évaluer au mieux le taux de dépôt de particules dans un écoulement non isotherme. Nous commençons par étudier la thermophorèse en milieu gazeux car elle est très bien expliquée par la théorie cinétique des gaz. L'étude de la thermophorèse en milieu liquide est plus ardue car elle ne dispose pas encore de théorie unifié ni de travaux expérimentaux qui étudient ce phénomène. La dernière partie de ce document est donc consacrée à la prise en compte du phénomène de thermophorèse en liquide dans le Modèle de Proche Paroi. Cette étude fournit de nouveaux résultats numériques de dépôt de particules en présence d'un gradient de température dans un milieu liquide.

La fermeture du bilan d'ozone troposphérique nécessite de mieux <br />connaître l'importance du transport vertical depuis les principaux <br />réservoirs: la couche limite atmosphérique et la stratosphère. Les <br />études de cas permettent de mettre en évidence les mécanismes entrant en <br />jeu et leur efficacité, mais leur impact sur le bilan ne peut-être <br />évalué qu'à partir d'analyses climatologiques. L'exposé est donc divisé <br />en deux grandes parties.<br />La première porte sur des études de cas à partir des mesures LIDAR ozone <br />effectuées pendant la campagne ESCOMPTE. La méthodologie repose sur une <br />analyse des données collectées par une approche hybride lagrangienne <br />associant des analyses météorologiques globales, et des modèles de <br />dynamique méso-échelle, de chimie-transport et de photochimie <br />lagrangienne. L'objectif étant d'identifier la part de la variabilité <br />observée imputable aux processus transport et, le cas échéant, à la <br />production photochimique d'ozone dans la troposphère.<br />Dans une deuxième partie, nous proposons une évaluation de l'importance <br />respective du transport depuis la couche-limite et depuis la région de <br />la tropopause sur le bilan d'ozone troposphérique à travers une analyse <br />climatologique. Un algorithme multivariable de classification des <br />couches d'ozone mesurées par sondage est présenté, et comparé à une <br />approche trajectographique. Cet algorithme est ensuite appliqué aux <br />données collectées depuis une trentaine d'années en Europe afin <br />d'analyser la variabilité saisonnière, géographique et interannuelle des <br />processus de transport vertical ainsi que leur influence sur le bilan <br />d'ozone troposphérique. La variabilité du mélange turbulent dans la <br />troposphère et ses conséquences sur la persistance des couches sera <br />aussi évoquée.

Les tortues marines, espèces emblématiques des écosystèmes marins, sont de plus en plus menacées par les effets directs et indirects des activités humaines. Leur cycle de vie est complexe, partagé entre divers habitats, souvent très éloignés les uns des autres. Leur conservation nécessite donc d'identifier les habitats occupés à chaque stade de vie et les routes migratoires empruntées entre ces différents habitats. Si l'écologie spatiale des tortues adultes est relativement bien connue, notamment grâce au suivi par satellite, il n'en va pas de même pour les juvéniles qui se développent plusieurs années en milieu pélagique sans pouvoir être suivis. Dans ce contexte, les simulations numériques constituent un outil adapté pour explorer la dispersion des tortues juvéniles à partir de leurs plages de naissance. Jusqu'à présent il a le plus souvent été supposé dans ces simulations que les juvéniles dérivaient passivement avec les courants marins. Dans ce travail de thèse nous présentons STAMM (Sea Turtle Active Movement Model), un nouveau modèle de dispersion active des tortues juvéniles qui s'attache à dépasser l'hypothèse initiale d'une dérive purement passive. Dans STAMM, les juvéniles simulés se déplacent sous l'influence de la circulation océanique et d'une nage motivée par la recherche d'habitats favorables. Ce modèle est appliqué ici à l'étude de la dispersion des juvéniles de trois populations de tortues marines : les tortues luths (Dermochelys coriacea) du Pacifique Ouest et de l'Atlantique Nord-Ouest puis les tortues caouannes (Caretta caretta) de l'ouest de l'océan Indien. Nos résultats montrent que, même si la circulation océanique détermine, à grande échelle, les zones de dispersion, la prise en compte des mouvements motivés par l'habitat augmente considérablement le réalisme des simulations et impacte profondément la distribution spatiale et temporelle des individus simulés à l'intérieur de leur zone de dispersion. Les mouvements motivés par l'habitat induisent notamment des migrations saisonnières en latitude qui réduisent la mortalité par hypothermie. Ces mouvements induisent également une concentration des individus simulés dans des zones productives (comme les upwellings de bord Est) inaccessibles en dérive passive. Ces résultats questionnent la vision classique des juvéniles circulant passivement autour des gyres océaniques et devraient rapidement être pris en compte pour la mise en place de mesures de conservation ciblées visant les tortues marines juvéniles. / Sea turtles are increasingly threatened by the direct and indirect effects of human activities. Their life cycle is complex, shared between various, and often very distant, habitats. Their conservation therefore requires identifying the habitats occupied at each stage of life and the migration routes between these different habitats. While the spatial ecology of adult turtles is relatively well known, particularly through satellite monitoring, the situation is not the same for juveniles which pelagic development phase remains largely unobserved. In that context, numerical simulation constitutes an appropriate tool to explore the dispersal of juvenile sea turtles from their natal beaches. Until now, simulations were mostly performed under the assumption that juveniles disperse passively with oceanic currents. In this PhD thesis we present STAMM (Sea Turtle Active Movement Model), a new model of active dispersal that aims to go beyond the initial hypothesis of passive drift. In STAMM, juvenile sea turtles move under the influence of ocean currents and swimming movements motivated by the search for favorable habitats. This model is applied here to the study of the dispersal of juveniles from three sea turtle populations: leatherback turtles (Dermochelys coriacea) of the Western Pacific and the Northwest Atlantic Oceans, and loggerhead turtles (Caretta caretta) of the Western Indian Ocean. Our results show that, although ocean currents broadly shape juvenile dispersal areas, simulations including habitat-driven movements provide more realistic results than passive drift simulations. Habitat-driven movements prove to deeply structure the spatial and temporal distribution of juveniles. In particular, they induce seasonal latitudinal migrations that reduce cold induce mortality. They also push simulated individuals to concentrate in productive areas that cannot be accessed through pure passive drift. These results challenge the classical view of juveniles circulating passively around oceanic gyres. They should rapidly be taken into account for the implementation of targeted conservation measures concerning juvenile sea turtles.

La 3D s'impose comme un nouveau média dont l'adoption généralisée passe par la conception d'outils, accessibles au grand public, de création et de manipulation de formes tridimensionnelles quelconques. Les outils actuels reposent fortement sur la modélisation sous-jacente des formes, généralement surfacique, et sont alors peu intuitifs ou limitatifs dans l'expressivité offerte à l'utilisateur.Nous souhaitons, dans ces travaux, définir une approche ne présentant pas ces défauts et permettant à l'utilisateur de se concentrer sur le processus créatif. En nous inspirant de l'utilisation séculaire de l'argile, nous proposons une approche modélisant la matière sous forme lagrangienne.Une forme est ainsi décrite par un système de particules, où chaque particule représente un petit volume du volume global.Dans ce cadre lagrangien, la méthode Smoothed Particle Hydrodynamics (SPH) permet l'approximation de grandeurs physiques en tout point de l'espace. Nous proposons alors une modélisation de matériaux à deux couches, l'une décrivant la topologie et l'autre décrivant la géométrie du système global.La méthode SPH permet, entre autres, d'évaluer la densité de matière. Ceci nous permet de définir une surface implicite basée sur les propriétés physiques du système de particules pour redonner un aspect continu à la matière.Ces matériaux peuvent alors être manipulés au moyen d'interactions locales reproduisant le maniement de la pâte à modeler, et de déformations globales. L'intérêt de notre approche est démontrée par plusieurs prototypes fonctionnant sur des stations de travail standard ou dans des environnements immersifs. / 3D is emerging as a new media. Its widespread adoption requires the implementation of userfriendly tools to create and manipulate three-dimensional shapes. Current softwares heavily rely on underlying shape modeling, usually a surfacic one, and are then often counter-intuitive orlimiting. Our objective is the design of an approach alleviating those limitations and allowing the user to only focus on the process of creating forms. Drawing inspiration from the ancient use of clay,we propose to model a material in a lagrangian description. A shape is described by a particles system, where each particle represents a small fraction of the total volume of the shape. In this framework, the Smoothed Particle Hydrodynamics method enables to approximate physical values anywhere in space. Relying on this method, we propose a modeling of material with two levels, one level representing the topology and the other one describing local geometry of the shape.The SPH method especially enables to evaluate a density of matter. We use this property todefine an implicit surface based on the physical properties of the particles system to reproduce the continuous aspect of matter. Those virtual materials can then be manipulated locally through interactions reproducing the handling of dough in the real world or through global shape deformation. Our approach is demonstrated by several prototypes running either on typical desktop workstation or in immersive environment system.