No description available.

Dry Powder Inhalers (DPIs) have great potential in pulmonary drug delivery; the granular powder, used as active ingredient in DPIs, is ozone friendly and the operation of DPIs ensures coordination between dose release and patient inhalation. However, the powder fluidisation mechanisms are poorly understood which leads to low efficiency of DPIs with 10-35 % of the dose reaching the site of action. The main aim of this thesis is to study the hydrodynamics of powder fluidisation in DPIs, using experimental and computational approaches. An experimental test rig was developed to replicate the process of transient powder fluidisation in an impinging air jet configuration. The powder fluidisation chamber was scaled up resulting in a two dimensional particle flow prototype, which encloses 3.85 mm glass beads. Using optical image processing techniques, individual particles were detected and tracked throughout the experimental time and domain. By varying the air flow rate to the test section, two particle fluidisation regimes were studied. In the first fluidisation regime, the particle bed was fully fluidised in less than 0.25 s due to the strong air jet. Particle velocity vectors showed strong convective flow with no evidence of diffusive motion triggered by inter-particle collisions. In the second fluidisation regime, the particle flow experienced two stages. The first stage showed strong convective flow similar to the first fluidisation regime, while the second stage showed more complex particle flow with collisional and convective flow taking place on the same time and length scales. The continuum Two Fluid Model (TFM) was used to solve the governing equations of the coupled granular and gas phases for the same experimental conditions. Sub-models for particle-gas and particle-particle interactions were used to complete the model description. Inter-particle interactions were resolved using models based on the kinetic theory of granular flow for the rapid flow regime and models based on soil mechanics for the frictional regime. Numerical predictions of the first fluidisation regime showed that the model should incorporate particle-wall friction and minimise diffusion, simultaneously. Ignoring friction resulted in fluidisation timing mismatch, while increasing the diffusion resulted in homogenous particle fluidisation in contrast to the aggregative convective fluidisation noticed in the experiments. Numerical predictions of the second fluidisation regime agreed well with the experiments for the convection dominated first stage of flow up to 0.3 s. However, later stages of complex particle flow showed qualitative discrepancies between the experimental and the computational approaches suggesting that current continuum granular models need further development. The findings of the present thesis have contributed towards better understanding of the mechanics of particle fluidisation and dense multiphase flow in DPI in particular, and particle bed fluidisation using impinging air jet in general. The use of TFM for predicting high speed convective granular flows, such as those in DPIs, is promising. Further studies are needed to investigate the form of particle-particle interactions within continuum granular flow models.

