Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Um código computacional para escoamentos bifásicos incorporando metodologia híbrida entre oMétodo dos Elementos Finitos e a descrição Lagrangeana-Euleriana Arbitrária do movimento é usado para simular a dinâmica de um jato transversal de gotas na zona primária de quebra. Os corpos dispersos são descritos por meio de um método do tipo front-tracking que produz interfaces de espessura zero através de malhas formadas pela união de elementos adjacentes em ambas as fases e de técnicas de refinamento adaptativo. Condições de contorno periódicas são implementadas de modo variacionalmente consistente para todos os campos envolvidos nas simulações apresentadas e uma versão modificada do campo de pressão é adicionada à formulação do tipo um-fluido usada na equação da quantidade de movimento linear. Simulações numéricas diretas em três dimensões são executadas para diferentes configurações de líquidos imiscí veis compatíveis com resultados experimentais encontrados na literatura. Análises da hidrodinâmica do jato transversal de gotas nessas configurações considerando trajetórias, variação de formato de gota, espectro de pequenas perturbações, além de aspectos complementares relativos à qualidade de malha são apresentados e discutidos. / A two-phase flow computational code taking a hybrid Arbitrary Lagrangian-Eulerian desciption of movement along with the Finite Element Method is used to simulate the dynamics of an incompressible drop jet in crossflow in the primary breakup zone. Dispersed entities are described by means of a front-tracking method which produces zero-thickness interfaces through contiguous element meshing and adaptive refinement techniques. Periodic boundary conditions are implemented in a variationally consistent way for all the scalar fields involved in the presented simulations and amodified version of the pressure field is added to the one-fluid formulation employed in the momentum equation. Three-dimensional direct numerical simulations for different flow configurations of immiscible liquids pertinent to experimental results found in literature. Analyses of the hydrodynamics of the drop jet in crossflow in these configurations considering trajectories, drop shape variations, spectrum of small disturbances, besides additional aspects relating to mesh quality are presented and discussed.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Um código computacional para escoamentos bifásicos incorporando metodologia híbrida entre oMétodo dos Elementos Finitos e a descrição Lagrangeana-Euleriana Arbitrária do movimento é usado para simular a dinâmica de um jato transversal de gotas na zona primária de quebra. Os corpos dispersos são descritos por meio de um método do tipo front-tracking que produz interfaces de espessura zero através de malhas formadas pela união de elementos adjacentes em ambas as fases e de técnicas de refinamento adaptativo. Condições de contorno periódicas são implementadas de modo variacionalmente consistente para todos os campos envolvidos nas simulações apresentadas e uma versão modificada do campo de pressão é adicionada à formulação do tipo um-fluido usada na equação da quantidade de movimento linear. Simulações numéricas diretas em três dimensões são executadas para diferentes configurações de líquidos imiscí veis compatíveis com resultados experimentais encontrados na literatura. Análises da hidrodinâmica do jato transversal de gotas nessas configurações considerando trajetórias, variação de formato de gota, espectro de pequenas perturbações, além de aspectos complementares relativos à qualidade de malha são apresentados e discutidos. / A two-phase flow computational code taking a hybrid Arbitrary Lagrangian-Eulerian desciption of movement along with the Finite Element Method is used to simulate the dynamics of an incompressible drop jet in crossflow in the primary breakup zone. Dispersed entities are described by means of a front-tracking method which produces zero-thickness interfaces through contiguous element meshing and adaptive refinement techniques. Periodic boundary conditions are implemented in a variationally consistent way for all the scalar fields involved in the presented simulations and amodified version of the pressure field is added to the one-fluid formulation employed in the momentum equation. Three-dimensional direct numerical simulations for different flow configurations of immiscible liquids pertinent to experimental results found in literature. Analyses of the hydrodynamics of the drop jet in crossflow in these configurations considering trajectories, drop shape variations, spectrum of small disturbances, besides additional aspects relating to mesh quality are presented and discussed.

