En tribologie, la modélisation numérique est aujourd'hui un outil indispensable pour étudier un contact afin de pallier les limites expérimentales. Pour comprendre de mieux en mieux les phénomènes mis en jeu, les modèles ne se situent plus à une seule échelle, mais en font intervenir plusieurs, rendant plus que jamais le concept de triplet tribologique incontournable. Travaillant avec cette philosophie et en se basant sur l'approche Non Smooth Contact Dynamics, dont nous rappelons les grandes lignes, nous proposons de franchir deux cas: proposer des modèles offrant des résultats quantitatifs et mettre en place les premières pièces d'une homogénéisation au niveau du contact (VER). Dans le premier cas, l'étude du couplage éléments finis/éléments discrets au sein d'une même simulation a pour but de proposer des modèles plus "réalistes". Même si l'interface utilisée est déjà présente au coeur du contact et ne va pas évoluer, elle permet de mettre en évidence l'utilisation d'outil de mesure permettant de lier le mouvement des particules aux instabilités dynamiques et permet d'avoir des résultats qualitatifs mais aussi quantitatifs puisque la comparaison avec les taux de contraintes expérimentaux sont en très bonne adéquation. Dans le second cas, le VER sous sollicitations tribologiques est étudié afin d'étendre les techniques d'homogénéisation aux problèmes de contact afin de s'affranchir de la description des interfaces aux grandes échelles en trouvant un moyen d'homogénéiser le comportement hétérogène de l'interface et de le faire dialoguer avec le comportement continu des corps en contact en faisant remonter, dans un sens, des grandeurs moyennées à l'échelle microscopique à l'échelle macroscopique des premiers corps et dans l'autre sens, se servir des données locales à l'échelle macroscopique comme conditions limites à l'échelle microscopique.

Nous constatons que les universités sont confrontées à un double mouvement de la société sur leur organisation provoquant une nécessité de changement. Le premier concerne l’environnement institutionnel changeant depuis la loi Enseignement Supérieur et de la Recherche de 2013 : passage à la responsabilité et compétences élargies, contraintes budgétaires fortes, évolutions sociales des apprenants, nouvelle pratique managériale dans le secteur public. Ce premier mouvement provoque pour nous un paradoxe : la naissance d’un phénomène d’isomorphisme institutionnel, les comportements des universités convergent quand bien même chacune d’elles tente d’être plus attractive que l’autre sur le marché de la formation tout au long de la vie. Par ailleurs, un second mouvement, qui concerne les membres et les usagers des universités, apparaît au regard des usages des Technologies de l'information et de la communication (TIC) dans la société : la facilité d’accès à la connaissance au moyen de nouvelles plateformes, l’adoption des TIC dans la société catalysent cela. En réaction à ces deux mouvements, nous estimons que les universités sont en changement. Et nous souhaitons étudier en quoi l’évolution des usages du numérique conduit à la nécessité d’une gouvernance du numérique universitaire. Nous tenterons de démontrer que l’instabilité provoquée par l’arrivée du numérique dans les universités est une occasion de faire évoluer le système universitaire, au point de créer un nouveau point d’équilibre s’appuyant sur une nouvelle gouvernance : la gouvernance du numérique universitaire. Celle-ci permettrait de développer les usages du numérique de tous les acteurs, et placerait ainsi l’établissement dans un meilleur positionnement dans un environnement devenu de plus en concurrentiel, tout en améliorant son fonctionnement. Pour cela nous construirons les analyses à partir de plusieurs terrains complémentaires, non pas dans un souci de comparaison de terrains, mais bien de vérification partielle d’une ou plusieurs hypothèses au sein de chacun de ces 5 terrains d’observation. Sur la base des analyses de ces expériences et de questionnement sur cette notion de gouvernance du numérique universitaire, cette thèse souhaite contribuer à la réflexion sur les gouvernances d’universités, afin qu’elles puissent mieux appréhender l’éducation supérieure du XXIe siècle. / We see that universities are confronted with a double movement of society on their organization, necessitating a change. The first concerns the changing institutional environment since the Law Higher Education and Research in 2013 : shift to responsibility and broadened competencies, strong budgetary constraints, social evolutions of learners, new managerial practice in the public sector. This first movement provokes for us a paradox: the birth of a phenomenon of institutional isomorphism, the behavior of universities converge even though each of them tries to be more attractive than the other on the market of training all the way long life. In addition, a second movement, involving the players and users of universities, appears in relation to the uses of Information and Communication Technologies (ICT) in society: the ease of access to knowledge through new platforms, the adoption of ICTs in society catalyze this. In response to these two movements, we believe that universities are changing. And we want to study how the evolution of digital uses leads to the need for digital university governance. We will try to show that the instability brought about by the arrival of digital in the universities is an opportunity to evolve the university system to the point of creating a new equilibrium point based on a new governance : digital university governance . This would allow the development of the digital uses of all the players, thus placing the establishment in a better position in a more competitive environment, while improving its functioning. For this we will construct the analyzes from several complementary fields, not for the sake of comparing land, but rather a partial verification of one or more hypotheses within each of these five observations. On the basis of the analyzes of these experiences and the questioning of this notion of digital university governance, this thesis wishes to contribute to the reflection on university governance so that they can better understand higher education in the 21st century.

