L'utilisation de stents actifs (DES) a révolutionné le traitement de l'athérosclérose. Le relargage contrôlé de médicaments anti-prolifératifs dans la paroi artérielle (PA) a permis de réduire fortement le taux de resténose intra-stent. Mais le risque de thromboses intra-stents tardives demeure un enjeu majeur des DES en partie lié au retard de cicatrisation de la PA endommagée lors de l'implantation. Cette thèse présente une méthode d'optimisation du design des DES afin d'inhiber la resténose sans affecter la cicatrisation. Pour quantifier la performance des différents designs, un modèle numérique décrivant l'écoulement sanguin et le transport de médicaments dans les artères stentées a été développé. Il prend en compte la structure multi-couches de la PA et les interactions du médicament avec les cellules. Un algorithme d'optimisation est couplé au modèle afin d'identifier les DES optimaux. L'optimisation du temps de relargage ainsi que de la concentration initiale du médicament dans le revêtement du DES ont un effet significatif sur la performance. Lorsque le médicament utilisé est le paclitaxel, les solutions optimales consistent à relarguer le produit à des concentrations nettement inférieures à celles des DES actuels soit pendant quelques heures, soit pendant une durée d'un an. Pour le sirolimus, un relargage lent est nécessaire. Les formes optimales des spires du DES sont toujours allongées mais profilées seulement lorsque le relargage est rapide. Ces résultats permettent d'expliquer en partie les performances des différents DES récents et révèlent un fort potentiel d'amélioration dans la conception des DES par rapport aux dispositifs commerciaux actuels.

L'aérodynamique latérale des véhicules automobiles suscite de nos jours de plus en plus d'intérêt de la part des constructeurs. L'automobiliste est en effet soumis quotidiennement à de forts courants d'air latéraux, que ce soit lors du dépassement d'un autre véhicule, ou alorsen passant dans un couloir de vent du à la topographie du terrain (passage devant un espace entre deux immeubles par exemple). Les efforts aérodynamiques mis en jeu dans ces situations peuvent provoquer des mouvements non désirés du véhicule, pouvant avoir des conséquences dramatiques si le conducteur se laisse surprendre. Des études expérimentales reproduisant les effets d'un dérapage dynamique ont mis en évidence des phénomènes transitoires importants mettant à défaut les modèles stationnaires généralement pratiqués par les constructeurs pour qualifier le comportement de leurs véhicules en présence de dérapage. Les mécanismes responsables de ces phénomènes transitoires sont encore mal connus de la communauté scientifique. Ce travail propose d'approfondir ce sujet au travers de l'étude de l'aérodynamique d'un véhicule terrestre fixe soumis à un vent longitudinal et à une rafale de vent latéral. Le but principal est d'identifier les structures tourbillonnaires au moyen de mesures PIV et de calculs numériques des champs de vitesse autour d'une maquette automobile et de les corréler aux efforts aérodynamiques. Un accord entre l'ISAT, composante de l'Université de Bourgogne, et l'Institut Supérieurde l'Aéronautique et de l'Espace (ISAE) de Toulouse a permis de mener l'étude avec les ressources de cet établissement. Le moyen d'essai principal, créé par l'ISAE, est le banc" rafale latérale ", constitué d'une soufflerie principale et d'une soufflerie secondaire, dont la sortie à volet déferlants (" Mexican Wave ") est inspirée de l'approche proposée par Ryan et Dominy (2000). L'analyse expérimentale a été effectuée à l'aide de la PIV résolue en temps et stéréoscopique, et d'une balance dard instationnaire à cinq composantes. Un banc" numérique " identique a été constitué à l'aide du logiciel FLUENT©, pour des calculs 3D. De plus, un modèle 2D annexe, basé sur la méthode " meshless ", a été développé pour de futures investigations, en raison de sa robustesse pour des problèmes à fortes discontinuités et sa bonne adaptabilité aux problèmes avec frontières mobiles.Une première phase de ce travail a consisté en la mise au point des bancs expérimental et numérique, avec génération d'un champ de dérapage homogène, de 21° dans la zone de mesure. L'évolution du dérapage en chaque point respecte bien la forme d'un créneau imposé par la rafale. Pour l'analyse des efforts, deux géométries de maquette ont été étudiées, à savoir un corps de Windsor à culot droit générant, pour un écoulement longitudinal, des structures de sillage bidimensionnelles, et son homologue à culot incliné de 25°, générant des tourbillons " cigare ". Des pics d'efforts ont été observés à l'arrivée de la rafale, tout comme la littérature le prédit. Pour ce qui est du coefficient du moment de lacet, les sursauts sont de 29 % et 19 % respectivement pour la maquette à culot droit et celle à culot incliné, par rapport aux valeurs stationnaires. Concernant le coefficient de force de dérive, ils sont de 10 % et 14 %, respectivement. Lors de nos essais, ces efforts se sont établis après 5.5 longueurs de maquette. Afin d'expliquer la différence de comportement entre les deux maquettes en termes d'efforts, l'évolution temporelle des tourbillons nommés, dans ce mémoire, ΓA, ΓB, ΓC et Γ1 à été discutée. Il en est ressorti une forte corrélation entre la circulation du tourbillon ΓA, le plusénergétique, naissant à l'avant du flanc sous le vent de la maquette, et les efforts latéraux, de sorte que ce tourbillon serait le meilleur témoin des phénomènes instationnaires mis en jeu dans l'étude de l'effet du vent latéral. [...]

