Ce travail traite des méthodes de calcul numérique pour les simulations hydrodynamiques appliquées principalement sur des produits développés par ANDRITZ HYDRO. Il s’agit ici de mettre en place des schémas d’ordre élevé pour des simulations CFD en utilisant le code de calcul ASPHODEL développé et utilisé par ANDRITZ HYDRO. Les principales motivations sont l’augmentation de la fiabilité des résultats de calculs numériques avec un coût de calcul raisonnable. Cette fiabilité s’exprime à travers l’augmentation de la précision et de la robustesse des schémas numériques. Le code de calcul ASPHODEL est basé sur la méthode sans maillage SPH-ALE. Mélange entre les volumes finis et la méthode SPH (Smoothed Particle Hydrodynamics), la méthode SPH-ALE emploie un ensemble de points appelés particules servant à la discrétisation du domaine fluide. Elle permet en particulier de par son caractère sans maillage, un suivi des surfaces libres sans effort de calcul supplémentaire. Cet aspect est véritablement attrayant pour bon nombre d’applications industrielles notamment la simulation des écoulements à surface libre se produisant dans une turbine Pelton, mais également le remplissage d’une turbine Francis. Cependant, le bémol à cette méthode est son manque de précision spatiale. En effet les points de calcul étant mobiles, les opérateurs spatiaux doivent être en mesure de conserver leur précision et leur robustesse au cours du temps. La qualité des résultats en est du coup impactée, en particulier le champ de pression souvent excessivement bruité. La montée en ordre et l’amélioration de la consistance des opérateurs pour un vaste panel de configurations géométriques sont donc les enjeux de ce travail. En utilisant des outils inspirés par les volumes finis non-structurés, il est possible d’améliorer les opérateurs spatiaux. En effet, la montée en ordre ou p-raffinement peut notamment se faire avec des reconstructions d’ordres élevés pour évaluer les états aux interfaces des problèmes de Riemann. La sommation des flux numériques résolus par un solveur de Riemann est ensuite retravaillée pour obtenir un schéma numérique d’ordre global cohérent. Le même soucis de cohérence avec les schémas en temps doit d’ailleurs être pensé. Le gain de précision apporté par les schémas numériques d’ordre élevé est comparé avec un raffinement spatial, c’est à dire une augmentation du nombre des particules de taille plus petite, aussi appelé h-raffinement. La méthode SPH-ALE améliorée est ensuite testée sur des cas représentatifs des applications visées. En conclusion, les développements effectués dans cette étude ont été guidés par l’application en turbine Pelton principalement mais il va de soi qu’ils sont applicables à des écoulements sans surface libre dans les turbines Francis par exemple. Ce travail montre les possibilités d’une méthode sans maillage pour des cas d’écoulements complexes autour de géométrie tournantes. / This work deals with numerical methods for hydrodynamic testing applied mainly on products developed by ANDRITZ HYDRO. This is to put in place high order schemes for CFD simulations using the ASPHODEL calculation code developed and used by ANDRITZ HYDRO. The main reasons are the increased reliability of the results of numerical calculations with a reasonable computational cost. This reliability is expressed through increasing the accuracy and robustness of numerical schemes. The ASPHODEL computer code is based on the meshfree method SPH-ALE. Mix between finite volume method and SPH (Smoothed Particle Hydrodynamics), the SPH-ALE method uses a set of points called particles serving as the fluid domain discretization. It allows track free surfaces without additional computational effort. This is truly attractive for many industrial applications including the simulation of free surface flows occurring in a Pelton turbine, but also filling a Francis turbine. However, the downside of this method is its lack of spatial accuracy. Indeed calculation points are mobile, space operators must be able to keep their accuracy and robustness over time. The quality of results is impacted especially the pressure field is often excessively noisy. The rise in order and improving the consistency of the operators for a wide range of geometric configurations are the challenges of this work. Using tools inspired by the unstructured finite volume, it is possible to improve the spatial operators. Indeed, the increasing order or p-refinement particular can be done with reconstructions of high order to assess the conditions at the interfaces of Riemann problems. The summation of discret fluxes solved by Riemann solver is then reworked to obtain a coherent global order scheme. The same concern for consistency with temporal schemes should also be considered. The precision gain provided by numerical schemes of higher orders is compared with a spatial refinement ie an increase in the number of smaller particles ; also called h -refinement . Improved SPH -ALE method is then tested on representative cases of intended applications. In conclusion, the developments made in this study were guided in accordance mainly with the Pelton turbine but it goes without saying that they are applicable to non- free surface flows in Francis turbines for example. This work shows the possibilities of a free mesh method for cases of complex flow around rotating geometry.