La modélisation mathématique est un outil largement employé dans plusieurs disciplines des sciences appliquées. En hydrologie, en biologie, en économie ainsi que d'autres domaines des sciences naturelles, sociales et humaines, le recours à la modélisation mathématique est une démarche de plus en plus fréquente. Par exemple, en hydrologie, plusieurs modèles mathématiques sont conçus pour décrire ou prédire la relation existante entre les hauteurs d'eau et les débits des rivières. Dans le cadre de cette thèse nous nous sommes intéressés au développement de nouveaux modèles permettant de modéliser les phénomènes de croissance qui nécessitent la présence d'une asymptote linéaire croissante ou curviligne. Pour atteindre cet objectif, l'idée de base a été d'utiliser quelques modèles parmi les plus répandus en pratique et de les modifier judicieusement (et simplement) de façon à introduire une asymptote soit linéaire soit curviligne tout en conservant leur unique point d'inflexion. La modification que nous avons introduite conserve aussi le caractère simple et continue de ces modèles ainsi que la forme lisse et croissante de leurs courbes. Nous obtenons ainsi des modèles qui répondent aux besoins de la modélisation lorsque les modèles standards échouent.

Ce mémoire de maîtrise en milieu industriel présente les travaux réalisés chez Maya Heat Transfer Technologies afin d'implémenter la technique de la simulation des grandes échelles dans le solveur de dynamique des fluides parallèle NIECE. Les modèles de sous-maille [sic] algébriques de Smagorinsky, WALE (Wale Adapting Local Eddy-Viscosity) et le modèle de Vreman ont été retenus et validés.Ce projet de maîtrise se divise [entre] quatre volets. La première partie présente l'implémentation et l'analyse d'un nouveau schéma de discrétisation spatial et temporel dans le contexte d'un algorithme de résolution des équations de Navier-Stokes par une approche mixte [de] volumes éléments finis sur maillage non structuré. Une étude comparative sur la stabilité, la diffusion et [la] dipersion [sic] numérique des nouveaux schémas et ceux déjà existants dans le code permet d'identifier le schéma centré de second ordre en espace, couplé au schéma de Crank-Nicolson en temps comme un compromis adéquat pour la simulation des grandes échelles dans un contexte industriel. Le second volet détaille les résulats [i.e. résultats] de la simulation de la décroissance d'une turbulence homogène isotrope. La décroissance des spectres d'énergie est présentée pour chacun des modèles de sous-maille, puis comparée à la solution sans modèle. L'effet du raffinement du maillage est investigué.Ce volet détaille également l'implémentation et la validation de l'algorithme permettant l'initialisation d'une turbulence isotrope incompressible. Une discussion sur les formulations analytiques des spectres d'énergie de Passot-Pouquet et de von Karman-Pao est finalement présentée. Le problème de la répresentation [i.e. représentation] de l'écoulement à l'entrée du domaine de calcul est ensuite abordé.Ce troisième volet détaille la méthode de génération de turbulence, appelée méthode des tourbillons synthétiques, servant à générer un champ de vitesse cohérent et représentatif d'une turbulence réelle. Une étude de l'influence des paramètres numériques de la méthode permet de dériver un critère automatique de sélection pour l'utilisateur. Pour terminer, une étude comparative de la méthode des tourbillons synthétiques à une méthode d'injection aléatoire est présentée. En dernier lieu, la réalisation du cas de validation du canal plan périodique permet d'investiguer la capacité du code à simuler les propriétés d'une turbulence cisaillée. Les résultats des modèles de sous-maille sont comparés au cas sans modèle de turbulence. L'effet du raffinement du maillage est étudié sur trois différentes grilles. Les statistiques de premier ordre (profil de vitesse moyenne) et les satistiques [i.e. statistiques] de second ordre (profils des composantes du tenseur de Reynolds) sont détaillées et analysées. Finalement, la contribution de l'énergie cinétique résolue et de l'énergie de sous-maille permet de quantifier la contribution des modèles de turbulence sur la solution.