Ce travail de thèse, réalisé dans le cadre d'une convention CIFRE entre TURBOMECA et le CERFACS et en partenariat avec l'IVK, se place dans un contexte d'amélioration des performances des turbines axiales équipant les turboréacteurs d'hélicoptère. Un des points critiques du dimensionnement de tels moteurs est la maitrise de la durée de vie des pales de la turbine haute pression qui font face à de très hautes températures provenant de la chambre de combustion. Les prédictions numériques de l'environnement aérothermique des pales (écoulements dans la veine et système de refroidissement) sont réalisées aujourd'hui dans le milieu industriel à l'aide de la modélisation Reynolds Averaged Navier-Stokes (RANS). Grâce à des capacités de calculs grandissantes, l'approche Simulation aux Grandes Echelles (SGE) offre désormais un nouveau potentiel de prédictions d'écoulements. Les travaux de cette thèse s'intéressent ainsi à la capacité de la SGE à prédire l'écoulement du circuit de refroidissement interne d'une pale de turbine. Pour simplifier l'analyse de ce problème ou plusieurs phénomènes physiques sont en jeu, une progression en trois parties est proposée. La première s'intéresse à l'étude aérothermique de géométries simplifiées de canaux de refroidissement (coude à 180° et canal avec promoteurs de turbulence) en configuration statique. Aux régimes d'écoulement considérés, une approche résolue en paroi avec maillage non-structuré hybride est proposée et validée en vue d'une application industrielle facilitée. La seconde partie étend l'analyse de l'écoulement à un cas de canal avec promoteurs de turbulence en rotation utilisant une méthode de résolution numérique dans un repère absolu. Les investigations des résultats de la SGE fournissent des prédictions moyennes et instationnaires en bon accord avec les expériences disponibles et les travaux précédents aussi bien pour la dynamique de l'écoulement que les transferts de chaleur. Enfin, une troisième partie présente une application de la méthode sur un cas de pale réelle avec couplage thermique entre le circuit de refroidissement et le solide de la pale. Cette dernière partie classée confidentielle n'est pas présente dans le manuscrit disponible publiquement. Les résultats de l'approche résolue en paroi et de la rotation dans le repère absolu comparés aux résultats RANS disponibles pour le cas applicatif montrent d'importante différences locales et ainsi le potentiel de la méthode proposée.