La cavitation peut avoir lieu dans les turbines hydrauliques. Ce phénomène se produit lorsque les bulles de vapeur s’effondrent à proximité de la surface de la machine. Ceci entraîne des conséquences négatives, telles que l’érosion, affectant ainsi les performances de la machine. L’effondrement d’une bulle de gaz non-condensable dans l’eau est simulé en utilisant la méthode sans maillage SPH-ALE, qui intègre un modèle pour simuler les écoulements compressibles et multiphases. Le modèle résout les équations de conservation de masse, de quantité de mouvement et d’énergie du système d’Euler, en utilisant l’équation d’état de Stiffened Gas pour l’eau et l’équation d’état de gaz parfait pour le gaz non-condensable à l´ıntérieur de la bulle. Les deux phases sont modélisées comme compressibles et le changement de phase n’est pas considéré. La caractéristique sans maillage de la méthode SPH-ALE permet le calcul des écoulements diphasiques où l’interface est nettement définie. Pour les applications de cavitation, où le nombre de Mach atteint des valeurs de 0.5, la distribution de particules doit être corrigée. Cela est réalisé grâce à la fonctionnalité ALE. Le modèle compressible a été validé à l’aide de configurations monodimensionnelles, comme le cas du tube à choc pour des écoulements monophase et multiphases. L’effondrement de la bulle près d´une paroi a été abordé comme le mécanisme fondamental qui produit des dégâts. Son comportement général se caractérise par la formation d’un micro jet d’eau et par l’effondrement de la bulle sur elle-même. Le phénomène est analysé en tenant compte des principaux paramètres qui le régissent, comme la distance initiale entre le centre de la bulle et la paroi (H0), la taille de la bulle (R0) et le taux de pression qui entraîne l’effondrement (pw/pb). Il est démontré que l’intensité de l’effondrement dépend principalement du rapport de pression entre le liquide et la bulle (pw/pb). De plus, quatre indicateurs, comme la pression en paroi, l’impulsion, la pression du coup de bélier et la vitesse du micro jet d’eau, servent à déterminer le chargement. Cette analyse indique qu’une bulle initialement située à une distance inférieure à H0/R0 = 2 présente un haut potentiel d’endommagement. Afin de prédire cet endommagement, la mécanique du solide est analysée à l’aide de simulations d’interaction fluide-structure. On obtient que le matériau réagit aux charges hydrauliques en ayant des zones de compression et de traction. Ceci suggère qu’un mécanisme de fatigue entraîne le phénomène d’endommagement. En plus, on constate que les contraintes les plus importantes sont situées sous la surface du matériau, indiquant que cette zone peut être sujette à une déformation plastique. / Hydraulic turbines can experience cavitation, which is a phenomenon occurring when vapor bubbles collapse in the vicinity of the machine’s surface. This phenomenon can lead to negative consequences, such as erosion, that affect the machine’s performance. The compression of a non-condensable gas bubble in water is simulated with the Smoothed Particle Hydrodynamics method following the Arbitrary Lagrange Euler approach (SPHALE), where a compressible and multiphase model has been developed. The model solves the mass, momentum and energy conservation equations of the Euler system using the Stiffened Gas EOS for water and the ideal gas EOS for the non-condensable gas inside the bubble. Both phases are modeled as compressible and the phase change is not considered. The meshless feature of the SPH-ALE method allows the calculation of multiphase flows where the interface is sharply defined. For cavitation applications, where the Mach number reaches values of 0.5, the distribution of particles must be corrected, which is achieved by the ALE feature. The compressible model was validated through monodimensional configurations, such as shock tube test cases for monophase and multiphase flows. The bubble compression close to the wall has been addressed as the fundamental mechanism producing damage. Its general behavior is characterized by the formation of a water jet and by the collapse of the bubble by itself. The phenomenon is analyzed by considering the major parameters that govern the bubble collapse dynamics, such as the initial distance between the bubble center and the wall (H0), the bubble size (R0), and the collapse driven pressure ratio (pw/pb). It is shown that the intensity of the collapse depends mainly on the pressure ratio between the liquid and the bubble (pw/pb). As well, four indicators, such as the pressure at the wall, the impulse, the water-hammer pressure and the water jet velocity, are used to determine the loading. This analysis gives that the bubble initially located at a distance lower than H0/R0 = 2 presents high potential to cause damage. In order to predict the damage due to the bubble collapse, the solid mechanics is analyzed through fluid-structure interaction simulations. It is obtained that the material reacts to the hydraulic loads by having compression and traction zones, suggesting that a fatigue mechanism drives the damage phenomenon. Additionally, it is found that the highest stresses are located below the material surface, indicating that this zone may reach plastic deformation.