L’amélioration des sols par la technique de vibrosubstitution est habituellement utilisée pour les sols cohésifs, mais peut toutefois être utilisée dans les sols pulvérulents avec potentiel de liquéfaction. Cette technique permet la création d’un réseau de colonnes ballastées (CB) qui réduit le tassement et augmente la capacité portante d’un sol. L'intérêt est certes bien affiché. En effet, Geopac Inc. (entreprise spécialisée) se propose d'évaluer sa méthode d'amélioration des sols par CB en vue de réduire (voire éliminer) le potentiel de liquéfaction des sols grenus. Geopac a fait recours à une collaboration avec l'UdeS pour analyser expérimentalement et numériquement l'influence de la mise en place des CB sur les propriétés d'un sol pulvérulent. Un programme expérimental a été élaboré dans le but d’étudier l’effet de la mise en place des CB sur le sol encaissant par cette technique. À cet effet, un modèle réduit unique a été conçu, instrumenté et calibré à l’UdeS pour permettre dans un premier temps, de simuler l’incorporation de pierre nette (ballast) dans un sol lâche non saturé tout en maintenant une contrainte verticale et en mesurant la contrainte latérale générée durant l’incorporation. Le modèle réduit permet l’étude de l’expansion de cavité telle qu’obtenue durant la construction des CB par vibrosubstitution. Les colonnes sont construites à l’intérieur d’une cellule de PVC de 600 mm de haut et de 382 mm de diamètre. Un total de trois colonnes ballastées ont été construites sous une contrainte de 60, 80 et 100 kPa respectivement représentant un taux d’incorporation de près de 4 %. L’incorporation démontre une augmentation de la pression latérale caractérisée par le coefficient de pression latérale K jusqu’à 1.37 et une augmentation de la densité de 44 % à 55 % en moyenne. Une modélisation numérique a ensuite permis de confirmer les résultats expérimentaux et de vérifier le comportement du montage expérimental. La modélisation numérique a été réalisée avec le logiciel PLAXIS 2D 2012 en axisymétrie de révolution et considérant une loi comportementale de type Hardening Soil Model pour le sol encaissant et de type Mohr-Coulomb pour la colonne de pierre. La modélisation numérique a permis de déterminer que la forme de la contrainte verticale sur le sol est parabolique, c’est-à-dire maximale au centre et nulle à l’interface sol-paroi. L’expansion du sol a été simulée par l’imposition d’une déformation volumique anisotropique en cinq couches dans le sens radial et circonférentiel équivalent aux dimensions de la colonne obtenue expérimentalement. Les résultats numériques ont démontré une erreur de < 1 % pour la CB #1, de pour la CB #2 et de < 30 % pour la CB #3. L’erreur rencontrée sur la dernière colonne effectuée sous une contrainte verticale de 60 kPa peut-être associés à une mesure erronée lors de l'expérimentation. L’obtention des résultats expérimentaux par modélisation numérique a confirmé la fonctionnalité du montage expérimental et de la procédure de construction des CB en laboratoire. Cela ouvre la porte à des études subséquentes de l’effet de la mise en place des CB sur le sol encaissant en incorporant la vibration et la saturation du sable reproduisant la technique de vibrosubstitution.