Mes activités de recherche se sont centrées sur les procédés de transformation des polymères, essentiellement le procédé d'injection-moulage mais aussi le soufflage et maintenant la peinture. J'ai participé à plusieurs projets de recherche appliquée pluridisciplinaire donc à finalité industrielle. Dans ces projets ma contribution porte sur la modélisation du procédé, que ce soit par l'utilisation de code de calculs ou le développement de modèles numériques à l'aide de logiciels libres. L'originalité de mon travail réside dans la pluridisciplinarité des projets qui mêlent des aspects théoriques et numériques aux aspects expérimentaux; la science des polymères à la thermique, à la mécanique et à la rhéologie. Les modélisations s'appuient sur des mesures expérimentales qui sont, dans les procédés de transformation, particulièrement délicats à isoler et à interpréter. Ces dispositifs expérimentaux représentent une part considérable de l'effort de recherche dans ces projets. Durant ces travaux, j'ai tenté de répondre à des problématiques distinctes : Dans les procédés de plasturgie, la viscoélasticité est finalement assez absente. Lorsqu'elle se manifeste, c'est sous la forme de défauts (instabilités, recirculations) qu'il faut éviter. C'est pourquoi, la plupart des écoulements impliqués dans les procédés sont de type lubrification~: une dimension au moins est petite devant les autres; ce qui atténue de façon effective les effets viscoélastiques. Une illustration frappante en est donnée par les recirculations viscoélastiques créées par la courbure du chenal d'une vis : lorsque le rapport d'aspect de la section augmente, elles s'atténuent et se décalent vers les parois latérales, précisément là où dans la section, aucune dimension n'est petite devant l'autre. Seul le procédé d'extrusion soufflage s'appuie totalement sur les propriétés viscoélastiques du polymère en écoulement, au travers des phénomènes de gonflement et de fluage. Cet écoulement est particulièrement discriminant et doit permettre à terme de déterminer les paramètres matériels d'une loi viscoélastique intégrale en comparant les prédictions des simulation numériques (par le code FEM Polyflow ou la méthode des tubes de courant) aux mesures expérimentales du gonflement et de fluage, devenues réalisables désormais avec notre technique de visualisation par illumination laser. En ce qui concerne le phénomène de plastification dans les monovis d'injection, ce thème de recherche est celui qui a nécessité le plus d'efforts et c'est celui qui est le plus ardu. La difficulté première, celle des observations, est en passe d'être surmontée avec le fourreau à fenêtres qui permet une visualisation inégalée des phénomènes de plastification. Nous allons pouvoir reprendre et vérifier toutes les hypothèses usuelles de la plastification : accélération du lit solide, fonte contiguë ou dispersée, position et épaisseur des films fluides. En ce qui concerne la modélisation numérique du convoyage solide du polymère dans le chenal de la vis. Une autre problématique, très importante dans le domaine industriel de l'automobile ou de la cosmétique, concerne les défauts de surfaces des pièces plastiques injectées. De nombreux facteurs peuvent contribuer à l'apparition de tels défauts. Nous avons étudié les effets thermiques en montrant que même des revêtements métalliques extrêmement fins pouvaient avoir une influence thermique. Qu'en est-il de l'aspect dynamique de l'écoulement et du contact (tension de surface) polymère-métal ? Ce sujet est d'importance pour l'injection moulage classique mais aussi la micro-plasturgie. La simulation numérique de la projection de peinture représente une direction nouvelle de recherche, active en Europe, mais inédite en France, et pour laquelle nous sommes bien placés à l'insa-lyon avec la création récente d'une cabine de peinture pilote. Cela reste toujours de la mécanique des fluides et des polymères mais avec des outils différents (code volumes finis, OpenFOAM, StarCCM). Encore une fois, une analyse des phénomènes prépondérants, notamment dans l'aérodynamique et le champ électrostatique, associée à des mesures expérimentales, devraient aboutir à un modèle simplifié des trajectoires des gouttes et donc de la projection de peinture.