Dans de nombreuses applications aérospatiales, on peut trouver des films liquides cisaillés, c'est-à-dire une fine couche liquide qui ruisselle sur une paroi entrainée par le gaz. Par exemple, une couche de liquide peut se développer sur la voilure des avions, givrer et dégrader les performances. Des vagues peuvent se développer à l'interface liquide-gaz, et l'analyse correcte de ces instabilités devient très importante pour modéliser ce phénomène physique. En effet, la présence d'instabilités modifie les échanges liquide-gaz, notamment les transferts de masse et chaleur. Le but de cette thèse est de développer une technique permettant de coupler la phase gazeuse afin de reproduire les interactions à l'interface. La couche de liquide étant beaucoup plus mince que celle du gaz, une approche intégrale sur l'épaisseur est utilisée pour la modélisation. Enfin, deux cas d'un écoulement diphasique se développant dans une conduite confinée et dans une conduite plus large sont étudiés. les résultats sont ensuite comparés à des autres méthodes de référence, plus coûteuses en temps de calcul. / In many aerospace applications one can find liquid films sheared by a gas flow. In example, these liquid sheets can develop on aircraft wings, freeze and then destroy the aerodynamics performances. Waves can develop at the liquid-gas interface, and the correct analysis of such instabilities becomes very important to model this physical phenomenon. Indeed, instabilities mdify liquid-gas exchanges, such as mass and heat transfers. The aim of the present work consists in developing a technique to couple the liquid phase to the gas phase in order to reproduce the interactions at the interface. Since the liquid layer is much thinner then the gas, anintegral approach is used for modeling. Finally, two cases of a two-phase flow developing in a strictly confined channel and in a large channel are studied. Results are then compared to other reference methods which are more expensive in terms of computational cost.

Made available in DSpace on 2014-07-29T14:57:51Z (GMT). No. of bitstreams: 1 THIAGO_VAZ_DISSERTACAO - TEXTO PARTE 1.pdf: 3764953 bytes, checksum: ca2c1d973d33ea6e46037d52fa0b3640 (MD5) Previous issue date: 2011-05-19 / RFID has brought to the virtual world a high volume of events that previously existed only in the physical world, boosting the development of technologies for Complex Event Processing (CEP). In this work, we investigated the CEP for RFID and found that the main standards in the area, defined by EPCglobal, do not address adequately the application needs for handling complex events. So we developed the CEPFID Middleware, a software platform that simplifies the development of RFID applications providing the discovery of complex events in the EPC Network. Due to need to adapt the middleware to different application scenarios of RFID technology, we designed an extensible architecture that enables the inclusion of new services for the processing and discovery of complex events. / O RFID trouxe ao mundo virtual um grande volume de eventos que anteriormente só existiam no mundo físico, impulsionando o desenvolvimento das tecnologias de Processamento de Eventos Complexos (CEP). Neste trabalho nós investigamos o CEP para RFID e descobrimos que os principais padrões na área, definidos pela EPCglobal, não contemplam satisfatoriamente as necessidades das aplicações em manipular eventos complexos. Portanto, desenvolvemos o Middleware CEPFID, uma plataforma de software que facilita o desenvolvimento de aplicações RFID proporcionando a descoberta de eventos complexos na Rede EPC. Devido à real necessidade de adequação do middleware aos diferentes cenários de aplicação da tecnologia RFID, projetamos uma arquitetura extensível que possibilita a inclusão de novos serviços para o processamento e descoberta de eventos complexos.