Thesis (PhD (Electrical and Electronic Engineering))--University of Stellenbosch, 2011. / ENGLISH ABSTRACT: The multiphase induction machine drive has been under investigation for the last half century. Although it offers several attractive advantages over the conventional three-phase induction machine drive, it is restricted to highly specialised applications. One aspect of the multiphase induction machine drive is the complexity of the control algorithm for decoupled flux and torque control. The complexity, arising from the required coordinate transformations, increases with increase in the number of phases of the machine. Recently, a method that allows the control of a six-phase induction machine drive without any coordinate transformations was developed and tested. This new control technique allows the control of the machine to be similar to that of dc machines through the use of special trapezoidal-shaped stator current waveforms. These stator phase current waveforms consist of field (flux) and torque current components, with flat-topped amplitudes allowing a stator phase to act alternately in time as either a flux or a torque producing phase. The idea is to have a number of stator phases acting as flux producing phases, whilst the remaining phases act as torque producing phases at each time instance. This dissertation takes a further step in the research on this particular control technique. As the control method relates directly to the brush-dc machine operation, in this dissertation, the control method is defined as a “brush-dc equivalent” (BDCE) control method. First, in this dissertation, a simple analytical method is developed to determine a defined optimal ratio of the number of field to the number of torque phases of a multiphase induction machine that utilises trapezoidal stator current waveforms. The method is applied to induction machines with up to fifteen stator phases. Finite element analysis is used to verify the validity of the developed criterion and to verify the square-like air gap flux density. Secondly, in this dissertation, an analytical method for predicting and evaluating the rotor bar current waveform of a cage multiphase induction machine is proposed. The method is based on the Fourier transform and the winding function theory under linear condition assumptions. The method also allows for the calculation of the electromagnetic torque and rotor bar losses. Skin effect is considered in the calculation of the rotor bar resistance of the machine. Again, finite element analysis is used to verify the analytically calculated results. The developed method can be expanded and used to evaluate the rotor current waveform of any multiphase induction machine supplied with any stator current waveforms. The BDCE control method is implemented on a prototype nine-phase cage-rotor induction machine drive. A nine-phase inverter and control system are developed for supplying the nine-phase induction machine with the trapezoidal stator current waveforms. Rotor current waveform measurements are taken on a specially designed rotor to verify the analytically predicted waveform. The linear relationship of the developed torque and torque current of the proposed BDCE control method is verified through measurements. Through the comparison of analytical calculated results with finite element calculated and measured results, it is shown in this dissertation that the developed analytical techniques can be used in the design and performance analysis of multiphase induction machines. Also, from the results, it is clear that the new control technique works remarkably well even in the flux weakening region. However, outstanding aspects, such as efficiency and generated torque quality of the proposed drive still need to be investigated further. / AFRIKAANSE OPSOMMING: Die multifase induksiemasjien aandryfstelsel word in die laaste halwe eeu al in navorsing ondersoek. Alhoewel dit verskeie aantreklike voordele bied bo die konvensionele driefase induksiemasjien aandryfstelsel, is dit beperk tot hoogs gespesialiseerde aanwendings. Een aspek van die multifase induksiemasjien aandryfstelsel is die kompleksiteit van die beheer algoritme vir ontkoppelde vloed en draaimoment beheer. Die kompleksiteit, wat voortspruit uit die vereiste koördinaat transformasies, neem toe met toename in die aantal fases van die masjien. Onlangs is 'n metode wat die beheer van ’n sesfase induksiemasjien sonder enige koördinaat transformasies doen, ontwikkel en getoets. Hierdie nuwe beheertegniek maak die beheer van die masjien soortgelyk aan dié van GS masjiene deur die gebruik van spesiale trapezium-vormige statorstroom golfvorms. Hierdie stator fasestroom golfvorms bestaan uit veld- (vloed-) en draaimoment-stroom komponente met plat amplitudes, sodat 'n statorfase om die beurt in tyd optree as óf' ’n vloed of 'n draaimoment genereerde fase. Die idee is om 'n aantal statorfases te hê wat as vloed genereerde fases dien, terwyl die oorblywende fases as draaimoment genereerde fases optree op enige tydstip. Hierdie tesis neem 'n verdere stap in die navorsing op hierdie spesifieke beheertegniek. Met die beheermetode wat direk verband hou met borsel-GS masjien werking, word in hierdie proefskrif die beheermetode as 'n "borsel-GS ekwivalente" ["brush-DC equivalent" (BDCE)] beheermetode gedefinieer. In die eerste plek word in hierdie proefskrif 'n eenvoudige analitiese metode ontwikkel om ’n gedefinieerde optimale verhouding van die aantal veld tot die aantal draaimoment fases van 'n multifase induksiemasjien te bepaal, wat van trapesoïdale statorstroom golfvorms gebruik maak. Die metode word toegepas op induksiemasjiene met tot vyftien statorfases. Eindige element analise is gebruik om die geldigheid van die ontwikkelde kriterium te verifieer en om die vierkantvormige luggaping vloeddigtheid te verifieer. In die tweede plek word in hierdie proefskrif 'n analitiese metode vir die voorspelling en evaluering van die rotorstaafstroom golfvorm van 'n kourotor multifase induksiemasjien voorgestel. Die metode is gebaseer op die Fourier transform en die wikkelingsfunksie teorie onder lineêre-toestand aannames. Die metode wend hom ook daartoe tot die berekening van die elektromagnetiese draaimoment en rotorstaafverliese. Die huideffek word in ag geneem in die berekening van die rotorstaafweerstand van die masjien. Weereens is eindige element analise gebruik om die analitiese berekende resultate te verifieer. Die ontwikkelde metode kan uitgebrei en gebruik word om die rotorstroom golfvorm van van enige multifase induksiemasjien te evalueer wat gevoer word met enige statorstroom golfvorms. Die BDCE beheermetode is toegepas op 'n prototipe negefase kourotor induksiemasjien. 'n Negefase omsetter en beheerstelsel is ontwikkel vir die toevoer van die trapesoïdale statorstroom golfvorms aan die negefase induksiemasjien. Die rotorstroomgolfvorm metings is geneem op 'n spesiaal ontwerpte rotor om die analitiese voorspelde golfvorm te verifieer. Die lineêre verwantskap tussen die ontwikkelde draaimoment en draaimomentstroom van die voorgestelde BDCE beheermetode is geverifieer deur metings. Deur die analitiese berekende resultate met die eindige element berekende en gemete resultate te vergelyk, wys hierdie proefskrif dat die ontwikkelde analitiese tegnieke gebruik kan word in die ontwerp en werkverrigting analise van ’n multifase induksiemasjien. Vanuit die resultate is dit ook duidelik dat die nuwe beheertegniek besonder goed werk, selfs in die vloedverswakking spoedgebied. Egter, uitstaande aspekte soos effektiwiteit en genereerde draaimoment kwaliteit van die voorgestelde aandryfstelsel moet nog verder ondersoek word.