Améliorer la gestion et l’organisation des blocs opératoires est une tâche critique dans les hôpitaux modernes, principalement à cause de la diversité et l’urgence des activités impliquées. Contrairement à l’aviation civile, qui a su optimiser organisation et sécurité, le management de bloc opératoire est plus délicat. Le travail ici présenté abouti au développement et à l’installation de nouvelles technologies assistées par ordinateur résolvant les problèmes quotidiens des blocs opératoires. La plupart des systèmes existants modélisent le flux chirurgical et sont utilisés seulement pour planifier. Ils sont basés sur des procédés stochastiques, n’ayant pas accès à des données sûres. Nous proposons une structure utilisant un modèle multi-agent qui comprend tous les éléments indispensables à une gestion efficace et au maintien de la sécurité dans les blocs opératoires, allant des compétences communicationnelles du staff, au temps nécessaire à la mise en place du service de nettoyage. Nous pensons que la multiplicité des ressources humaines engagées dans cette structure cause des difficultés dans les blocs opératoires et doit être prise en compte dans le modèle. En parallèle, nous avons construit un modèle mathématique de flux d’air entre les blocs opératoires pour suivre et simuler la qualité de l’environnement de travail. Trois points sont nécessaires pour la construction et le bon fonctionnement d’un ensemble de bloc opératoire : 1) avoir accès au statut du système en temps réel grâce au placement de capteurs 2) la construction de modèles multi-échelles qui lient tous les éléments impliqués et leurs infrastructures 3) une analyse minutieuse de la population de patients, du comportement des employés et des conditions environnementales. Nous avons développé un système robuste et invisible qui permet le suivi et la détection automatique d’événements dans les blocs. Avec ce système nous pouvons suivre l’activité à la porte d’entrée des blocs, puis l’avancement en temps réel de la chirurgie et enfin l’état général du bloc. Un modèle de simulation numérique de mécanique des fluides de plusieurs blocs opératoires est utilisé pour suivre la dispersion de fumée chirurgicale toxique, ainsi qu’un modèle multi-domaine qui évalue les risques de propagation de maladie nosocomiale entre les blocs. La combinaison de ces trois aspects amène une nouvelle dimension de sensibilisation à l’environnent des blocs opératoires et donne au staff un système cyber-physique capable de prédire des événements rares impactant la qualité, l’efficacité, la rentabilité et la sécurité dans l’hôpital. / Improving operating room management is a constant issue for modern large hospital systems who have to deal with the reality of day to day clinical activity. As opposed to other industrial sectors such as air civil aviation that have mastered the topic of industry organization and safety, progress in surgical flow management has been slower. The goal of the work presented here is to develop and implement technologies that leverage the principles of computational science to the application of OR suite problems. Most of the currently available models of surgical flow are used for planning purposes and are essentially stochastic processes due to uncertainties in the available data. We propose an agent-based model framework that can incorporate all the elements, from communication skills of the staff to the time it takes for the janitorial team to go clean an OR. We believe that human factor is at the center of the difficulty of OR suite management and should be incorporated in the model. In parallel, we use a numerical model of airflow at the OR suite level to monitor and simulate environment conditions inside the OR. We hypothesize that the following three key ingredients will provide the level of accuracy needed to improve OR management : 1) Real time updates of the model with ad hoc sensors of tasks/stages 2) Construction of a multi-scale model that links all key elements of the complex surgical infrastructure 3) Careful analysis of patient population factors, staff behavior, and environment conditions. We have developed a robust and non-obtrusive automatic event tracking system to make our model realistic to clinical conditions. Not only we track traffic through the door and the air quality inside the OR, we can also detect standard events in the surgical process. We propose a computational fluid dynamics model of a part of an OR suite to track dispersion of toxic surgical smoke and build in parallel a multidomain model of potential nosocomial contaminant particles flow in an OR suite. Combining the three models will raise the awareness of the OR suite by bringing to the surgical staff a cyber-physical system capable of prediction of rare events in the workflow and the safety conditions.