Un dispositif expérimental de dépôt électrochimique en couches atomiques (EC-ALE) a été conçu et réalisé dans cette thèse. Sa mise en oeuvre a permis de valider cet équipement. Grâce à la flexibilité de cet équipement, un nouveau point de vue du dépôt en sous potentiel (UPD) du cobalt sur un substrat d'or a été proposé. Les résultats montrent aussi que le dépôt alternatif de monocouches de Co et de Sb est possible. Pour la première fois, l'électrodéposition de films minces Sb2Se3 par la méthode EC-ALE sur des électrodes poly-cristallines d'or a été obtenue et étudiée. Les paramètres de dépôt ont été déterminés et le film mince déposé a été caractérisé par analyses MEB et RAMAN. L'adsorption irréversible et le phénomène UPD réversible de Sb sur une électrode de Pt a aussi été étudié. Les résultats montrent qu'après la réduction en Sb des espèces adsorbées SbO+, des atomes de Sb peuvent être déposés sur cette électrode de Pt modifiée au Sb en utilisant le dépôt en sous potentiel pour augmenter le recouvrement du substrat du Pt par le Sb. / An electrochemical atomic layer epitaxy (EC-ALE) experiment platform was designed and constructed in this thesis, and this platform was proved to be qualified for EC-ALE experiments.Benefiting from the flexibility of the EC-ALE equipment, a new viewpoint about the UPD behavior of cobalt on the gold substrate has been put forward in this work. The results also show that the subsequent alternate deposition of Co and Sb monolayers is feasible.For the first time the electrodeposition of Sb2Se3 thin films by EC-ALE method on polycrystalline Au electrodes has been obtained and investigated. The deposition parameters were determined and the deposit was characterized by SEM and Raman analysis.The irreversible adsorption and reversible UPD behaviour of Sb on Pt electrode have also been studied. The results show that after the irreversibly adsorbed SbO+ species are reduced to metallic Sb, Sb atoms can be further deposited onto this Sb-modified Pt electrode in the way of UPD to increase the coverage of Sb on the Pt substrate.