La compréhension et la modélisation des phénomènes physiques qui sont mis en jeu au cours du procédé d'extrusion monovis de mélange d'élastomère constituent l'objectif du travail de recherche. Après formulation d'un mélange clair, celui-ci a fait l'objet de plusieurs campagnes d'essais rhéologiques. Le même mélange a été rendu plus glissant par l'ajout de stéarate de zinc. Nous avons pu montrer que la rhéologie de volume restait proche de celle du mélange précédent et proposer différentes lois de glissement, dépendantes de la géométrie, par l'intermédiaire des dépouillements classiques de la littérature. Une étude expérimentale systématique a été entreprise, tout d'abord, sur une extrudeuse industrielle (5 vis différentes représentatives des vis utilisées en production), puis sur une extrudeuse de laboratoire (profils simples et modulaires). Des informations ont été obtenues par le relevé de paramètres tels que la pression, la température, le débit, le taux de remplissage (démontage)... Ainsi, à partir des résultats expérimentaux et de leurs interprétations, nous avons pu comparer les performances des différentes géométries et mélanges, et extrapoler les résultats d'une machine à l'autre. L'approche de modélisation 1D a mis en évidence la grande importance des effets de bord, liée au caractère très pseudoplastique du produit utilisé. En conséquence, le modèle 1D est inutilisable. On s'est alors dirigé vers un calcul en deux dimensions par éléments finis. Il intègre différents couplages qui ont été réalisés et validés au fur et à mesure. La comparaison à l'expérience montre un certain nombre de désaccords. La comparaison entre un vrai calcul 3D et le calcul à deux dimensions couplé montre un écart grandissant lorsque le coefficient de pseudoplasticité diminue. Comme pour l'approche 1D, le couplage longitudinal/transversal est pris en défaut pour des caoutchoucs à cause d'un indice de pseudoplasticité très faible. L'avenir semble donc passer par un vrai calcul tridimensionnel. Nous en avons montré la faisabilité en simulant l'écoulement dans de petits tronçons de vis. L'accord avec l'expérience est satisfaisant.

