Le prototypage virtuel des circuits pneumatiques de puissance, par exemple les circuits de freinage des véhicules industriels, constitue un enjeu important en raison de la complexité des écoulements en régime transsonique et des couplages entre les échelles locales et macroscopiques. Ces problèmes sont rencontrés lors de la conception, de la synthèse des commandes et de l'analyse des performances statiques et dynamiques de ces circuits et l'analyse. La mise au point des modèles numériques de ces systèmes induit des coûts et des temps importants par rapport à d'autres systèmes. La démarche proposée dans cette thèse repose sur la construction numérique de bases de données permettant de caractériser le comportement local et macroscopique d'un composant de circuit en fonction de la variation de certains paramètres physiques ou géométriques par rapport à un point de fonctionnement de référence. Les bases de données résultent de l'extrapolation de la solution des équations de Navier Stokes moyennées (RANS) pour le point de référence considéré obtenu à l'aide d'un logiciel de paramétrisation en mécanique des fluides (Turb'Opty). La contribution de cette thèse repose pour l'essentiel dans un travail d'analyse des solutions issues de la paramétrisation dans deux contextes différents: la tuyère De Laval et un élément "coude", des composants élémentaires de circuit. Nous avons montré que ces exemples "simples" conduisent déjà à des difficultés importantes en termes de paramétrisation du problème et du calcul des dérivées des champs aérodynamiques en raison de la taille du problème. Pour pallier cette difficulté, nous avons proposé de déraffiner le maillage et nous avons alors montré que cette démarche conduit parfois à déplacer ou à atténuer certains phénomènes (chocs). La deuxième contribution de ce travail repose sur l'évaluation de la qualité des solutions extrapolées, de leur domaine de validité et la construction des liens entre grandeurs locales et macroscopiques. Nous avons enfin proposé une démarche permettant de reconstruire la caractéristique en débit d'un composant à partir de la détermination de la solution extrapolée pour un nombre limité de points de référence.

Cette étude s’intéresse à la modélisation de la ségrégation des impuretés lors des procédés de solidification dirigée du silicium pour l’industrie photovoltaïque. Il s’agit d’un problème multi-physiques et multi-échelles qui nécessite des modèles efficaces pour pouvoir traiter des configurations industrielles en 3D. La première partie de l’étude porte sur le développement de fonctions de paroi solutales dans le but d’estimer la ségrégation sans résoudre numériquement la couche limite solutale. Un modèle analytique est alors utilisé pour estimer le paramètre convecto-diffusif à partir de la contrainte de frottement à l’interface solide/liquide. L’approche proposée consiste à coupler une simulation hydrodynamique de la convection dans la phase liquide et un calcul analytique de la ségrégation. Cette démarche est validée sur un cas de référence en 2D. Le modèle développé fournit une estimation pertinente de la concentration dans un lingot solidifié avec une vitesse imposée, y compris sur un maillage ne résolvant pas la couche limite solutale. La deuxième partie de l’étude concerne l’utilisation d’un système de brassage mécanique dans le but de favoriser la ségrégation des impuretés. Des simulations transitoires de brassage sont réalisées avec le logiciel FLUENT. Les résultats sont comparés à des mesures de champ de vitesse par PIV effectuées sur une maquette en eau, afin de valider le modèle hydrodynamique. La simulation de brassage est ensuite couplée à un calcul de ségrégation en régime quasi-permanent qui permet d’analyser l’influence de l’écoulement sur la couche limite solutale. Dans la dernière partie de l’étude, une simulation de solidification dirigée 3D instationnaire en régime de convection forcée, intégrant la thermique du four, est réalisée. Un modèle empirique de forces volumiques est alors proposé pour décrire l’écoulement lié au brassage dans le silicium liquide. Une première tentative d’estimation des ségrégations par l’approche analytique est ensuite mise en œuvre. / The present study focuses on solute segregation during photovoltaic silicon directional solidification. This multi-physics problem involves various spatial and temporal scales. The numerical simulation of this process requires efficient models, especially for 3D industrial configurations. In the first part of the study, solute wall functions are derived from a scaling analysis in order to estimate the segregations without numerical resolution of the solute boundary layer. The method is based on the coupling of an hydrodynamic simulation of convection in the liquid phase and an analytical segregation computation. The developed analytical model provides an estimation of the convecto-diffusive parameter from the wall shear-stress at the solid/liquid interface. A reference case in 2D with imposed solidification rate is used for validation purposes. The developed model provides a meaningful estimation of concentration fields in the ingots. In a second part, we focus on segregation optimization by a mechanical stirrer. Transient stirring simulations, using a sliding mesh technique, are achieved with FLUENT commercial software. Results are compared to PIV velocity field measurements performed on an experimental setup using water. A segregation computation in a quasi-steady regime is then implemented in the stirring simulation. The effect of the stirring parameters are directly observed on the solute boundary layer at the solid/liquid interface. In a third part, a transient solidification simulation, including furnace thermal conditions, is performed in a 3D configuration with forced convection. The flow generated by the impeller is described thanks to an empirical model based on body forces. A first attempt is finally made to retrieve segregations in the ingot with the developed analytical method.