Cette thèse exploratoire pose les bases scientifiques et méthodologiques d’une approche transversale visant à étudier l’énergétique urbaine et le bio-climatisme. Elle fait appel à des concepts et des outils de l’architecture et l’urbanisme, et à la physique du bâtiment et de la ville. Cette thèse étudie les relations entre la morphologie urbaine et les processus aérodynamiques qui se développent dans la canopée urbaine et leurs effets sur la demande énergétique des bâtiments induite par les infiltrations d’air et les échanges thermiques convectifs. Les spécificités de l’aérodynamique et de la physique urbaines sont d’abord synthétisées et la morphologie de tissus urbains réels est analysée. Une typologie générique de bâtiments isolés et une autre d’îlots urbains en sont déduites. Le modèle CFD est ensuite validé par comparaison des prédictions du modèle avec des résultats expérimentaux et numériques, et des expérimentations numériques sont réalisées sur les différents types morphologiques. Les écoulements moyens sont analysés dans leurs rapports avec la morphologie bâtie, et la distribution des coefficients de pression sur les façades des bâtiments est analysée. Ensuite, les échanges thermiques sont couplés aux processus aérodynamiques. L’amélioration des estimations des échanges convectifs des bâtiments grâce à la CFD est vérifiée par comparaison des résultats de simulation avec des données expérimentales et numériques, ainsi qu’avec les valeurs standard. Une adaptation des fonctions de paroi relatives au transfert thermique est proposée sur la base d’études existantes, et la distribution des échanges convectifs sur les façades de bâtiments est analysée. Enfin, la demande énergétique des bâtiments due aux infiltrations d’air et à la transmission de chaleur au travers de leur envelope est estimée pour différents types morphologiques, et comparée avec les valeurs estimées suivant une approche réglementaire. Les résultats de cette thèse mettent en évidence les effets des propriétés topologiques et métriques des bâtiments et ensembles bâtis sur le développement de recirculations d’air dans la canopée urbaine. Celles-ci induisent une distribution et intensité hétérogènes des coefficients de pression et d’échange convectif sur les façades des bâtiments, qui influent sur le comportement thermique des bâtiments non isolés et perméables à l’air. Par ailleurs, l’estimation de leur demande énergétique diffère suivant si celle-ci est basée sur les valeurs simulées ou standard des coefficients de pression et d’échange convectif. Cependant, l’influence relative de la structure bâtie sur la demande énergétique des bâtiments apparaît plus importante pour les bâtiments isolés thermiquement. La différence entre la demande énergétique par unité de surface de plancher, due aux infiltrations d’air et pertes thermiques au travers de l’enveloppe peut varier de 18% à 47% suivant si le bâtiment est isolé ou situé dans un environnement bâti. / This thesis is an exploratory study that lays the scientific and methodological foundations of a transverse approach for studying urban energy and bio-climatic issues. This approach involves concepts and tools of building and urban physics as well as urban planning and architecture. It addresses the relations between urban morphology and aerodynamic processes, and studies their effects on the building energy loads due to infiltration and convective heat losses. This thesis is divided into three main parts. The first part synthesizes the specificities of urban aerodynamics and urban physics, and analyzes existing urban fabrics from a morphological point of view. Generic typologies of isolated buildings and urban blocks for small scale aerodynamic studies are deduced. The second part validates the computational fluid dynamics (CFD) model (steady RANS RSM) against detailed experimental and numerical data, and presents the numerical experiments performed on the different morphological types. Mean flow structures that develop according to the construction shape and built environment, as well as pressure distribution on the building outer walls are examined. The last part couples heat and air fluxes to evaluate the contribution of urban air flows on the building energy loads. The improvement brought by CFD to the assessment of building convective heat transfers is verified by comparing numerical results to experimental data, detailed numerical studies and standard correlations. An enhanced temperature wall-function adapted for forced convection problems is adjusted to the model based on existing studies, and the convective heat transfers distribution on building facades is analyzed. Finally, the building energy loads due to air infiltration and heat transmission are estimated for typical constructions and compared to standard values. The results of this thesis show strong effects of the topology and dimensionality of constructions and urban structures on the development of recirculation phenomena within the urban canopy layer. The related aerodynamic conditions yield heterogeneous pressure and convective heat transfer intensities and distributions on building facades, which depend upon the considered built morphology. Their effects on building energy loads are logically particularly important in absolute value for buildings that are neither insulated nor air tight. Nonetheless, the estimates of the building energy needs based on standard or simulated pressure and convective heat transfer coefficients often show substantial deviation. Focusing on the relative contribution of the built structure, the effects of the aerodynamic context appear more influential for insulated buildings. Essentially, switching from an exposed to a sheltered building may decrease the energy needs per surface unit of floor due to air infiltration and heat transmission through outer walls by 18% up to 47% according to simulation.