Thesis (PhD)--Stellenbosch University, 2012. / ENGLISH ABSTRACT: See full text for abstract / AFRIKAANSE OPSOMMING: Sien volteks vir opsomming

The dynamic solution of multiphase flow through porous media is of special interest to several fields of science and engineering, such as petroleum, geology and geophysics, bio-medical, civil and environmental, chemical engineering and many other disciplines. A natural application is the modeling of the flow of two immiscible fluids (phases) in a reservoir. Others, that are broadly based and considered in this work include the hydrodynamic dispersion (as in reactive transport) of a solute or tracer chemical through a fluid phase. Reservoir properties like permeability and porosity greatly influence the flow of these phases. Often, these vary across several orders of magnitude and can be discontinuous functions. Furthermore, they are generally not known to a desired level of accuracy or detail and special inverse problems need to be solved in order to obtain their estimates. Based on the physics dominating a given sub-region of the porous medium, numerical solutions to such flow problems may require different discretization schemes or different governing equations in adjacent regions. The need to couple solutions to such schemes gives rise to challenging domain decomposition problems. Finally, on an application level, present day environment concerns have resulted in a widespread increase in CO₂capture and storage experiments across the globe. This presents a huge modeling challenge for the future. This research work is divided into sections that aim to study various inter-connected problems that are of significance in sub-surface porous media applications. The first section studies an application of mortar (as well as nonmortar, i.e., enhanced velocity) mixed finite element methods (MMFEM and EV-MFEM) to problems in porous media flow. The mortar spaces are first used to develop a multiscale approach for parabolic problems in porous media applications. The implementation of the mortar mixed method is presented for two-phase immiscible flow and some a priori error estimates are then derived for the case of slightly compressible single-phase Darcy flow. Following this, the problem of modeling flow coupled to reactive transport is studied. Applications of such problems include modeling bio-remediation of oil spills and other subsurface hazardous wastes, angiogenesis in the transition of tumors from a dormant to a malignant state, contaminant transport in groundwater flow and acid injection around well bores to increase the permeability of the surrounding rock. Several numerical results are presented that demonstrate the efficiency of the method when compared to traditional approaches. The section following this examines (non-mortar) enhanced velocity finite element methods for solving multiphase flow coupled to species transport on non-matching multiblock grids. The results from this section indicate that this is the recommended method of choice for such problems. Next, a mortar finite element method is formulated and implemented that extends the scope of the classical mortar mixed finite element method developed by Arbogast et al [12] for elliptic problems and Girault et al [62] for coupling different numerical discretization schemes. Some significant areas of application include the coupling of pore-scale network models with the classical continuum models for steady single-phase Darcy flow as well as the coupling of different numerical methods such as discontinuous Galerkin and mixed finite element methods in different sub-domains for the case of single phase flow [21, 109]. These hold promise for applications where a high level of detail and accuracy is desired in one part of the domain (often associated with very small length scales as in pore-scale network models) and a much lower level of detail at other parts of the domain (at much larger length scales). Examples include modeling of the flow around well bores or through faulted reservoirs. The next section presents a parallel stochastic approximation method [68, 76] applied to inverse modeling and gives several promising results that address the problem of uncertainty associated with the parameters governing multiphase flow partial differential equations. For example, medium properties such as absolute permeability and porosity greatly influence the flow behavior, but are rarely known to even a reasonable level of accuracy and are very often upscaled to large areas or volumes based on seismic measurements at discrete points. The results in this section show that by using a few measurements of the primary unknowns in multiphase flow such as fluid pressures and concentrations as well as well-log data, one can define an objective function of the medium properties to be determined, which is then minimized to determine the properties using (as in this case) a stochastic analog of Newton’s method. The last section is devoted to a significant and current application area. It presents a parallel and efficient iteratively coupled implicit pressure, explicit concentration formulation (IMPEC) [52–54] for non-isothermal compositional flow problems. The goal is to perform predictive modeling simulations for CO₂sequestration experiments. While the sections presented in this work cover a broad range of topics they are actually tied to each other and serve to achieve the unifying, ultimate goal of developing a complete and robust reservoir simulator. The major results of this work, particularly in the application of MMFEM and EV-MFEM to multiphysics couplings of multiphase flow and transport as well as in the modeling of EOS non-isothermal compositional flow applied to CO₂sequestration, suggest that multiblock/multimodel methods applied in a robust parallel computational framework is invaluable when attempting to solve problems as described in Chapter 7. As an example, one may consider a closed loop control system for managing oil production or CO₂sequestration experiments in huge formations (the “instrumented oil field”). Most of the computationally costly activity occurs around a few wells. Thus one has to be able to seamlessly connect the above components while running many forward simulations on parallel clusters in a multiblock and multimodel setting where most domains employ an isothermal single-phase flow model except a few around well bores that employ, say, a non-isothermal compositional model. Simultaneously, cheap and efficient stochastic methods as in Chapter 8, may be used to generate history matches of well and/or sensor-measured solution data, to arrive at better estimates of the medium properties on the fly. This is obviously beyond the scope of the current work but represents the over-arching goal of this research. / text