L'utilisation croissante des sources d'énergie renouvelable avec une grande volatilité de production, comme l'énergie éolienne et solaire, conduit à des fluctuations dans le réseau électrique qui doivent être compensées. Pour cette raison les machines hydrauliques, turbines et pompes, sont plus souvent opérées dans les régimes de fonctionnement hors fonctionnement nominal et la fréquence des phases de démarrage et arrêt augmente. Ce type de fonctionnement peut avoir des conséquences importantes sur le cycle de vie des machines. Il est donc essentiel de prendre en compte l'écoulement dans les phases transitoires lors de la conception de la machine et la simulation numérique des écoulements est un outil adapté pour cela. La présente étude a pour objectif de développer une méthode de couplage flexible qui combine la méthode à maillage volumes finis (VF) et la méthode sans maillage Smoothed Particle Hydrodynamics - Arbitrary Lagrange Euler (SPH-ALE). Cette méthode couplée peut être utilisée comme outil pour l'investigation des phénomènes transitoires dans les machines hydrauliques. SPH-ALE est particulièrement bien adapté aux simulations des écoulements fortement dynamiques avec des géométries mobiles mais elle a des difficultés pour calculer des forts gradients de pression et vitesse. Un raffinement de particules est difficile à implémenter, surtout si les particules doivent être raffinées de manière anisotrope. Les méthodes volumes finis (VF) sont établies pour les simulations numériques d'écoulements grâce à leur stabilité et précision. Par contre, elles peuvent être lourdes pour les simulations avec des géométries mobiles et demandent souvent une interface entre des parties mobiles et statiques du maillage ce qui génère des erreurs supplémentaires. Pour combiner les deux approches complémentaires, une méthode de couplage a été développée qui décompose le domaine de calcul en zones où la vitesse et la pression sont calculées par la méthode VF, en zones où elles sont obtenues par SPH-ALE et en zones de recouvrement où les informations sont transférées de la zone VF à la zone SPH et inversement. Dans les zones de recouvrement les points de calcul VF sont utilisés comme voisins pour l'intégration en espace des particules SPH. Aux limites du maillage VF la vitesse et la pression sont interpolées des particules SPH, similairement aux méthodes Chimére des maillages recouvrants. Un logiciel SPH-ALE existant du groupe ANDRITZ est utilisé pour cette étude. Un solveur VF faiblement compressible est implémenté dans ce logiciel. Le solveur discrétise la même forme des équations de Navier-Stokes que le solveur SPH-ALE. Des solveurs de Riemann avec des états reconstruits par la méthode MUSCL sont employés. En outre, le solveur SPH-ALE est amélioré et adapté aux écoulements internes. Pour cette raison des conditions à l'entrée et à la sortie du type subsonique sont implémentées. Du plus, une méthode de correction du gradient de la fonction kernel est présentée qui améliore la précision du champ de pression, notamment si les particules ne sont pas distribuées régulièrement. La méthode couplée est validée à l'aide des cas test académiques en unidimensionnel et en bidimensionnel, comme le cas de tube à choc, les tourbillons de Taylor-Green et l'écoulement autour d'une aube symétrique du type NACA avec des particules en description eulérienne. En outre, le couplage offre la possibilité d'imposer des conditions à la sortie aux particules lagrangiennes. La méthode est appliquée aux simulations d'écoulement transitoire en 2D avec des particules qui se déplacent en suivant les géométries mobiles. / The increased use of intermittent forms of renewable energy like wind and solar energy produces fluctuations in the electric grid that have to be compensated. For this reason, hydraulic machines like turbines and pumps are more often operated under non-conventional operating conditions and are submitted to frequent starts and stops. This type of operating conditions has important consequences on the life cycle of the machines. It is thus of paramount importance that transient flows at off-design conditions are properly taken into account in the design phase and numerical simulation is an appropriate way to do so. The present study aims at developing a flexible coupling method of the meshbased Finite Volume Method (FVM) and the meshless Smoothed Particle Hydrodynamics - Arbitrary Lagrange Euler (SPH-ALE) method, which can be used as a tool for the investigation of transient phenomena in hydraulic machines. SPH-ALE is very well adapted for the simulation of highly dynamic flows with moving geometries but has difficulties to correctly represent rapidly changing gradients of the field variables. Particle refinement is difficult to implement, especially if particles are refined in an anisotropic way. FV methods are well established in CFD because of their accuracy and stability. However, they can be tedious for simulations with moving geometries and often necessitate an interface between moving and static parts of the mesh which introduces additional errors. To overcome the shortcomings of both methods, a coupling method is developed that uses a decomposition of the computational domain into regions where the physical field variables are computed by the FV method, regions where they are computed by SPH-ALE and overlapping regions where the information is transferred from the FV domain to the SPH domain and vice versa. In the overlapping regions FV calculation points are used as neighbors for the SPH integration in space. At the boundaries of the FV mesh, velocity and pressure are interpolated from the SPH particles by means of scattered data interpolation techniques, similarly to Chimera methods for overlapping grids. For this study, an existing SPH-ALE software of the ANDRITZ Group is used. A weakly compressible FV solver is implemented into this software that discretizes the same form of the Navier-Stokes equations than the SPH-ALE solver. Similar to the present SPH-ALE method, Riemann solvers with reconstructed states, obtained by a MUSCL scheme, are employed. Moreover, adaptations and improvements of the SPH-ALE solver itself are made, which are important for the coupling and for the simulation of internal flows in general. Thus, subsonic inlet and outlet conditions are implemented. Furthermore, a correction method of the kernel gradient is presented that ensures zeroth order consistency of the SPH-ALE approximation of the divergence of the convective fluxes. The correction improves greatly the SPH pressure field on non-uniform particle distributions. The implemented coupled method is successfully validated by means of inviscid academic one-dimensional and two-dimensional testcases like a shock tube case, Taylor-Green vortices and the flow around a symmetric NACA airfoil with particles in Eulerian description. Furthermore, the coupling provides a possibility to implement outlet boundary conditions to Lagrangian moving SPH particles. It is then applied to the simulation of transient flows in rotor stator systems in 2D with moving particles.