L'ablation par radiofréquence (ARF) de tumeurs abdominales est rendue difficile par l’influence des vaisseaux sanguins et les variations de la conductivité thermique, compliquant la planification spécifique à un patient donné. En fournissant des outils prédictifs, les modèles biophysiques pourraient aider les cliniciens à planifier et guider efficacement la procédure. Nous introduisons un modèle mathématique détaillé des mécanismes impliqués dans l’ARF des tumeurs du foie comme la diffusion de la chaleur et la nécrose cellulaire. Il simule l’étendue de l’ablation à partir d’images médicales, d’après lesquelles des modèles personnalisés du foie, des vaisseaux visibles et des tumeurs sont segmentés. Dans cette thèse, une nouvelle approche pour résoudre ces équations basée sur la méthode de Lattice Boltzmann est introduite. Le modèle est d’abord évalué sur des données de patients qui ont subi une ARF de tumeurs du foie. Ensuite, un protocole expérimental combinant des images multi-modales, anatomiques et fonctionnelles pré- et post-opératoires, ainsi que le suivi de la température et de la puissance délivrée pendant l'intervention est présenté. Il permet une validation totale du modèle qui considère des données les plus complètes possibles. Enfin, nous estimons automatiquement des paramètres personnalisés pour mieux prédire l'étendu de l’ablation. Cette stratégie a été validée sur 7 ablations dans 3 cas cliniques. A partir de l'étude préclinique, la personnalisation est améliorée en comparant les simulations avec les mesures faites durant la procédure. Ces contributions ont abouti à des résultats prometteurs, et ouvrent de nouvelles perspectives pour planifier et guider l’ARF. / The outcome of radiofrequency ablation (RFA) of abdominal tumors is challenged by the presence of blood vessels and time-varying thermal conductivity, which make patient-specific planning extremely difficult. By providing predictive tools, biophysical models may help clinicians to plan and guide the procedure for an effective treatment. We introduce a detailed computational model of the biophysical mechanisms involved in RFA of hepatic tumors such as heat diffusion and cellular necrosis. It simulates the extent of ablated tissue based on medical images, from which patient-specific models of the liver, visible vessels and tumors are segmented. In this thesis, a new approach for solving these partial differential equations based on the Lattice Boltzmann Method is introduced. The model is first evaluated against clinical data of patients who underwent RFA of liver tumors. Then, a comprehensive pre-clinical experiment that combines multi-modal, pre- and post-operative anatomical and functional images, as well as the interventional monitoring of the temperature and delivered power is presented. This enables an end-to-end validation framework that considers the most comprehensive data set for model validation. Then, we automatically estimate patient-specific parameters to better predict the ablated tissue. This personalization strategy has been validated on 7 ablations from 3 clinical cases. From the pre-clinical study, we can go further in the personalization by comparing the simulated temperature and delivered power with the actual measurements during the procedure. These contributions have led to promising results, and open new perspectives in RFA guidance and planning.

De récentes études expérimentales ont montré que les valeurs usuelles du coefficient d’échange convectif sont différentes en présence de matériaux à changement de phase. Cette thèse de doctorat porte sur l'étude numérique des échanges convectifs fluide/paroi dans une cavité ouverte en régime dynamique. Plus précisément, les parois étudiées sont une paroi avec une capacité thermique et une paroi qui contient des matériaux à changement de phase. Trois modèles distincts ont été développés. Dans un premier temps un modèle (modèle 1) qui concerne l’interaction fluide-paroi à la surface d’une paroi résistive (temperature imposée) en régime laminaire stationnaire a été développé et validé. Les résultats ont été confrontés avec la littérature. Ensuite, les échanges convectifs à la surface d’une plaque capacitive (modèle 2) soumise à une rampe de température d’air ont été étudiés. Finalement, un troisième modèle (modèle 3) a été développé, à la suite du modèle 2. Ce dernier modèle concerne l’interaction fluide-paroi à la surface d’une paroi contenant des matériaux à changement de phase en régime dynamique. Les résultats obtenus révèlent des pics locaux du flux de chaleur au cours du temps. Ce fait témoigne du changement d’état à l’intérieur de la paroi qui contient le materiau à changement de phase. De plus, les courbes des coefficients d’échanges convectifs moyens révèlent la dépendance du coefficient d’echange convectif à la capacité thermique du materiau. Par conséquent, la présence des matériaux à changement de phase à l’intérieur d'une paroi influence l’évolution et la forme de la couche limite thermique. / Recent experimental studies have shown that the usual values ​​of the convective heat transfer coefficient h are no longer valid in the presence of phase change materials. Three separate models were developed. Initially a model 1 which treats the fluid-wall (constant temperature) interaction in steady laminar flow has been developed and validated. Then, the wall with heat capacity (model 2) subjected to an air temperature ramp were studied. Finally, a third model (3) has been developed which treats the interaction fluid-wall which contains a phase change material. The results show local peaks of heat flow over time. This fact reflects the phase change inside the wall. Moreover, the curves of the convective heat transfer coefficient indicate the dependence of the coefficient h to the wall’s energy storage capacity. Therefore, the presence of the phase change materials within a wall effect and changes the shape of the thermal boundary layer.