Ce travail a concerné les anodes tournantes légères utilisées pour la production de rayons X dans les scanners médicaux. Ces anodes sont constituées de graphite revêtu de tungstène par projection plasma à basse pression. Nous nous sommes intéressés au comportement mécanique des matériaux constituants (particulièrement le tungstène) dans les conditions de sollicitations thermomécaniques intenses rencontrées durant la vie en service des anodes. La première étape de ce travail a été une phase d'identification du problème posé. Nous avons ainsi cerné les conditions de sollicitation en service en termes de fréquence d'excitation, de température, de déformation et de vitesse de déformation. Des expertises d'anodes usagées ont permis de dégager les principaux modes d'endommagement en service. Pour caractériser le comportement mécanique des matériaux, un dispositif d'essai existant a été adapté pour pouvoir travailler sur des éprouvettes de flexion de petite taille. Les procédures d'essai ont, quant à elles, été choisies pour générer des sollicitations proches de celles associées au fonctionnement. Dans toute la plage de température investiguée (ambiante à 1800°C), notre graphite s'est avéré posséder un comportement élastique fragile. Le comportement du tungstène s'est révélé plus complexe: élastique-fragile jusqu'à 400°C, il devient ensuite plastique pour devenir sensible au fluage à partir de 1200°C. Des trajets de chargement originaux ont permis de mettre en évidence l'existence d'une contrainte interne de rappel, ainsi qu'un couplage entre déformations plastiques et déformations visqueuses. Une loi phénoménologique originale de comportement, à double déformation inélastique, plastique et viscoplastique avec un terme d'interaction entre les deux mécanismes d'écoulement, a donc été nécessaire pour décrire le comportement mécanique du tungstène. L'écrouissage par chacun des écoulements est traduit par des variables cinématiques. L'identification numérique des paramètres a été réalisée par un optimiseur, couplé à un code de calcul par éléments finis simulant l'essai de flexion. La loi ainsi établie a été finalement validée en confrontant ses prévisions au comportement observé expérimentalement pour des trajets de chargement complexes. Cette loi de comportement a finalement été utilisée pour simuler les conditions d'utilisation d'une anode réelle. On a ainsi calculé, à l'aide d'un maillage bidimensionnel axisymétrique, les contraintes générées par le refroidissement post-recuit, une puis plusieurs séries de clichés radiographiques, et enfin un refroidissement total après utilisation. La répétition de séries de clichés conduit rapidement à des cycles stabilisés. Contrairement aux simples simulations thermoélastiques réalisées préalablement à ce travail, les niveaux de contraintes calculés sont réalistes et demeurent inférieurs à la résistance à la rupture des matériaux. Cette simulation peut d'ores et déjà être utilisée industriellement pour évaluer l'influence d'une modification de la géométrie de l'anode et/ou des conditions de sollicitation sur les contraintes en service.

L'objectif principal de cette thèse est la modélisation de matériaux supraconducteurs, et plus particulièrement, de ceux dits à haute température critique. Ce travail est divisé en deux parties distinctes : dans un premier temps, nous nous intéressons à la reconstruction de la distribution dynamique du courant dans la section d'un ruban supraconducteur à l'aide d'un modèle décrivant son comportement global. Ensuite, nous nous focalisons sur la modélisation électromagnétique et thermique par éléments finis des supraconducteurs en régime sur-critique, c'est à dire au-delà du courant critique Ic. Une introduction à la supraconductivité est faite dans le premier chapitre de ce travail, présentant les propriétés physiques principales que possède un matériau supraconducteur. Ensuite, nous décrirons les différents modèles existants permettant de décrire le comportement d'un échantillon supraconducteur. En se basant sur le modèle de Bean, nous proposons une nouvelle méthode permettant d'évaluer la distribution du courant dans la section d'un ruban monofilamentaire à partir de la mesure du profil magnétique à sa surface. Les résultats obtenus sont comparés avec ceux données par une modélisation certes plus précise, mais aussi plus complexe : l'utilisation de la méthode des éléments finis. Cette dernière sera largement utilisée dans la seconde partie de ce travail, qui traite de la modélisation de matériaux en régime sur-critique. Ceci nécessite d'introduire une formulation différente de l'habituelle loi de puissance Ec(J/Jc)^n pour décrire le comportement électrique du matériau dans une plage de courant plus importante, incluant les régimes sur-critiques. Pour cela, nous proposons une méthode originale, permettant d'approcher des mesures effectuées sur des échantillons d'YBCO et autorisant la description de la transition de l'état supraconducteur vers l'état normal. La dépendance en température a aussi été introduite afin d'étudier les phénomènes thermiques qui se produisent. Il est ainsi possible de résoudre de manière couplée les équations électromagnétiques et thermiques du problème : la résolution d'un pas de temps électromagnétique conduit à la connaissance des pertes locale. Ces pertes sont injectées dans le modèle thermique, et permet de calculer l'augmentation de température correspondante. Ce résultat est alors utilisé pour modifier les paramètres électriques en vue de la résolution du prochain pas de temps. Finalement, cette méthode a été utilisée afin de simuler le comportement d'un limiteur de courant supraconducteur. Le comportement global du système peut être reproduit par le modèle numérique implémenté, autorisant ainsi l'étude des variables locales, telles que la densité de courant ou le profil de température. Les résultats obtenus peuvent être utilisés afin d'optimiser les performances du dispositif selon des critères spécifiques. C'est ainsi que nous avons proposé une modification de la géométrie permettant d'éviter un possible emballement thermique qui pourrait conduire à la destruction du système. Sur la base de cette proposition, des modifications ont été apportées au design du dispositif, et sont actuellement testées.