Depuis l'invention du soudage par frottement malaxage (FSW) en 1991 par The Welding Institute (TWI), l'industrie aéronautique envisage de remplacer l'assemblage par rivetage des tôles par ce procédé. L'objectif est d'alléger les structures, plus particulièrement les panneaux des structures d'aéronefs qui sont souvent raidis avec des lisses et des profilés assemblés à la peau extérieure. La qualification du procédé FSW dans l'aéronautique requiert de produire des joints de très haute qualité avec une longue durée de vie des outils. Ceci nécessite un réglage fin des paramètres procédé tels que la géométrie de l'outil (dimensions de l'épaulement et du pion, filetage) et les vitesses d'avance et de rotation. Le but de ce travail est d'accompagner la conception du procédé par la simulation numérique. Il comporte trois parties: (1) le développement de la formulation arbitrairement lagrangienne ou eulérienne (ALE) dans la version parallèle du logiciel Forge® pour l'adapter au FSW, (2) la modélisation du procédé FSW et (3) la calibration du modèle EF et la validation des résultats numériques suite à la réalisation d'essais FSW en configuration par recouvrement.Ce travail est basé sur la formulation ALE parallèle développée dans Forge® afin de modéliser les différents défauts caractéristiques du FSW (bavures et trous). L'algorithme de transport des champs continus, indispensable au caractère eulérien, est affiné afin de considérer les particularités de l'écoulement du procédé (nœud du plan d'entrée, en contact avec l'outil ou dans les bavures et sortant du domaine d'étude). Un nouveau schéma d'intégration temporelle basée sur les coordonnées cylindriques plus adapté à ce procédé fortement rotationnel est introduit. Enfin, un filetage sur l'outil est modélisé par l'introduction d'une nouvelle loi de frottement à l'interface outil / matière.Cette formulation ALE générale et parallèle a montré sa robustesse pour modéliser le procédé FSW. Les résultats thermomécaniques sont en accord avec de précédents résultats validés par l'expérience. De plus l'outil numérique obtenu a montré ses capacités pour la modélisation des bavures ou des pertes de contact derrière le pion.Finalement la configuration de soudage par recouvrement industrielle a été étudiée. Pour cela des essais expérimentaux ont été réalisés pour divers paramètres du procédé et divers types de joints de recouvrement. Puis les mesures de forces et du couple dans l'outil et de la température dans l'outil et la tôle ont permis de calibrer les paramètres du modèle (frottement, loi de comportement, coefficients thermiques) et de valider les résultats EF. / Since the invention of Friction Stir Welding (FSW) in 1991 by The Welding Institute (TWI), aeronautics industry has been investigating the possibility to use this process instead of riveting with the objective to lighten its structures and more particularly the aircrafts structure panels. Aircraft panels are often straightened with stringers and profiles which are joined to the outer skin. The qualification of FSW process in aeronautics requires producing very high quality joints with the longevity of tools, which requires fine tuning of process parameters such as the geometry of the welding tool (shoulder and pin dimensions, threads on pin and shoulder) and the advancing and rotating speeds. The aim of this work is to support the design of the process by numerical simulation. It has three parts: (1) developing an efficient and accurate Arbitrary Lagrangian or Eulerian (ALE) formulation within the parallel version of Forge® software, (2) modelling the FSW process and (3) calibrating the F.E. model and validating simulation results thanks to FSW experiments on lap joints.This work is based on the parallel ALE formulation developed in Forge® to model the different possible defects taking place in FSW (flashes and worm holes). The transport algorithm of continuous fields, required by the Eulerian frame, is enhanced to take the special characteristics of the FSW's flow into account (nodes located in input plan or flashes or in contact with the tool). A new time integration scheme based on cylindrical coordinates, which are best suited for this process, is introduced. Finally, the pin and shoulder threads are modelled by introducing a new friction law at the tool / material interface.This general and parallel ALE formulation is robust enough to model the FSW process. Thermomechanical results obtained are in agreement with previous results validated by experiences. And the numerical tool demonstrated its ability to model flashes formation and losses of contact behind the pin. Finally industrial welding lap joints configuration was studied. Experimental tests were conducted with several process parameters and type of lap joint. And measure of torque and forces in tool, and temperature in tool and sheet metal allowed to calibrate model parameters (friction, behavior law, thermal coefficients) and to validate FE results.

Obsahem této práce je numerická simulace dvoudimenzionálního proudění nestlačitelné vazké kapaliny. Uvažujeme rotující elipsu soustředně umístěnou v kružnici. Prostor mezi elipsou a kružnicí je vyplněn kapalinou. Cílem je popsat proudění kapaliny vyvolané otáčející se elipsou, tzn. stanovit rychlostní pole a rozložení tlaku. Dále pak chceme stanovit přídavné silové účinky kapaliny působící na elipsu. Tyto výsledky získáme řešením Navierových-Stokesových rovnic metodou konečných prvků. Důraz je kladen na odvození numerického schématu v maticové formě vhodné pro numerickou implementaci. Časově závislá výpočetní síť je popsána pomocí Arbitrary Lagrangian-Eulerian (ALE) formulace. Pro obdržení relevantních výsledků je nutná stabilizace metody konečných prvků. Uvedené výsledky naznačují, že odvozená metoda je dostatečně přesná.