<p align="justify">Une approche multidimensionelle pour la résolution numérique des équations d'Euler et de Navier-Stokes sur maillages non-structurés est proposée. Dans une première partie, un exposé complet des schémas de distribution, dits de "fluctuation-splitting" ,est décrit, comprenant une étude comparative des schémas décentrés, positifs et de 2ème ordre, pour résoudre l'équation de convection à coefficients constants, ainsi qu'une étude théorique et numérique de la précision des schémas sur maillages réguliers et distordus. L'extension à des lois de conservation non-linéaires est aussi abordée, et une attention particulière est portée au problème de la linéarisation conservative. Dans une deuxième partie, diverses discrétisations des termes visqueux pour l'équation de convection-diffusion sont développées, avec pour but de déterminer l'approche qui offre le meilleur compromis entre précision et coût. L'extension de la méthode aux systèmes des lois de conservation, et en particulier à celui des équations d'Euler de la dynamique des gaz, représente le noyau principal de la thèse, et est abordée dans la troisième partie. Contrairement aux schémas de distribution classiques, qui reposent sur une extension formelle du cas scalaire, l'approche développée ici repose sur une décomposition du résidu par élément en équations scalaires, modélisant le transport de variables caracteristiques. La difficulté vient du fait que les équations d'Euler instationnaires ne se diagonalisent pas, et admettent une infinité de solutions élémentaires (ondes simples) se propageant dans toutes les directions d'espace. En régime stationnaire, en revanche, les équations se diagonalisent complètement dans le cas des écoulements supersoniques, et partiellement dans le cas des écoulements subsoniques. Ainsi, les équations sous forme conservative peuvent être remplacées par un système équivalent comprenant deux équations totalement découplées, exprimant l'invariance de l'entropie et de l'enthalpie totale le long des lignes de courant, et deux autres équations, modélisant les effets purement acoustiques. En régime supersonique, celles-ci se découplent aussi, et expriment la convection le long des lignes de Mach d'invariants de Riemann généralisés. La discrétisation de ces équations par des schémas scalaires décentrés permet de simuler des écoulements continus et discontinus avec une grande précision et sans oscillations. Finalement, dans une dernière partie, l'extension aux équations de Navier-Stokes est abordée, et la discrétisation des termes visqueux par une approche éléments finis est proposée. Les résultats numériques confirment la précision et la robustesse de la méthode.</p> / Doctorat en sciences appliquées / info:eu-repo/semantics/nonPublished