The knowledge of the behavior of non-spherical particles suspended in turbulent flows covers a wide range of applications in engineering and science. Dispersed two-phase flows and turbulence are the most challenging subjects in engineering, and when combined it gives rise to more complexities as the result of the inherent stochastic nature of the turbulence of the carrier-phase together with the random distribution of the dispersed phase. Moreover, for anisotropic particles the coupling between the translation and rotation of particle increases the complication. Because of the practical importance of prolate particleladen turbulent flows, the plenty of numerical and experimental works have been conducted to study such suspensions. Numerical approaches have given valuable insight of turbulent suspension flows, although the computation has been only carried out at the macro scale and models, not including flow distortion around the particle, comprise the detail of the flow in the order of a particle size. In addition, the model of the forces imposed on the particle by the fluid and mass point treatment are strictly valid for infinitely small particle having size less than all scales of the fluid turbulence. Fully resolved solution at the scale of the dispersed phase in turbulent flows for high Reynolds number has been recently performed but is still a challenge. On the other hand, the presence of particle as the dispersed phase makes experimental measurements much more complicated than those with single phase as a result of particles interference. The area of considerable difficulty with this type of experiments is the measurement of the fluid-phase velocity remarkably close to the particle surface. Generally, experimental researches have been concentrated on measuring the mean velocity and Reynolds stresses of the carrier-phase, and the mean velocity, fluctuations, orientation and accumulation of the non-spherical particles. Higher-order quantities, including Lagrangian particle velocity correlations, the carrier-phase turbulence modulation, and two-particle and particlefluid velocity correlations are also of interest. It has been found that the rotational and translational movements of the fibershaped particle depend on the nature of carrier-phase field and fiber characteristics such as aspect ratio, fiber Stokes number, fiber Reynolds number, and the ratio of fiber to flow length scale. With the development of PIV (Particle Image Velocimetry) and PTV (Particle Tracking Velocimetry) techniques, it has been appeared that combined PIV/PTV will be the best available choice for the experimental study of dispersed two-phase flows. The purpose of combined PIV/PTV measurement of two-phase systems is simultaneous measurements of fluid and suspended objects, where the PIV measurement of the fluid phase are combined with PTV measurement of the dispersed phase. The objective of this doctoral thesis is to study the behavior of rod-like particles suspended in wall-bounded turbulent flow through simultaneous PIV/PTV measurements of the velocity of the flow field and particle motion. As a representative of rod-like particles, I have employed cellulose acetate fibers with the length to diameter ratio (aspect ratio) larger than one. Here, It has been considered only dilute suspensions with no flocculation; thus fiber-fiber interaction is negligible. The measurements have been conducted within the parallel planes (2D view) illuminated by laser in the streamwise direction in thin film suspension flowing on the water table setup at Linné FLOW Centre, KTH Mechanics Lab. It is shown that this setup is a well-behaved experimental model of half channel flows often used in Direct Numerical Simulation (DNS) investigations. Therefore, the experimental results are comparable to their DNS counterpart where it is convenient. A single camera PIV technique has been used to measure flowing suspension. Therefore, it has been needed to preprocess images using a spatial median filter to separate images of two phases, tracer particles as representative of fluid and fibers suspended. The well-known PIV processing algorithms have been applied to the phase of fluid. I have also introduced a novel algorithm to recognize and match fibers in consecutive images to track fibers and estimate their velocity. It is not feasible to study all relevant aspects of particle-laden turbulent flows in a single study. In this study, I present the statistics of the rotational and translational motion of fiber-like particles and the surrounding fluid velocity. To the author’s knowledge, remarkably little experimental work has been published to date on simultaneous measurement of fiber motion and turbulence field in a turbulent fiber suspension flow to reveal dynamics of fibers in this regime. Therefore, the results of this work will be profitable in better understanding of such multiphase flows. The statistical analysis of the translational motion of fibers shows that the size of fiber is a significant factor for the dynamical behavior of the fiber near the wall. It has been observed that, in the region near the wall, the probability of presence of the long fibers is high in both the high-speed and low-speed streaks of flow, and the mean velocity of fibers almost conforms to the mean velocity of flow; whereas the short fibers are mostly present in the low-speed areas, and the fiber mean velocity obey the dominant flow velocity in these areas. In the far-wall regions, the translation of fibers is practically unaffected by the aspect ratio, whereas it depends crucially on the wall-normal distance. Moreover, it was found that in the case of long fibers near the wall, the low speed fibers mostly are orientated in streamwise direction. On the other hand, there is no preferential orientation for fast long fibers. Although wall-normal velocities were not measured in this study, it is hypothesized that this behavior is a result of fibers being affected by the sweep and ejection events known to occur in wall-bounded turbulent flow. The fast fibers are in sweep environment and comes from the upper layer. The low speed fibers are into ejection areas in the vicinity of the wall, and the wall has a stabilizing effect on them. The short fibers are still oriented mostly in streamwise direction for a certain range of low velocity. Furthermore, since a considerable change of the fiber behavior is observed in a certain ratio of the fiber length to the fiber distance from the solid wall, it is supposed that this ratio is also a prominent parameter for the behavior of fiber near the wall. The results presented are in terms of viscous wall units wherever are denoted by superscript “+”.