<p align="justify">Les fonctions de Lagrange sont des fonctions indéfiniment dérivables qui s'annulent en tous les points d'un réseau sauf un. Ces fonctions sont utilisées comme fonctions de base d'un calcul variationnel. Les éléments de matrice de ce calcul sont évalués à l'aide de la règle de quadrature de Gauss définie par le réseau de points. Les équations à résoudre prennent ainsi la forme d'équations sur réseau.</p> <p align="justify">La méthode des réseaux de Lagrange allie simplicité et précision. La matrice représentant le potentiel est diagonale et ne dépend que des valeurs prises par le potentiel aux points du réseau. Contrairement à la méthode des différences finies, une expression analytique est obtenue pour la solution. Nous cherchons clans cette thèse à cerner les avantages et inconvénients de la méthode des réseaux de Lagrange, ainsi qu'à étendre son champ d'application en mécanique quantique. Nous montrons notamment que cette méthode peut être reliée à d'autres méthodes sur réseau, telles que les méthodes de la variable discrétisée (DVR) ou du réseau de Fourier, qui sont fort utilisées en physique atomique et moléculaire.</p> <p align="justify">Dans les problèmes à deux corps, nous appliquons la méthode à l'étude des états liés et nous l'étendons au cas des collisions, c'est-à-dire aux états libres. Une nouvelle technique de calcul de la longueur de diffusion et de la portée effective est également considérée. Dans certains cas, la solution exacte du problème à deux corps existe sous forme analytique, ce qui permet une étude de la précision de la méthode en ce qui concerne les valeurs propres et les vecteurs propres de la matrice hamiltonienne. L'extension de la méthode aux problèmes à deux corps régis par une dynamique semi-relativiste est également examinée.</p> <p align="justify">Dans le cas des problèmes à trois corps, nous effectuons une comparaison entre plusieurs systèmes de coordonnées auxquels sont couplés différents réseaux de Lagrange. Les résultats de cette comparaison dépendent de la présence de singularités dans les potentiels, celles-ci pouvant limiter fortement la précision de la méthode.</p> <p align="justify">En physique nucléaire, nous comparons deux approches sur réseaux de Lagrange lors de l'étude de l'état fondamental du noyau 6He. Il s'agit d'un noyau à halo de neutrons, pour lequel il existe une forte probabilité de trouver deux des neutrons loin des autres nucléons. Le noyau 6He peut ainsi être traité comme un système à trois corps, constitué d'une particule alpha et de deux neutrons. Nous étendons également le modèle à trois corps pour ce noyau au cas d'interactions à deux corps plus générales, c'est-à-dire contenant différents opérateurs agissant sur les spins des nucléons.</p> <p align="justify">En physique atomique et moléculaire, où les interactions sont, en première approximation, purement coulombiennes, nous nous sommes intéressé aux états S et P des principaux systèmes à trois corps que sont l'atome d'hélium He, les ions hydrogène H-et positronium Ps-, l'ion moléculaire d'hydrogène HZ et la molécule muonique dt"mu". Les fonctions d'onde approchées obtenues lors de la détermination des états liés sont utilisées pour évaluer des rayons quadratiques moyens et les rayons de masse de ces systèmes.</p>

.