L’aérodynamique latérale des véhicules automobiles suscite de nos jours de plus en plus d’intérêt de la part des constructeurs. L’automobiliste est en effet soumis quotidiennement à de forts courants d’air latéraux, que ce soit lors du dépassement d’un autre véhicule, ou alorsen passant dans un couloir de vent du à la topographie du terrain (passage devant un espace entre deux immeubles par exemple). Les efforts aérodynamiques mis en jeu dans ces situations peuvent provoquer des mouvements non désirés du véhicule, pouvant avoir des conséquences dramatiques si le conducteur se laisse surprendre. Des études expérimentales reproduisant les effets d’un dérapage dynamique ont mis en évidence des phénomènes transitoires importants mettant à défaut les modèles stationnaires généralement pratiqués par les constructeurs pour qualifier le comportement de leurs véhicules en présence de dérapage. Les mécanismes responsables de ces phénomènes transitoires sont encore mal connus de la communauté scientifique. Ce travail propose d’approfondir ce sujet au travers de l’étude de l’aérodynamique d’un véhicule terrestre fixe soumis à un vent longitudinal et à une rafale de vent latéral. Le but principal est d’identifier les structures tourbillonnaires au moyen de mesures PIV et de calculs numériques des champs de vitesse autour d’une maquette automobile et de les corréler aux efforts aérodynamiques. Un accord entre l’ISAT, composante de l’Université de Bourgogne, et l’Institut Supérieurde l’Aéronautique et de l’Espace (ISAE) de Toulouse a permis de mener l’étude avec les ressources de cet établissement. Le moyen d’essai principal, créé par l’ISAE, est le banc« rafale latérale », constitué d’une soufflerie principale et d’une soufflerie secondaire, dont la sortie à volet déferlants (« Mexican Wave ») est inspirée de l’approche proposée par Ryan et Dominy (2000). L’analyse expérimentale a été effectuée à l’aide de la PIV résolue en temps et stéréoscopique, et d’une balance dard instationnaire à cinq composantes. Un banc« numérique » identique a été constitué à l’aide du logiciel FLUENT©, pour des calculs 3D. De plus, un modèle 2D annexe, basé sur la méthode « meshless », a été développé pour de futures investigations, en raison de sa robustesse pour des problèmes à fortes discontinuités et sa bonne adaptabilité aux problèmes avec frontières mobiles.Une première phase de ce travail a consisté en la mise au point des bancs expérimental et numérique, avec génération d’un champ de dérapage homogène, de 21° dans la zone de mesure. L’évolution du dérapage en chaque point respecte bien la forme d’un créneau imposé par la rafale. Pour l’analyse des efforts, deux géométries de maquette ont été étudiées, à savoir un corps de Windsor à culot droit générant, pour un écoulement longitudinal, des structures de sillage bidimensionnelles, et son homologue à culot incliné de 25°, générant des tourbillons « cigare ». Des pics d’efforts ont été observés à l’arrivée de la rafale, tout comme la littérature le prédit. Pour ce qui est du coefficient du moment de lacet, les sursauts sont de 29 % et 19 % respectivement pour la maquette à culot droit et celle à culot incliné, par rapport aux valeurs stationnaires. Concernant le coefficient de force de dérive, ils sont de 10 % et 14 %, respectivement. Lors de nos essais, ces efforts se sont établis après 5.5 longueurs de maquette. Afin d’expliquer la différence de comportement entre les deux maquettes en termes d’efforts, l’évolution temporelle des tourbillons nommés, dans ce mémoire, ΓA, ΓB, ΓC et Γ1 à été discutée. Il en est ressorti une forte corrélation entre la circulation du tourbillon ΓA, le plusénergétique, naissant à l’avant du flanc sous le vent de la maquette, et les efforts latéraux, de sorte que ce tourbillon serait le meilleur témoin des phénomènes instationnaires mis en jeu dans l’étude de l’effet du vent latéral. [...] / The automotive manufacturers are nowadays more and more interested in crosswind aerodynamics. Indeed, the driver is subjected every day to strong side air flows, for example when overtaking another vehicle or when passing through a lateral wind wall, generated by terrain topography (as in the case of the passage near the empty space between two buildings).The aerodynamic efforts generated in these situations can lead to undesired lateral deviations,which can be dramatic if the driver is surprised. Different experimental studies, reproducing the effects of a dynamic yaw angle, pointed out the issues of the steady methods, commonly used to qualify the vehicle crosswind behaviour. Little is still known about the physics behind these unsteady phenomena. This is the main topic of this work, by studying the aerodynamics of a fixed vehicle subjected to both a longitudinal flow and a side wind gust. The goal is the identification of the near-vehicle vortex structures, by means of PIV measurements and CFD calculations, and their correlation with the evolution of the efforts. An agreement between the ISAT, a department of the University of Burgundy, and the ISAE of Toulouse, permitted to carry out this research with the resources of the latter laboratory. The work focuses on the use of the “rafale latérale” (side gust) test bench, made up with a main wind tunnel connected with an auxiliary one by means of a shutter system,whose opening is held by a “Mexican Wave” law. This approach is inspired by the work of Dominy and Ryan (2000). The experimental analysis was carried out by means of Time-Resolved and stereoscopic PIV, and by a five components unsteady balance as well. Anidentical test bench was numerically reproduced with the 3D CFD software FLUENT©.Moreover, an additional 2D CFD model, based on the meshless method, has been developed for future studies. This kind of approximation method has been chosen for its robustness innon-continuous problems and because of its adaptability when moving boundaries are needed.The first phase of this work consisted on wind tunnels set-up, both for the real test bench and for the CFD model. The yaw angle field is homogeneous, 21° in the measurement region. The yaw angle evolution, at given point, respects the step wise behaviour, imposed by the gust passage. As far as the efforts are concerned, two versions of the Windsor body car model were studied, that is a squareback geometry, generating, for longitudinal flows, 2D wakestructures, and a fastback geometry (rear window inclined by 25°), producing cone-liketrailing vortices. Force overshoots were seen after the gust arrival, as seen in literature. In particular, the yaw moment coefficient overshoots are 29% and 19% higher than the steady yaw angle tests, for the squareback and the fastback geometries, respectively. If the side forceis concerned, the entities of these overshoots are 10% and 14%, respectively. Our testspointed out that efforts establish after the vehicle has driven 5.5 times its length in thecrosswind. In order to explain the different behaviour of the two geometries, it is discussed about the unsteady evolution of the vortices called, in this report, ΓA, ΓB, ΓC et Γ1. A strong correlation between the side efforts and the circulation of the most energetic vortex, ΓA,having its origin in the front leeward side of the vehicle. The ΓA vortex is so the best index for the study of the crosswind unsteady phenomena. The coupled analysis between vortex structures and efforts confirmed the presence of a higher side force for the squareback geometry. The inverted effect has been observed for the yaw moment