No description available.

Ce travail est consacré à l'étude expérimentale des écoulements diphasiques et triphasiques d'inclusions (bulles, gouttes, sphères) en milieux tant newtoniens que non newtoniens à l'échelle microscopique et mésoscopique, en utilisant la visualisation par une caméra rapide, la vélocimétrie par images des particules (PIV) ainsi que la micro-vélocimétrie par images des particules. Des bulles et des gouttes ont été étudiées expérimentalement depuis leur formation, en passant par leur déformation jusqu'à leur coalescence. La formation de bulles dans des micro-mélangeurs a été étudiée et caractérisée par l'obtention de champs de vitesses. Différents paramètres, tels que le cisaillement, la géométrie de la zone de formation, les débits ou encore les propriétés physiques ont été testés afin de développer des lois d'échelles. La traversée d'une interface liquide-liquide par une inclusion a été abordée par des expériences originales, permettant de décrire la dynamique du phénomène, de définir des nombres adimensionnels et de mettre en évidence des instabilités interfaciales. L'effet Weissenberg a également été étudié aux différentes échelles afin de comprendre les phénomènes conduisant à son amplification. Enfin, en milieu viscoélastique et rhéofluidifiant, nous avons caractérisé l'écoulement autour d'une inclusion isolée solide par l'obtention de champs de vitesses. Ces résultats ont permis de confirmer l'origine viscoélastique du sillage négatif et de prédire ses caractéristiques / The present work was devoted to the experimental study of the multiphase flow around inclusions in both Newtonian and non-Newtonian media at respectively microscopic and mesoscopic scales, by means of the Particle Image Velocimetry (both PIV and µ-PIV) and fast camera visualization. Bubbles and drops were experimentally studied starting from their formation, the rising and up to their coalescence and fragmentation. Bubble formation in micro-mixers was also investigated and characterized by measuring the liquid velocity fields. Different parameters, such as the shear rate, the geometry of the micro-mixer, the flow rates or the physical properties were tested to develop correlations of power-law kind. The deformation of a liquid-liquid interface due to the passage of an inclusion was investigated by original experiments in order to observe and describe the dynamics of the phenomenon, to define the reliable dimensionless numbers and also to highlight several interfacial instabilities. The Weissenberg effect was also studied at different scales in various non-Newtonian fluids in comparison with Newtonian fluids to understand the amplification phenomena under the combining effects of instability and normal forces. Finally, the flow around an isolated solid inclusion was characterized by performing the measurements of velocity fields in viscoelastic and shear-thinning fluids. These results were used to confirm the viscoelastic origin of the negative wake behind the solid sphere and to model its characteristics