Dans les approches numériques pour l'usinage à grande vitesse, le comportement rhéologique des matériaux usinés est généralement décrit par une loi de Johnson Cook et le frottement à l’interface par un coefficient constant de type Coulomb. Une première approche propose de déterminer une loi de frottement à paramètres multiples pour des basses températures combinées à des vitesses de glissement élevées au moyen d'un tribomètre. Les résultats expérimentaux sont comparés à un modèle numérique et une méthode inverse est utilisée pour minimiser l'erreur entre des simulations numériques et expérimentales sur les forces tangentielles et normales. Cette méthode permet de récupérer un coefficient de frottement de type Coulomb qui est associé à la pression locale, la température et la vitesse de glissement. La réalisation de plusieurs essais fournit une loi de frottement à paramètres multiples pour des vitesses de glissement élevées et des basses températures. Une seconde étude est menée sur les phénomènes microstructuraux intervenant à l’interface outil-copeau. De nombreuses études ont montré que les phénomènes de recristallisation dynamique apparaissent lors de l'usinage dans l'interface outil-copeau. La loi de Johnson Cook ne comprend pas de tels phénomènes. Ainsi, les modèles rhéologiques spécifiques basés sur la métallurgie sont introduits pour tenir compte de ces phénomènes de recristallisation dynamique. Un modèle éléments finis de la coupe orthogonale à deux dimensions est développé avec le logiciel Abaqus Explicit en utilisant une formulation ALE. Ce modèle éléments finis peut alors prédire la formation des copeaux, les températures d'interface, les longueurs de contact et les forces de coupe. De nombreux essais spécifiques sont réalisés sur un banc d'essai de coupe orthogonale et un tribomètre grande vitesse sur un acier AISI 1045 et avec un outil en carbure non revêtu. Les résultats expérimentaux sont ensuite comparés aux simulations numériques. Ces deux approches montrent qu’il est nécessaire, afin d’optimiser les modèles de coupe orthogonale, de prendre en compte une loi de frottement à paramètres multiples qui tient compte des pressions, températures et vitesses de glissement locales et d’intégrer une loi de comportement rhéologique à base métallurgique. / In numerical approaches for high speed machining, the rheological behavior of machined materials is generally described by Johnson Cook law and the friction at the interface by a constant coefficient of Coulomb. A first approach proposes to determine a multiparameter friction law for low temperatures in combination with high sliding speeds by means of a tribometer developed by Meresse et Al. [Mer11]. The experimental results are compared to a numerical model and a inverse method is used to minimize the error between the numerical and experimental simulations on tangential and normal forces. This method allows to recover a Coulomb friction coefficient. This one is associated with the local pressure, temperature and sliding velocity. Several tests provide a multiparameter friction law for higher sliding speeds and low temperatures. A second study is conducted on the microstructural phenomena occurring at the tool-chip interface. Numerous studies have shown that the dynamic recrystallisation phenomena appear during machining at the interface. Johnson Cook law excludes such phenomena. Thus, specific rheological models based on metallurgy are introduced to take into account the dynamic recrystallization phenomena. Two empirical models proposed by Kim et al. [Kim03] and Lurdos [Lur08] are studied. A two-dimensional finite element model of the orthogonal cutting is developed with Abaqus Explicit software using an ALE formulation. This finite element model can predict chip formation, interface temperatures, contact lengths and cutting forces. Many specific tests are performed with an orthogonal cutting test bench and with an high speed tribometer on an AISI 1045 steel. The experimental results are then compared with numerical simulations. Both approaches show that it is necessary to optimize the orthogonal cutting model, to take into account a multi-parameter friction law that considers the local pressures, temperatures, and sliding velocities and to integrate a rheological behavior law based on mettalurgy.

New approach for computational modeling of the dynamic behavior of elastic body immersed in incompressible viscous stagnant fluid is described in this work. It is based on determination of added effects (added mass and added damping). This effects are inserted into computational model and it replace influence of the fluid. Commonly used commercial computational software may be used. Approach is based on assumption appropriate for the linear flow. Two pressure field are determined. One for movement of the unite acceleration of the fluid boundary and the second for unite velocity. Nonlinear model (Navier-Stokes equation in ALE form) had to be used for determination of the added damping, hence results are valid only for pre-selected amplitude of vibration.