Numerical solution of differential equations having multitude of scales in the solution field is one of the most challenging research areas, but highly demanded in scientific and industrial applications. One of the natural approaches for handling such problems is to separate the scales and approximate the solution of the segregated scales with appropriate numerical method. Variational multiscale method (VMS) is a predominant method in the paradigm of scale separation schemes. In our work we have used the VMS technique to develop a numerical scheme for computations of turbulent flows in time-dependent domains. VMS allows separation of the entire range of scales in the flow field into two or three groups, thereby enabling a different numerical treatment for the different groups. In the context of computational fluid dynamics(CFD), VMS is a significant new improvement over the classical large eddy simulation (LES). VMS does away with the commutation errors arising due to filtering in LES. Further, in a three-scale VMS approach the model for the subgrid scale can be contained to only a part of the resolved scales instead of effecting the entire range of resolved scales. The projection based VMS scheme that we have developed gives a robust and efficient method for solving problems of turbulent fluid flows in deforming domains, governed by incompressible Navier {Stokes equations. In addition to the existing challenges due to turbulence, the computational complexity of the problem increases further when the considered domain is time-dependent. In this work, we have used an arbitrary Lagrangian-Eulerian (ALE) based VMS scheme to account for the domain deformation. In the proposed scheme, the large scales are represented by an additional tensor valued space. The resolved large and small scales are computed in a single unified equation, and the effect of unresolved scales is confined only to the resolved small scales, by using a projection operator. The popular Smagorinsky eddy viscosity model is used to approximate the effects of unresolved scales. The used ALE approach consists of an elastic mesh update technique. Moreover, a computationally efficient scheme is obtained by the choice of orthogonal finite element basis function for the resolved large scales, which allows to reformulate the ALE-VMS system matrix into the standard form of the NSE system matrix. Thus, any existing Navier{Stokes solver can be utilized for this scheme, with modifications. Further, the stability and error estimates of the scheme using a linear model of the NSE are also derived. Finally, the proposed scheme has been validated by a number of numerical examples over a wide range of problems.

Este trabalho consiste no desenvolvimento de ferramentas computacionais para análise não linear geométrica de interação fluido-casca utilizando o Método dos Elementos Finitos (MEF). O algoritmo para dinâmica dos fluidos é explícito e a integração temporal é baseada em linhas características. O código computacional é capaz de simular as equações de Navier-Stokes para escoamentos compressíveis tanto na descrição Euleriana como na descrição Lagrangeana-Euleriana arbitrária (ALE), na qual é possível prescrever movimentos para a malha do fluido. A estrutura é modelada em descrição Lagrangeana total através de uma formulação de MEF para análise dinâmica não linear geométrica de cascas baseada no teorema da mínima energia potencial total escrito em função das posições nodais e vetores generalizados e não em deslocamentos e rotações. Essa característica evita o uso de aproximações de grandes rotações. Dois modelos de acoplamentos são desenvolvidos. O primeiro modelo, ideal para problemas onde a escala de deslocamentos não é muito grande comparada com as dimensões do domínio do fluido, é baseado na descrição ALE e o acoplamento entre as duas diferentes malhas é feito através do mapeamento das posições locais dos nós do contorno do fluido sobre os elementos de casca e vice-versa, evitando a necessidade de coincidência entre os nós da casca e do fluido. A malha do fluido é adaptada dinamicamente usando um procedimento simples baseado nas posições e velocidades nodais da casca. O segundo modelo de acoplamento, ideal para problemas com grande escala de deslocamentos tais como estruturas infláveis, considera a casca imersa na malha do fluido e consiste em um procedimento robusto baseado em curvas de nível da função distância assinalada do contorno, o qual integra o algoritmo Lagrangeano de casca com o Fluido em descrição Euleriana, sem necessidade de movimentação da malha do fluido, onde a representação computacional do fluido se resume a uma malha não estruturada maior ou igual ao domínio inicial do fluido e a interface fluido-casca dentro da malha do fluido é identificada por meio de curvas de nível da função distância assinalada do contorno. Ambos os modelos são testados através de exemplos numéricos mostrando robustez e eficiência. Finalmente, como uma sugestão para o futuro desenvolvimento desta pesquisa, iniciaram-se estudos relativos a funções B-splines. O uso desse tipo de funções deverá resolver problemas de estabilidade relativos a oscilações espúrias devidas ao uso de polinômios de Lagrange para a representação de descontinuidades. / This work consists of the development of computational tools for nonlinear geometric fluid-shell interaction analysis using the Finite Element Method (FEM). The fluid solver is explicit and its time integration based on characteristics. The computational code is able to simulate the Navier-Stokes equations for compressible flows written in the Eulerian description as well as in the arbitrary Lagrangian-Eulerian (ALE) description, enabling movements prescription for the fluid mesh. The structure is modeled in a total Lagrangian description, using a FEM formulation to deal with geometrical nonlinear dynamics of shells based on the minimum potential energy theorem written regarding nodal positions and generalized unconstrained vectors, not displacements and rotations, avoiding the use of large rotation approximations. Two partitioned coupling models are developed. The first model, ideal for simulations where the displacements scale is not very large compared to the fluid domain, is based on the ALE description and the coupling between the two different meshes is done by mapping the fluid boundary nodes local positions over the shell elements and vice-versa, avoiding the need for matching fluid and shell nodes. The fluid mesh is adapted using a simple approach based on shell nodal positions and velocities. The second model, ideal for problems with large scales of displacements such as inflatable structures, is based on immersed boundary and consists of a robust level-set based approach that integrates the Lagrangian shell finite and the Eulerian finite element high speed fluid flow solver, with no need for mesh adaptation, where the fluid representation relies on a fixed unstructured mesh larger or equal to the initial fluid domain and the fluid-shell interface inside the fluid mesh is tracked with level sets of a boundary signed distance function. Both models are tested with numerical examples, showing efficiency and robustness. Finally, as a suggestion for future development of this research, we started studies relatives to B-Spline functions. The use of this kind of functions should solve stability problems related to spurious oscillations due to the use of Lagrange polynomials for representing discontinuities.

A interação entre fluidos e estruturas caracteriza um problema multi-físico não linear e está presente numa grande variedade de áreas da engenharia. Este trabalho apresenta o desenvolvi mento de ferramentas computacionais com base no Método dos Elementos Finitos (MEF) para a análise de interação fluido-estrutura (IFE) considerando escoamentos com baixas velocidades. Dada a interdisciplinaridade do tema, se faz necessário o estudo em três diferentes assuntos: a dinâmica das estruturas computacional, a dinâmica dos fluidos computacional, e o problema de acoplamento. No caso da dinâmica das estruturas empregar-se um elemento finito que seja adequado para a simulação de problemas de IFE, que claramente demandam uma análise não linear geométrica, optando-se pelo emprego de uma formulação descrita em posições, a qual evita problemas relativos à aproximação de rotações finitas. Quanto à dinâmica dos fluidos computacional, é empregado um método estável e ao mesmo tempo sensível à movimentação da estrutura, utilizando a descrição Lagrangeana-Euleriana Arbitrária (ALE). Os casos considerados neste trabalho, assim como muitos dos problemas de engenharia, ocorrem com escoamentos em baixas velocidades, implicando na incompressibilidade do fluido, o que demanda, para um método estável, a utilização de elementos que atendam à condição de Ladyzhenskaya-Babuska-Brezzi (LBB). Além disso, é necessário também o emprego de métodos que consigam neutralizar as variações espúrias decorrentes da não-linearidade de possíveis escoamentos com convecção dominante e que surgem com a aplicação do processo clássico de Galerkin. Para superar esse problema, é aplicado o método Streamline-Upwind/Petrov-Galerkin (SUPG), que adiciona difusividade artificial na direção do escoamento, controlando a amplitude dos termos convectivos. No que se refere ao acoplamento fluido-casca, buscam-se modularidade e versatilidade adotando-se o modelo particionado. O modelo de acoplamento implementado garante ainda a utilização de malhas do fluido e da estrutura sem a necessidade de coincidência de nós. / Interaction between fluids and structures characterizes a nonlinear multi-physics problem presente in a wide range of engineering fields. This works presets the development of computational tools based on finite element method (FEM) for fluid-structure interaction (FSI) analysis considering low speed flows (incompressible), as a great part of the engineering problems. Given the topic multidisciplinary nature, it is necessary to study three different subjects: the computational structural dynamics, the computational fluid mechanics and the coupling problem. Regarding structural mechanics, we seek to employ a finite element adequate to FSI simulation, what clearly demands a geometric nonlinear analysis. We chose to employ shell elements with formulation in terms of positions, which avoids problems related to finite rotations approximations. Concerning computational fluid dynamics, we employ a stable method, at same time sensible o structural movements, which is written in the arbitrary Lagrangian-Eulerian (ALE) description. The flow incompressibility demands, for a stable method, the use of elements according to the Ladyzhenskaya-Bbuska-Brezzi (LBB) condition. It is also necessary to employ methods able to neutralize the spurious variations that appears from convection dominated flows when applying the standard Galerking method. In order to overcome this problem, we apply the Streamline-Upwind/Petrov-Galerkin (SUPG) method, which adds artificial diffusivity to the streamline direction, controlling spurious variations. Considering the fluid-shell coupling, we seek modularity and versatility, adopting the partitioned model. The developed coupling model ensure the use of fluid and structure meshes with no need for matching nodes.

Um método lagrangeano-euleriano arbitrário para a resolução de escoamentos dominados por tensão superficial é apresentado neste trabalho. Tais escoamentos são importantes em muitas aplicações, especialmente em canais capilares que frequentemente aparecem em escoamentos em microescala. A resolução deste tipo de escoamento apresenta vários desafios que são abordados neste trabalho. O escoamento é resolvido somente para a fase líquida, com condições de contorno apropriadas para a superfície livre que delimita o líquido e o gás, que é representada por arestas e vértices da malha computacional. Esta se move e se deforma, sendo que sua qualidade é mantida sob controle para não degradar a solução numérica. As equações de Navier-Stokes são discretizadas pelo método de elementos finitos em um referencial arbitrário. O método de incorporação dos efeitos de tensão superficial e linha de contato é explicado em detalhes. Validações comprovam a precisão do método proposto, com comparações através de soluções pseudo-analíticas para casos simples. Finalmente alguns resultados sobre escoamentos em capilares são apresentados / An arbitrary lagrangian-eulerian finite element method to solve surface tension dominated flows is presented. Such flows are important in many applications, particularly in capillary channels, that appear in microscale flows. The resolution of such flows presents several challenges that are addressed in this work. The flow is solved only in the liquid phase, and proper boundary conditions are applied on the free-surface, bounding the liquid and gas, which is explicitly represented by vertices and edges of the computational mesh. The mesh is moved and deformed, but its quality is kept under control in order to control errors in the numerical solution. The Navier-Stokes equations are discretized by standard Galerkin finite element method in an arbitrary reference. Details of the computation of surface tension and contact line effects are presented. The methodology is validated for a number of simple test cases against known pseudo-analytical solutions, and numerical results are presented, showing the robustness and accuracy of the methodology. Finally, some results about surface-tension-driven flows in capillaries are presented

Submitted by Geandra Rodrigues (geandrar@gmail.com) on 2018-04-18T13:51:58Z No. of bitstreams: 1 marcusviniciusdesouzaferraz.pdf: 4103901 bytes, checksum: e4adcd64380c6ba8941b29bcc9d0abfd (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2018-04-19T17:48:41Z (GMT) No. of bitstreams: 1 marcusviniciusdesouzaferraz.pdf: 4103901 bytes, checksum: e4adcd64380c6ba8941b29bcc9d0abfd (MD5) / Made available in DSpace on 2018-04-19T17:48:41Z (GMT). No. of bitstreams: 1 marcusviniciusdesouzaferraz.pdf: 4103901 bytes, checksum: e4adcd64380c6ba8941b29bcc9d0abfd (MD5) Previous issue date: 2018-02-27 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O conhecimento da topografia das superfícies e uma compreensão da interação entre elas é essencial para qualquer estudo que envolva os fenômenos de atrito, desgaste e lubrificação. O estudo da relação entre o atrito e os parâmetros de rugosidade é um problema difícil e de interesse tanto industrial como acadêmico e trabalhos experimentais e teóricos têm mostrado que uma película de fluido entre duas superfícies rugosas em movimento relativo impede o contato sólido - sólido e pode proporcionar atrito muito baixo e desgaste desprezível. A modelagem matemática utilizada neste trabalho é baseada em modelos clássicos, tais como a equação de Reynolds para a descrição dos fenômenos hidrodinâmicos e as formulações de Hertz (1896) e Greenwood e Williamson (1966) para a modelagem do contato das asperezas entre as superfícies rugosas. Para tratar a complexidade das interações entre o fluido e os pares sólidos contactados, a descrição Lagrangiana-Euleriana Arbitrária é apresentada nesta pesquisa. Através do Método dos Elementos Finitos um modelo tridimensional é gerado no Abaqus ®, a fim de identificar as pressões de contato, as tensões tangenciais e normais resultantes e os coeficientes de atrito decorrrentes do deslizamento entre uma superfície texturizada e lubrificada e um plano rígido (em analogia aos modelos de contato clássicos), cujos perfis de rugosidade são construídos a partir de informações da rugosidade média quadrática de superfícies dentárias. São avaliados também a sensibilidade de alguns parâmetros do lubrificante na determinação do coeficiente de atrito e são propostos modelos com condições de contorno distintas. Entretanto, para a verificação destes últimos busca-se reproduzir qualitativamente o resultado encontrado por Lorentz (2013) na investigação numérica de sistemas tribológicos no regime misto de lubrificação. A metodologia aqui proposta emerge como uma alternativa eficaz no campo da Tribologia, na predição do coeficiente de atrito e outras variáveis pertinentes a um fenômeno ainda pouco compreendido. Realiza-se uma análise de sensibilidade dos parâmetros de modelagem, a fim de identificar como os mesmos afetam consideravelmente o comportamento mecânico na interface de contato. / The knowledge of the topography of surfaces and an understanding of the interaction between them is essential for any study involving the phenomena of friction, wear and lubrication. The study of the relationship between friction and roughness parameters is a difficult problem of both industrial and academic interest and experimental and theoretical works have shown that a fluid film between two rough surfaces in relative motion prevents solid - solid contact and can provide very low friction and negligible wear. The mathematical modeling used in this paper is based on classical models, such as the Reynolds equation for the description of the hydrodynamic phenomena and the formulations of Hertz (1896) and Greenwood e Williamson (1966) of the contact between the asperities of rough surfaces. To address the complexity of the interactions between the fluid and the contacted solid pairs, the Lagrangian-Eulerian Arbitrary description is presented in this research. Through the Finite Element Method, a three-dimensional model is generated in Abaqus ®R to identify contact pressures, resulting tangential and normal stresses, and friction coefficients resulting from sliding between a textured and lubricated surface and a rigid plane (in analogy to classic contact models), whose roughness profiles are constructed from information on the quadratic roughness of dental surfaces. The sensitivity of some lubricant parameters in the determination of the coefficient of friction is also evaluated and models with different boundary conditions are proposed. However, for the vefrification of the latter, it is sought to qualitatively reproduce the result found by Lorentz (2013) in the numerical investigation of tribological systems without mixed lubrication regime. A methodology proposed here emerges as an effective alternative in the field of Tribology, in the prediction of the coefficient of friction and other relevant variables to a phenomenon still little understood. A sensitivity analysis of the modeling parameters is performed, in order to identify how they considerably affect the mechanical behavior at the contact interface.

The numerical modelling of the VKI-Longshot facility remains a challeng-ing task as it requires multi-physical numerical methods in order to simulate all the components. In the current dissertation, numerical tools were developed in order to study each component of the facility separately and a deep investigations of each stage of the shot were performed. This helped to better understand the different processes involved in the flow development inside this hypersonic wind tunnel. However the numerical computation of different regions of the facility treated as independent from each others remains an approximation at best.The accuracy of the rebuilding code for determining the free stream conditions and the total enthalpy in the VKI-Longshot facility was investigated by using a series of unsteady numerical computations of axisymmetric hypersonic flow over a heat flux probe. Good agreement was obtained between the numerical results and the measured data for both the stagnation pressure and the heat flux dur- ing the useful test time.The driver-driven part of the Longshot facility was modelled using the quasi one-dimensional Lagrangian solver L1d2. The three main conditions used for the experiments —low, medium and high Reynolds number —were considered.The chambrage effect due to the junction between the driver and the driven tubes in the VKI-Longshot facility was investigated. The computation showed great ben- efit of the chambrage in increasing the speed of the piston and thus the final compression ratio of the test gas.Two dimensional simulations of the flow in the driver and the driven tube were performed using Arbitrary Lagrangian Eulerian (ALE) solver in COOLFLuiD. A parallel multi-domain strategy was developed in order to integrate the moving piston within the computational domain.The computed pressure in the reservoir is compared to the one provided by the experiment and good agreement was obtained for both con- editions.Finally, an attempt was made to compute the starting process of the flow in the contoured nozzle. The transient computation of the flow showed how the primary shock initiates the flow in the nozzle before reaching the exit plan at time of 1.5 [ms] after the diaphragm rupture. The complex interactions of the reflected shocks in the throat raise the temperature above 9500 [K] which was not expected. Chemical dissociation of Nitrogen was not taken into account during this transient investigation which may play a key role considering the range of temperature reached near the throat. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Le travail réalisé dans cette thèse introduit le concept de l'utilisation de la méthode Level-Set pour simuler l'interface Glace/Air au cours du temps lors du processus de givrage en vol des aéronefs. Pour cela, un code de givrage tri-dimensionnel multi-blocs et parallélisé a été implémenté au sein du solveur NSMB (Navier-Stokes-Multi-Blocks). Il comprend notamment un module de calcul des trajectoires des gouttelettes par une approche Eulérienne compatible avec l'utilisation de grilles chimères et un module thermodynamique pour le calcul des masses de glace incluant deux modèles différents : un modèle algébrique itératif et un modèle à dérivées partielles. Une attention particulière a été portée sur la vérification du code de givrage implémenté en comparant systématiquement, si possible, les résultats obtenus avec les données expérimentales et numériques existantes dans la littérature. Pour cette raison, le module de déformation de maillage existant dans NSMB a été intégré au code implémenté afin de pouvoir simuler le givrage par une méthode traditionnelle. Enfin, un nouveau principe pour le suivi de l'interface glace/air est introduit via l'utilisation d'une méthode Level-Set. Puisque dans ce travail de thèse nous nous intéressons particulièrement au concept, la méthode Level-Set développée est d'ordre un et est résolue implicitement. On montrera cependant que des résultats valides sont obtenus avec une telle approximation. / This thesis introduces the concept of the Level-Set method for simulating the evolution through time of the ice/air interface during the process of in-flight aircraft icing. For that purpose, a three-dimensionnal multi-block and parallelized icing code have been implemented in the NSMB flow solver (Navier-Stokes-Multi-Blocks). It includes a module for calculating the droplet trajectories by an Eulerian approach compatible with the use of chimera grids and a thermodynamic module to calculate the ice masses including two different models : an iterative algebraic model and a PDE model. Particular attention was paid to the validation of the icing code irnplemented by comparing results with existing experimental and numerical data in the literature. For this reason, the existing mesh deformation algorithm in NSMB was integrated into the code to simulate icing by a traditional method. Finally a new principle to track the ice/air interface is introduced using the Level-Set method. Since we are particularly interested in the concept, the Level-Set method developped is first order and solved implicitly. However it will be shown that valid results are obtained with such an approximation.

Au sein des opérations de retraitement du combustible usé, la colonne pulsée à garnissage est l'appareil d'extraction liquide-liquide principalement utilisé. Dans un contexte de compétitivité économique et de raréfaction des ressources, l'efficacité de ces appareils est devenue un enjeu pour l'industriel. Afin d'améliorer leur rendement à travers un meilleur dimensionnement, la connaissance de la vitesse de glissement entre les phases de l'émulsion est nécessaire.Les travaux menés et présentés dans ce manuscrit s'articulent autour de la modélisation physique et numérique de l'hydrodynamique de l'émulsion ainsi que de sa caractérisation expérimentale.Dans ce travail, une modélisation d'approche eulérienne, inspirée des travaux de D. Lhuillier, permet l'obtention d'un modèle bi-fluide couplé à une équation d'évolution de la surface d'échange (aire interfaciale). La résolution du modèle s'effectue par éléments finis sous le logiciel CAST3M. A l'issue des calculs, le modèle montre sa capacité à restituer le comportement de l'émulsion et permet l'obtention des vitesses de glissement. Dans une optique de validation expérimentale du modèle, des expérimentations sur deux installations sont menées. Celles-ci font notamment intervenir un couplage entre les méthodes de vélocimétrie par image de particules et de fluorescence induite par laser afin d'obtenir les vitesses de chaque phases et le taux de rétention de la phase dispersée. Un algorithme de détection et de suivi de gouttes est développé afin d'obtenir la vitesse de la phase dispersée et sa fraction volumique. La confrontation de ces aux résultats numériques permet une première qualification encourageante du modèle. / In the core of spent fuel reprocessing operations, the pulsed columns with packing are the liquid-liquid extraction apparati mainly used. The context of economical competiveness and scarce resources, industrials are driven to improve the efficiency of these processes. Pulsed column efficiency is bound to the amount of available exchange surface, which depends on geometrical parameters of the column and the operating conditions. A better design would improve the efficiency. In this aim the knowledge of the interphase slip velocity is necessary. The work presented in this thesis revolves around physical and numerical modelling of the hydrodynamics of the emulsion and its experimental characterization.In this work, a eulerian approach, based on the work of D.Lhuillier, allows to obtain a two-fluid model coupled with an evolution equation of the exchange surface (interfacial area). We use finite elements method to solve this model along with CAST3M software. Numerical simulations have shown the model abilities to correctly reproduce the emulsion behaviour and to obtain the slip velocity.In order to experimentally validate the model, we carried out two types of experimentation. Particles images velocimetry coupled to laser induced fluorescence are involved to obtain velocities of each phases and the dispersed phase volume fraction. We developed a tracking algorithm to obtain the dispersed phase velocity and the hold up. These results, such as velocities and strain rate tensor, have been used in a first validation of the model.