Le contexte général de cette thèse s'inscrit dans le cadre de la modélisation eulérienne instationnaire des écoulements turbulents anisothermes gaz - particules. La modélisation de ces écoulements est cruciale pour de nombreuses applications industrielles et pour la prédiction de certains phénomènes naturels. Par exemple, la combustion diphasique dans les moteurs automobiles et aéronautiques est précédée par l'injection et la dispersion de carburant liquide dans la chambre de combustion. Les phénomènes mis en jeu exigent alors une prédiction locale tenant compte du caractère instationnaire de l'écoulement turbulent et de la présence de géométries complexes. De plus, de nombreuses études expérimentales et numériques récentes ont mis en évidence le rôle prépondérant de l'inertie des particules sur les mécanismes de dispersion et de concentration préférentielle en écoulement turbulent. Ceci rend donc indispensable la prise en compte de ces mécanismes dans la modélisation diphasique. Au cours de ce travail de thèse, une approche eulérienne locale et instantanée a été développée pour prédire les écoulements gaz-particules anisothermes et turbulents. Elle est basée sur l'approche statistique du Formalisme Eulérien Mésoscopique (MEF) introduite par Février et al. (JFM, 2005). Cette approche a été ici étendue aux variables thermiques pour la prise en compte du caractère anisotherme de l'écoulement. Cette approche a été ensuite utilisée dans le cadre de la méthode des moments (Kaufmann et al., JCP, 2008), et un système d'équations locales et instantanées pour la phase dispersée a été proposé. La modélisation au premier ordre exige la fermeture des moments de second ordre apparaissant dans les équations de la quantité de mouvement et de l'énergie. La proposition de telles relations constitutives fait l'objet d'une partie de la thèse. Afin de fournir une méthode capable de prédire le comportement local, instantané et anisotherme de la phase dispersée dans des configurations `a une échelle réaliste, les équations pour la phase dispersée ont été filtrées et une modélisation aux grandes échelles (LES) est effectuée. Cette modélisation étends, par la prise en compte des variables thermiques, le travail de Moreau et al. (FtaC, 2010) sur l'approche LES Euler-Euler en conditions isothermes. L'approche complète est enfin appliquée aux résultats de simulation numérique d'un jet plan turbulent gazeux froid, chargé en particules, dans une turbulence homogène isotrope chaude monophasique. / The aim of this thesis is to provide an Eulerian modeling for the dispersed phase interacting with unsteady non-isothermal turbulent flows. The modeling of these flows is crucial for several industrial applications and for predictions of natural events. Examples are the combustion chambers of areo engines where the combustion is preceded by the injection and dispersion of liquid fuel. The prediction of such phenomena involves a local modeling of the mixture for taking into account the unsteady behavior of the turbulent flow and the presence of complex geometries. Moreover, many experimental and numerical studies have recently highlighted the significant role of the particle inertia on the mechanisms of dispersion and preferential concentration. Accounting for such mechanisms is therefore essential for modeling the particle-laden turbulent flows. In this thesis, a local and instantaneous Eulerian approach able to describe and to predict the local behavior of inertial particles interacting with non-isothermal turbulent flows has been developed. It is based on the statistical approach known as Mesoscopic Eulerian formalism (MEF) introduced by Février et al. (JFM, 2005). The statistical approach has been extended to the thermal quantities in order to account for the non-isothermal conditions into the modeling. This formalism is then used in the framework of the moment approach (Kaufmann et al., JCP, 2008) and a system of local and instantaneous equations for the non-isothermal dispersed phase has been suggested. The first order modeling requires to close second-order moments appearing in momentum and energy equations. The proposal of such constitutive relations makes the object of a part of this study. In order to provide an Eulerian approach usable in real configurations at industrial scale, the equations of the dispersed phase are filtered and the approach developed in the framework of the Large-Eddy Simulations. From the work of Moreau et al. (FTaC, 2010), the Eulerian-Eulerian LES approach is then extended to non-isothermal conditions. The whole modeling is then a priori tested against numerical simulations of a cold planar turbulent particle-laden jet crossing a homogeneous isotropic decaying hot turbulence.

An experimental study on impaction of a single drop on solid surfaces was conducted to show the effects of particles on the impact process. The parameters were: volume fraction of particles (0 to 0.3), particle size (0.47 to 250 micron), and ratio of particle size to drop size (0.00017 to 0.074). The effect of particle volume fraction on the spreading process depended on impact speed and substrate. At low impact speed, particles had little effect on the spreading except for surfaces where the equilibrium contact angle was low. For high impact speed, the influence of particles on spreading can be described by the effective viscosity. The effect of particle size on the spreading process also depended on impact speed and substrate. At low impact speed, the drop did not have enough kinetic energy to overcome the energy barrier associated with the large particles. For particle-laden liquids, retraction was affected by particle parameters. When pure liquid drops retracted from the maximum spreading ratio, the retraction appeared to be symmetric around the point of impaction while retraction of the particle-laden drop was sometimes asymmetric. Rebounding on the Teflon film depended on impact speed, particle volume fraction, and particle size. The impact speed must reach a critical value for rebounding to occur. Bouncing results suggested that the probability of bouncing decreased as viscosity increased, impact speed increased, and surface tension decreased. The non-wetting behavior and bouncing probably involved an air layer between the surface and the drop. When a low-velocity liquid drop impacts on a surface, ejection of a secondary drop from the top of the impacting drop was sometimes observed. When Renardy et al.'s (2003) criterion for the range of velocities for existence of a capillary wave was applied to for a 3.2-mm water drop; the range was between 0.2 to 1.5 m/s. However, drop ejection was observed at lower impact speed. When apparent viscosity of the particle-laden liquid obtained from Krieger's equation (1972) was used in the pure liquid models for predicting the maximum spreading ratio, good agreement between model predictions and experimental results was obtained when Park et al's model (2003) was used.

Le but de cette thèse est d'étudier le transport et le dépôt de particules solides dans les écoulements à travers les fractures. Dans un premier temps, l'écoulement monophasique à travers les fractures est étudié afin d'évaluer la validité de la loi cubique locale comme modèle de l'écoulement. Des canaux à parois sinusoïdales à géométrie variable sont utilisés pour représenter différents types de fractures. Un premier développement analytique montre que l'ouverture hydraulique de la fracture diffère de son ouverture moyenne lorsque la rugosité des parois est élevée. La méthode des éléments finis est ensuite utilisée pour résoudre les équations de continuité et de Navier-Stokes et comparer les solutions numériques aux prédictions théoriques de la loi cubique locale sur une gamme relativement étendue de nombres de Reynolds Re. Pour de faibles Re, typiquement inférieurs à 15, la loi cubique locale décrit raisonnablement l'écoulement, surtout lorsque la rugosité et le déphasage entre les parois sont relativement faibles. Dans un deuxième temps, les écoulements chargés de particules sont étudiés. Une approche analytique est d'abord développée pour montrer comment des particules distribuées dans un écoulement stationnaire et laminaire à travers une fracture peuvent être transportées sur de longues distances ou au contraire se déposer à l'intérieur. Plus précisément, une équation simple décrivant la trajectoire d'une particule est établie. Sur la base de cette équation, il est démontré que, quand l'inertie des particules est négligeable, leur comportement dépend directement de la géométrie de la fracture et d'un nombre adimensionnel W qui relie la vitesse de sédimentation des particules à la vitesse moyenne de l'écoulement. L'équation proposée est vérifiée en comparant ses prédictions à des simulations numériques de suivi de particules prenant en compte l'inertie des particules et résolvent complètement les équations de Navier-Stokes. Il est montré que l'équation est valide lorsque l'inertie du fluide est faible. Des diagrammes de régimes, permettant de prévoir le comportement des particules à travers la fracture sont proposés. Enfin, un appareil expérimental conçu dans le but d'effectuer une évaluation pratique du modèle analytique est présenté et les résultats préliminaires sont discutés. Les résultats expérimentaux préliminaires tendent valider le modèle analytique. De façon plus générale, les résultats obtenus à travers ce travail de thèse font progresser nos connaissances du comportement des petites particules transportées dans les écoulements de fractures. Potentiellement, ce travail devrait permettre d'améliorer notre prévision de la pollution souterraine, et peut avoir des applications dans le développement de nouvelles techniques de filtration de l'eau et de séparation des minéraux / The aim of the present thesis is to study the transport and deposition of small solid particles in fracture flows. First, single-phase fracture ow is investigated in order to assess the validity of the local cubic law for modeling ow in corrugated fractures. Channels with sinusoidal walls having different geometrical properties are considered to represent different fracture geometries. It is analytically shown that the hydraulic aperture of the fracture clearly deviates from its mean aperture when the walls roughness is relatively high. The finite element method is then used to solve the continuity and the Navier-Stokes equations and to simulate fracture ow in order to compare with the theoretical predictions of the local cubic law for Reynolds numbers Re in the range 0.067-67. The results show that for low Re, typically less than 15, the local cubic law can properly describe the fracture ow, especially when the fracture walls have small corrugation amplitudes. For Re higher than 15, the local cubic law can still be valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitude, and moderate phase lag between its walls. Second, particle-laden flows are studied. An analytical approach has been developed to show how particles sparsely distributed in steady and laminar fracture flows can be transported for long distances or conversely deposited inside the channel. More precisely, a rather simple particle trajectory equation is established. Based on this equation, it is demonstrated that when particles' inertia is negligible, their behavior is characterized by the fracture geometry and by a dimensionless number W that relates the ratio of the particles sedimentation terminal velocity to the ow mean velocity. The proposed particle trajectory equation is verified by comparing its predictions to particle tracking numerical simulations taking into account particle inertia and resolving the full Navier-Stokes equations. The equation is shown to be valid under the conditions that ow inertial effects are limited. Based on this trajectory equation, regime diagrams that can predict the behavior of particles entering closed channel flows are built. These diagrams enable to forecast if the particles entering the channel will be either deposited or transported till the channel outlet. Finally, an experimental apparatus that was designed to have a practical assessment of the analytical model is presented. Preliminary experimental results tend to verify the analytical model. Overall, the work presented in this thesis give new insights on the behavior of small particles in fracture flows, which may improve our prediction and control of underground contamination, and may have applications in the development of new water filtration and mineral separation techniques

Les écoulements chargés en particules sont présents dans de nombreuses applications industrielles telles que le transport de boues ou l’industrie chimique en général. Dans des mélanges constitués de particules solides immergées dans un fluide visqueux, les interactions entre particules jouent un rôle essentiel dans la viscosité globale du mélange.Le phénomène de suspension est causé par des interactions hydrodynamiques à courte distance, connues sous le nom de lubrification. Les forces de lubrification sont généralement sous-estimées en raison de leur nature et de la discrétisation spatiale du problème.Dans cette thèse, nous proposons un modèle de lubrification qui estime les forces et couples hydrodynamiques non résolues par un solveur couplant la résolution des équations de Navier-Stokes incompressible par une méthode de volumes pénalisés, à la résolution de la dynamique des particules par une méthode aux éléments discrets. Les corrections des contraintes hydrodynamiques sont faites localement sur la surface des particules en interaction sans aucune hypothèse sur la forme générale des particules. La version finale du modèle de lubrification proposée peut être utilisée pour des suspensions de particules convexes sans aucune tabulation. La méthode numérique a été validée avec des particules sphériques et des ellipsoïdes, en comparant des simulations à des données expérimentales.Dans le cas de particules sphériques, le modèle de lubrification est aussi précis que les modèles de lubrification existants qui sont limités à ce type de géométrie. La compatibilité du modèle avec des particules convexes a été validée en comparant des simulations,utilisant des ellipsoïdes, à des mesures expérimentales que nous avons réalisées. / Particle-laden flows can be found in many industrial applications such as slurry transport or the chemical industry in general. In mixtures made of solid particles emerged in a viscous fluid, particle interactions play an essential role in the overall mixture viscosity. The suspension phenomenon is caused by short-range hydrodynamic interactions, known as lubrication. Lubrication forces are usually underestimated due to their singularities and the spatial discretization of the numerical schemes. In this thesis, we propose a lubrication model for a coupled volume penalization method and discrete element method solver that estimates the unresolved hydrodynamic forces and torques in incompressible Navier-Stokes flows. Corrections are made locally on the surfaces of the interacting particles without any assumption on the global particle shapes. The final version of the local lubrication model can be used for suspension of convex particles without any tabulations. The numerical method has been validated against experimental data with spherical and ellipsoidal particles. With spherical particles, the lubrication model performs as well as existing numerical models that are limited to this specific particle shape. The model compatibility with convex particles has been validated by comparing simulations using ellipsoids to experimental measurements we made.

Objetivamos analisar a construção do conceito de inimigo representado nos discursos de Osama bin Laden, durante o período de 1996 a 2004. Estes inimigos são constituídos a partir da edificação de uma realidade que apresenta Estados Unidos, Israel e seus aliados, como opressores dos muçulmanos. A fim de melhor compreender esta construção simbólica, buscaremos estudá-la com base na dinâmica de dois eixos fundamentais para a constituição do pensamento político-religioso de bin Laden: a herança teórica do islamismo radical, conjugada com a sua percepção acerca dos principais eventos geopolíticos contemporâneos envolvendo o mundo muçulmano. Dentro desta perspectiva, empregaremos o construtivismo em sua corrente linguística dos autores Nicholas Onuf e François Debrix como arcabouço teórico. Nos guiaremos pela prerrogativa de que a compreensão do mundo e da realidade ocorre de acordo com o modo como nós nos referimos a eles a partir de nossa linguagem. Como referencial metodológico que orientará nosso exame da construção do conceito de inimigo nos discursos de bin Laden, utilizaremos as categorias de operação da ideologia propostas por John Thompson na obra “Ideologia e Cultura Moderna – Teoria Social Crítica na Era dos Meios de Comunicação”. / We aimed to analyze the construction of the concept of enemy represented in Osama bin Laden‟s speeches, during the period 1996 to 2004. These enemies are making from the building of a reality that presents U.S., Israel and its allies, as oppressors of muslims. To better understand this symbolic construction, we will seek to study it based on the dynamics of two pillars for the constitution of bin Laden's politico-religious thought: the theoretical heritage of radical islam, combined with the author's perception about the main events contemporary geopolitical involving the muslim world. From this perspective, we will employ the Linguistic Constructivism of the authors Nicholas Onuf and François Debrix, as theoretical framework. We will be guided by the prerogative of the understanding of the world and reality occurs according to the way we refer to them from our language. As a methodological framework that will guide our examination of the construction of the concept of enemy in bin Laden's speeches, we will use the categories of operation of ideology proposed by John Thompson in his work "Ideology and Modern Culture - Critical Social Theory in the Era of the Media."

Wetgewing wat onbestaanbaar is met die bepalings van Wet 200 van 1993 as die hoogste reg is nietig tensy geoorloof volgens die maatstawwe van art 33 (1) 1 en die Konstitusionele Hof en die provinsiale en plaaslike afdelings van die Hooggeregshof is bevoeg om dit ongeldig te verklaar. Die beoordeling van die geldigheid van wetgewing waardeur inbreuk gemaak word op 'n Fundamentele Reg of Vryheid in Hoofstuk 3 verleen 1 vereis in eerste instansie die vertolking van die betrokke bepaling van die Hoofstuk. In die lig van die besondere karakter van 'n grondwet en die gebrek aan omlyning van die verleende regte en vryhede 1 moet die vertolking nie volgens geykte metodes geskied nie I maar wel doelgerig aan die hand van die waardes onderliggend tot die Handves as geheel. Waardes wat nie herleibaar is tot die grondwet self nie en die persoonlike filosofiee van die regters behoort die howe egter nie te beinvloed nie. / Legislation inconsistent with the provisions of Act 200 of 1993 as the supreme law is of no force and effect and the Constitutional Court and the provincial and local divisions of the Supreme Court may declare it invalid, unless it complies with the criteria of a permissible limitation under sec 33(1). Before the validity of legislation impinging upon a Fundamental Right or Freedom conferred in Chapter 3 can be considered, the relevant provision of the Chapter must first be interpreted. In view of the special nature of a constitution and the generality of the conferred rights and freedoms, the ordinary methods of interpretation should not be applied. Chapter 3 must be interpreted purposively taking into account the values underlying the Bill of Rights as a whole. Values not founded in the Constitution and the personal philosophies of the judges ought not to influence the courts. / Private Law / LL.M.

abstract: This research is particularly concerned with organizations’ advocacy of value-based change aimed at improving consumers’ well-being. This work contributes to the Transformative Services Research area and presents a conceptualization of the value-laden service organization (VLSO), which I define as organizations that advocate for specific value-based behaviors from consumers both within and beyond the particular service setting. In a VLSO, consumers are expected to act in accordance with the values of the organization. If the consumer’s pre-existing value system is not aligned with the values of the service organization, the consumer may experience a sense of psychological disequilibrium, which can lead to unintended decrease in well-being. This research explores how value conflicts are managed by both the organization and by the consumers. This work emerges out of an interpretive study of a Catholic-based homeless shelter for pregnant women. From it, I identify the practices of consumers and the service organization and explored their interactions. This has resulted in a theoretical conceptualization of a Rescue Institution, which combines aspects of both a Total Institution and a Reinventive Institution in a unique way. Further, I conceptualize a cycle of agency and authenticity that maps the dynamics of the consumer in a VLSO as they negotiate the structure/agency duality. In gathering data, I used an interpretive approach over the course of three years’ of direct involvement with a service organization, St. Mary’s House. My methods included participant observation, collection of artifacts, and one-on-one in-depth interviews. I interviewed a total of 30 participants, whose transcribed interviews resulted in over 1500 pages of text. Analysis of themes and concepts occurred as a result of repeated examinations of both existing theory and data. My findings reveal key organizational and consumer practices that negotiate the tension between structure and agency. Organizational practices include rules and social norms, as well as two forms of hierarchy. Consumer practices, often in response to organizational practices, include a cycle of agency and authenticity and participation in a shadow structure. These practices collectively influence consumer’s interpretive drift, which is their adoption of the organization’s values that creates internalized change. I conclude with implications for theory and service organization management. First, value priorities mean that tradeoffs must be made, which can cause unexpected and painful conflict. The experience of change, from both the consumer and service provider perspective, can be very messy. This process includes a dynamic and individual negotiation of authenticity and agency, which will be of interest in future studies. The service providers must be open to this process, carefully navigating their responses to the consumer’s dynamic authenticity, agency and values. Service providers should expect and acknowledge the conflict in consumers’ experience in order to foster their long-term perspective and perseverance. / Dissertation/Thesis / Doctoral Dissertation Business Administration 2015

Objetivamos analisar a construção do conceito de inimigo representado nos discursos de Osama bin Laden, durante o período de 1996 a 2004. Estes inimigos são constituídos a partir da edificação de uma realidade que apresenta Estados Unidos, Israel e seus aliados, como opressores dos muçulmanos. A fim de melhor compreender esta construção simbólica, buscaremos estudá-la com base na dinâmica de dois eixos fundamentais para a constituição do pensamento político-religioso de bin Laden: a herança teórica do islamismo radical, conjugada com a sua percepção acerca dos principais eventos geopolíticos contemporâneos envolvendo o mundo muçulmano. Dentro desta perspectiva, empregaremos o construtivismo em sua corrente linguística dos autores Nicholas Onuf e François Debrix como arcabouço teórico. Nos guiaremos pela prerrogativa de que a compreensão do mundo e da realidade ocorre de acordo com o modo como nós nos referimos a eles a partir de nossa linguagem. Como referencial metodológico que orientará nosso exame da construção do conceito de inimigo nos discursos de bin Laden, utilizaremos as categorias de operação da ideologia propostas por John Thompson na obra “Ideologia e Cultura Moderna – Teoria Social Crítica na Era dos Meios de Comunicação”. / We aimed to analyze the construction of the concept of enemy represented in Osama bin Laden‟s speeches, during the period 1996 to 2004. These enemies are making from the building of a reality that presents U.S., Israel and its allies, as oppressors of muslims. To better understand this symbolic construction, we will seek to study it based on the dynamics of two pillars for the constitution of bin Laden's politico-religious thought: the theoretical heritage of radical islam, combined with the author's perception about the main events contemporary geopolitical involving the muslim world. From this perspective, we will employ the Linguistic Constructivism of the authors Nicholas Onuf and François Debrix, as theoretical framework. We will be guided by the prerogative of the understanding of the world and reality occurs according to the way we refer to them from our language. As a methodological framework that will guide our examination of the construction of the concept of enemy in bin Laden's speeches, we will use the categories of operation of ideology proposed by John Thompson in his work "Ideology and Modern Culture - Critical Social Theory in the Era of the Media."

Objetivamos analisar a construção do conceito de inimigo representado nos discursos de Osama bin Laden, durante o período de 1996 a 2004. Estes inimigos são constituídos a partir da edificação de uma realidade que apresenta Estados Unidos, Israel e seus aliados, como opressores dos muçulmanos. A fim de melhor compreender esta construção simbólica, buscaremos estudá-la com base na dinâmica de dois eixos fundamentais para a constituição do pensamento político-religioso de bin Laden: a herança teórica do islamismo radical, conjugada com a sua percepção acerca dos principais eventos geopolíticos contemporâneos envolvendo o mundo muçulmano. Dentro desta perspectiva, empregaremos o construtivismo em sua corrente linguística dos autores Nicholas Onuf e François Debrix como arcabouço teórico. Nos guiaremos pela prerrogativa de que a compreensão do mundo e da realidade ocorre de acordo com o modo como nós nos referimos a eles a partir de nossa linguagem. Como referencial metodológico que orientará nosso exame da construção do conceito de inimigo nos discursos de bin Laden, utilizaremos as categorias de operação da ideologia propostas por John Thompson na obra “Ideologia e Cultura Moderna – Teoria Social Crítica na Era dos Meios de Comunicação”. / We aimed to analyze the construction of the concept of enemy represented in Osama bin Laden‟s speeches, during the period 1996 to 2004. These enemies are making from the building of a reality that presents U.S., Israel and its allies, as oppressors of muslims. To better understand this symbolic construction, we will seek to study it based on the dynamics of two pillars for the constitution of bin Laden's politico-religious thought: the theoretical heritage of radical islam, combined with the author's perception about the main events contemporary geopolitical involving the muslim world. From this perspective, we will employ the Linguistic Constructivism of the authors Nicholas Onuf and François Debrix, as theoretical framework. We will be guided by the prerogative of the understanding of the world and reality occurs according to the way we refer to them from our language. As a methodological framework that will guide our examination of the construction of the concept of enemy in bin Laden's speeches, we will use the categories of operation of ideology proposed by John Thompson in his work "Ideology and Modern Culture - Critical Social Theory in the Era of the Media."

Le transport de particules par des écoulements turbulents est un phénomène présent dans de nombreux écoulements naturels et industriels, tels que la dispersion de polluants dans l'atmosphère ou du phytoplancton et plastiques dans et à la surface des océans. Les modèles prédictifs classiques ne peuvent prévoir avec précision la formation de larges fluctuations de concentrations. La première partie de cette thèse concerne une étude de la dispersion turbulente de traceurs émis à partir d'une source ponctuelle et continue. Les fluctuations spatiales de masse sont déterminées en fonction de la distance à la source et à l'échelle d'observation. La combinaison de plusieurs phénomènes physiques à l'origine du mélange limite la validité d'une caractérisation de géométrie fractale. Une approche alternative est proposée, permettant d'interpréter les fluctuations massiques en terme des différents régimes de séparation de pair dans des écoulements turbulents. La seconde partie concerne des particules ayant une inertie finie, dont la dispersion dans l'espace des vitesses requiert de développer des techniques de modélisation adaptées. Une méthode numérique originale est proposée pour exprimer la distribution des particules dans l'espace position-vitesse. Cette méthode est ensuite utilisée pour décrire la modulation de la turbulence bi- dimensionnelle par des particules inertielles. A grand nombres de Stokes, l'effet montré est analogue à celui d'une friction effective à grande échelle. Aux petits Stokes, le spectre de l'énergie cinétique du fluide et les transferts non-linéaires sont modifiées d'une manière non triviale. / The transport of particles by turbulent flows is ubiquitous in nature and industry. It occurs in planet formation, plankton dynamics and combustion in engines. For the dispersion of atmospheric pollutants, traditional predictive models based on eddy diffusivity cannot accurately reproduce high concentration fluctuations, which are of primal importance for ecological and health issues. The first part of this thesis relates to the dispersion by turbulence of tracers continuously emitted from a point source. Mass fluctuations are characterized as a function of the distance from the source and of the observation scale. The combination of various physical mixing processes limits the use of fractal geometric tools. An alternative approach is proposed, allowing to interpret mass fluctuations in terms of the various regimes of pair separation in turbulent flows. The second part concerns particles with a finite and possibly large inertia, whose dispersion in velocity requires developing efficient modelling techniques. A novel numerical method is proposed to express inertial particles distribution in the position-velocity phase space. Its convergence is validated by comparison to Lagrangian measurements. This method is then used to describe the modulation of two-dimensional turbulence by large-Stokes-number heavy particles. At high inertia, the effect is found to be analogous to an effective large-scale friction. At small Stokes numbers, kinetic energy spectrum and nonlinear transfers are shown to be modified in a non-trivial way which relates to the development of instabilities at vortices boundaries.