This thesis reports on investigations carried out to study of the effect of horizontal wall jets on rough, fixed and mobile beds in open channel flow. Experimental tests were carried out, using fixed and mobile sediment beds. Computer simulation models for the flow within the jet and resulting sediment transport were developed and their results analysed in this study. In the experimental phase, tests were carried out with both fixed and mobile sediment beds. The shape of the water surface, numerous point velocity measurements and measurements of the evolving scour hole shape were made. Detailed descriptions of the turbulent flow field over a fixed rough bed and for scour holes at equilibrium were obtained for a range of initial jet conditions. Fully turbulent, multiphase flow was modelled using the Fluent Computational Fluid Dynamics software. This was used to analyze the flow caused by a jet in a rectangle open-channel with a rough bed, and also the flow pattern in a channel with a local scour hole. The volume of fluid (VOF) multiphase method and K- model was used to model the fluid flow in both cases. The model predictions of velocity and shear stress were compared against experimental observations. The experimental data was used to develop new empirical relationships to describe the pattern of boundary shear stress caused by a wall jet over fixed beds and in equilibrium scour holes. These relationships were linked with existing bed-load transport rate models in order to predict the temporal evolution of scour holes. An analytical model describing the relationship between the wall jet flow and the development of a local scour hole shape was reported and its predictions compared with experimental data.

Ce manuscrit s'intéresse au devenir de contaminants organiques dans le sol, et plus précisément celui des composés organo-halogénés volatils (COHV).Il propose des outils d'aide à la décision en utilisant la modélisation numérique appliquée à des problématiques environnementales portant sur le traitement de sites et sols pollués par des COHV. Il présente, à ma connaissance, la première modélisation numérique en 3D, grâce au simulateur TMVOC, d'une technique de dépollution physique, à savoir le pompage réussi au sein d'un aquifère d'un liquide en phase non-aqueuse plus dense que l'eau ou Dense Non-Aqueous Phase Liquid (DNAPL). Les très bons résultats de cette simulation permettent d’envisager l’optimisation d’un système de pompage asservi pour dépolluer un site contaminé aux COHV. Cette thèse s'intéresse aussi à une méthode de d'évaluation, à la fois qualitative et quantitative, de l'efficacité du pompage : un test de traçage utilisant des traceurs bisolubles à coefficients de partage variables (partitioning interwell tracer test – PITT). Ce PITT permet de connaître à la fois la répartition spatiale, au sein d’un aquifère, d’une bulle de DNAPL et aussi d’en évaluer sa saturation et donc son volume. Ce manuscrit met à disposition des outils élaborés et adaptés au traitement d'un DNAPL dans un aquifère. La combinaison des deux techniques présentées, à savoir la modélisation d'un pompage de DNAPL et un PITT, sont parfaitement reproductibles dans des environnements similaires et à l'échelle industrielle. Enfin, ces méthodes permettent de réduire considérablement les coûts de caractérisation (PITT) et d'exploitation par l'optimisation de systèmes de pompage / This manuscript deals with the fate of organic contaminants in soil, more precisely of volatile organo-chlorinated compounds (VOHC) and offers some decision making techniques and tools using numerical modelling applied to environmental issues about the treatment of soils contaminated by VOHC. It presents, to the best knowledge of the auhor, the first 3D numerical modelling, with the simulator TMVOC, of a physical treatment technique, namely the successful pumping within an aquifer, of a Dense Non-Aqueous Phase Liquid (DNAPL). The very good results of this simulation give the opportunity to design a controlled and automated pumping system to remediate a polluted site. This thesis also deols with a tracer test assessment method, both qualitative and quantitative, for the efficiency of this pumping : a partitioning interwell tracer test (PITT). The PITT allows us to determine the spatial repartion of the DNAPL and also to assess its saturation, hence its volume. This manuscript offers elaborated tools adapted to the remediation of a DNAPL in an aquifer. The combination of these two techniques, namely the DNAPL pumping and the PITT, are perfectly reproducible in similar environments up to an industrial scale. Finally, exploitation and characterisation costs for DNAPL remediation can be extremely reduced by numerical modelling and optimisation