Le but du travail présenté est la compréhension puis la simulation numérique des mécanismes d'érosion des augets de turbine Pelton par impacts répétés de gouttes d'eau dans le but de prédire la durée de vie des composants. Tout d'abord, les phénomènes de propagation d'ondes dans les milieux fluide et solide sont étudiés. Cela permet de mettre en lumière l'évolution temporelle et la distribution spatiale des pressions de contact, et l'apparition de microjets par éjection supersonique du fluide au contact. Les études expérimentales de l'érosion par gouttes d'eau traduisent un dommage basé sur la fissuration par fatigue. Des simulations numériques en dynamique rapide couplées fluide-structure sont alors effectuées. Le domaine solide est discrétisé par la Méthode des éléments Finis (MEF), et le domaine fluide par la méthode Smoothed Particle Hydrodynamics (SPH), qui est une méthode particulaire (sans maillage) particulièrement adaptée aux grandes distorsions et au suivi des surfaces libres. L'analyse des états de contraintes vient corroborer la nature cyclique de l'endommagement. La simulation d'érosion est alors réalisée à l'aide de critères de fatigue multiaxiaux. Le choix se porte vers un premier critère général de l'American Society of Mechanical Engineers (ASME), utilisant les valeurs principales des différences de contraintes au cours du temps. Le second choix concerne un critère à plan critique : le critère de Dang Van 2. Il traite séparément la contrainte hydrostatique et le cisaillement alterné maximal local. Ces critères permettent de définir les régions érodées du solide au bout d'un nombre d'impact donné, ce qui fait de cette démarche une méthode prédictive. Une étude paramétrique pour différentes tailles de gouttes et vitesses d'impact est ensuite réalisée, puis on évalue l'influence de la présence d'une couche de coating. / The goal of this work is the comprehension and the numerical simulation of erosion caused by repeated droplet impact on Pelton turbine buckets, to predict the lifetime of these components. First, waves propagation phenomenon inside fluid and solid domains are presented, which allows determining the time evolution and spatial distribution of contact pressure, and the birth of lateral microjets by supersonic ejection of the fluid on the contact. Experimental studies of erosion by droplet impact highlight a fatigue cracking-based erosion mechanism. Then, coupled FSI computation are performed. The solid subdomain is discretized by the Finite Element Method (FEM), and the fluid subdomain by the Smoothed Particle Hydrodynamics (SPH), which is a particle method (meshless) effectively recommended for large distortions and free surface tracking. Stress analysis confirms the cyclic nature of the damage mechanism, and erosion simulation is performed using multiaxial fatigue criteria. The first selected criterion is a general one from the American Society of Mechanical Engineers (ASME) using principal values of stress differences over time. The second one is the Dang van 2 criterion, belonging to the family of critical plane criteria. This criterion considers separately the effects due to hydrostatic stress on one hand, and the ones induced by maximum local shear on the other. These two criteria are used to building the equivalent eroded zones of the solid subdomain for a given number of impacts, which allows to qualify this procedure as a predictive predictive. Finally, a parametric study for different droplet sizes and velocites is computed, and the effects of a coating layer are investigated.

Le prototypage virtuel des circuits pneumatiques de puissance, par exemple les circuits de freinage des véhicules industriels, constitue un enjeu important en raison de la complexité des écoulements en régime transsonique et des couplages entre les échelles locales et macroscopiques. Ces problèmes sont rencontrés lors de la conception, de la synthèse des commandes et de l'analyse des performances statiques et dynamiques de ces circuits et l'analyse. La mise au point des modèles numériques de ces systèmes induit des coûts et des temps importants par rapport à d'autres systèmes. La démarche proposée dans cette thèse repose sur la construction numérique de bases de données permettant de caractériser le comportement local et macroscopique d'un composant de circuit en fonction de la variation de certains paramètres physiques ou géométriques par rapport à un point de fonctionnement de référence. Les bases de données résultent de l'extrapolation de la solution des équations de Navier Stokes moyennées (RANS) pour le point de référence considéré obtenu à l'aide d’un logiciel de paramétrisation en mécanique des fluides (Turb’Opty). La contribution de cette thèse repose pour l'essentiel dans un travail d'analyse des solutions issues de la paramétrisation dans deux contextes différents: la tuyère De Laval et un élément "coude", des composants élémentaires de circuit. Nous avons montré que ces exemples "simples" conduisent déjà à des difficultés importantes en termes de paramétrisation du problème et du calcul des dérivées des champs aérodynamiques en raison de la taille du problème. Pour pallier cette difficulté, nous avons proposé de déraffiner le maillage et nous avons alors montré que cette démarche conduit parfois à déplacer ou à atténuer certains phénomènes (chocs). La deuxième contribution de ce travail repose sur l'évaluation de la qualité des solutions extrapolées, de leur domaine de validité et la construction des liens entre grandeurs locales et macroscopiques. Nous avons enfin proposé une démarche permettant de reconstruire la caractéristique en débit d'un composant à partir de la détermination de la solution extrapolée pour un nombre limité de points de référence. / Virtual prototyping pneumatic circuits for power transmission, for example braking circuits of trucks, is still a difficulty because of the complexity of the flow behavior in transonic conditions and of the coupling between local and macroscopic scales. These problems are met during system design, control synthesis and for static and dynamic performance analysis. Tuning accurate numerical models requires important costs and time when compared to other systems. The methodology proposed in this PhD thesis relies on numerically determining a data base that characterizes the local and macroscopic behavior of a circuit component according the variation from a reference point of some physical or geometrical parameters. The data bases are obtained from the extrapolation of the Mean Navier Stokes solution (RANS) for a given reference point with the help of a parametrization software dedicated to fluid mechanics (Turb'Flow). The main contribution of this thesis relies io the analysis of the solution obtained from the parametrization in two different cases: the De Laval nozzle and un "elbow" connecting element, which are elementary component in a circuit. We have shown that these two "simple" cases lead already to important difficulties in term of problem parametrization and calculation of the derivatives of the aerodynamic fields because of the problem dimension. In order to tackle this, we proposed to reduce the spatial discretization (mesh derefining) and we showed that this approach could sometimes lead to damp or move some phenomena (shocks). The second contribution of this work relies on evaluating the quality of the extrapolated solution and their validity domain, and on building links between local and macroscopic behavior. Finally, we proposed a method that allows the mass flow rate characteristic of a component to be determined from the calculation of the extrapolated solution issued from a limited number of reference points.