Quelques développements récents en traitement du signal

Ce travail porte sur le développement d'une méthode Galerkin discontinue (GDDT) d'ordre élevé pour la résolution numérique des équations de Maxwell instationnaires sur des maillages simplexes non-conformes. On présente tout d'abord une méthode GDDT reposant sur des fonctions de base nodales pour approcher le champ électromagnétique dans un simplexe, un schéma centré pour évaluer les flux numériques aux interfaces entre cellules voisines et un schéma saute-mouton du second ordre pour l'intégration temporelle. De plus, cette méthode autorise l'utilisation de maillages non-conformes présentant un nombre arbitraire de noeuds flottants. La méthode résultante est non-dissipative, stable sous une condition de type CFL, conserve un équivalent discret de l'énergie électromagnétique, et très peu dispersive. Afin de diminuer le coût de calcul de cette méthode, on propose une méthode GDDT de type /hp/, qui combine /h-/raffinement et /p/-enrichissement locaux tout en préservant la stabilité. On réalise ensuite une étude numérique détaillée des méthodes GDDT sur la base d'une série de problèmes de propagation d'ondes en milieux homogène et hétérogène. En particulier, on effectue une comparaison des méthodes Galerkin discontinues conformes et non-conformes en termes de précision, convergence et coûts de calcul.<br />Afin d'améliorer la précision et la vitesse de convergence des méthodes GDDT précédentes, on étudie une famille de schémas saute-mouton d'ordre<br />arbitrairement élevé. Ces schémas temporels nous assurent sur tout maillage la conservation d'un équivalent discret de l'énergie électromagnétique ainsi que la stabilité des méthodes GDDT résultantes sous une condition de type CFL. On réalise aussi une étude de convergence /hp a priori/ ainsi qu'une étude de convergence de l'erreur sur la divergence. Des expériences numériques montrent que pour un maillage donné, le schéma saute-mouton du quatrième ordre est moins coûteux en temps de calcul et plus précis que le schéma saute-mouton du second ordre, en dépit d'une complexité arithmétique accrue.<br />De plus, on obtient une convergence exponentielle avec le schéma saute-mouton du quatrième ordre.

These thesis review some recent results on the construction of very high order multidimensional upwind schemes for the solution of steady and unsteady conservation laws on unstructured triangular grids. We also consider the extension to the approximation of solutions to conservation laws containing second order dissipative terms. To build this high order schemes we use a sub-triangulation of the triangular Pk elements where we apply the distribution used for a P1 element. This manuscript is divided in two parts. The first part is dedicated to the design of the high order schemes for scalar equations and focus more on the theoretical design of the schemes. The second part deals with the extension to system of equations, in particular we will compare the performances of 2nd, 3rd and 4th order schemes. The first part is subdivided in four chapters: The aim of the second chapter is to present the multidimensional upwind residual distributive schmes and to explain what was the status of their development at the beginning of this work. The third chapter is dedicated to the first contribution: the design of 3rd and 4th order quasi non-oscillatory schemes. The fourth chapter is composed of two parts: We start by understanding the non-uniformity of the accuracy of the 2nd order schemes for advection-diffusion problem. To solve this issue we use a Finite Element hybridisation. This deep study of the 2nd order scheme is used as a basis to design a 3rd order scheme for advection-diffusion. Finally, in the fifth chapter we extend the high order quasi non-oscillatory schemes to unsteady problems. In the second part, we extend the schemes of the first part to systems of equations as follows: The sixth chapter deals with the extension to steady systems of hyperbolic equations. In particular, we discuss how to solve some issues such as boundary conditions and the discretisation of curved geometries. Then, we look at the performance of 2nd and 3rd order schemes on viscous flow. Finally, we test the space-time schemes on several test cases. In particular, we will test the monotonicity of the space-time non-oscillatory schemes and we apply residual distributive schemes to acoustic problems.

Le principe du contrôle par flux de fuite magnétique (Magnetic Flux Leakage MFL) consiste à aimanter la pièce à contrôler par un champ magnétique et à détecter à l'aide d'un capteur magnétique les fuites des lignes du champ qui résultent de la présence d'un défaut dans la pièce. Dans le but d'améliorer les performances d'un dispositif de détection, le CEA et la société Vallourec collaborent pour développer un modèle numérique dédié au contrôle virtuel des défauts longitudinaux dans les tubes ferromagnétiques. Le dispositif expérimental comprend un circuit magnétique tournant à une vitesse constante autour du tube qui défile. Dans le cadre de cette thèse, on débute le problème de la modélisation sans tenir compte des effets de la vitesse de rotation, il s'agit donc de résoudre un problème d'électromagnétisme en régime magnétostatique.Pour résoudre ce problème, on propose de comparer une approche semi-analytique basée sur le formalisme des équations intégrales (EI) et une approche purement numérique utilisant les éléments finis (EF).Dans la première partie de cette thèse, après avoir établi le formalisme théorique par EI, un premier modèle considérant des matériaux ferromagnétiques à perméabilité magnétique constante (régime linéaire) a été mis en œuvre en 2D. Ce modèle a été appliqué pour un exemple de système extrait de la littérature et validé numériquement par une comparaison des résultats EI/EF. Pour une meilleure détection, il est opportun de saturer magnétiquement la pièce. Le matériau ferromagnétique est alors caractérisé par une courbe B(H) non-linéaire. Par conséquent, la deuxième partie de la thèse a été consacrée à la mise en œuvre du modèle en régime non linéaire qui tient compte de cette caractéristique.Différentes méthodes de discrétisation ont été étudiées afin de réduire le nombre d'inconnues et le temps de calcul. L'originalité de la thèse réside dans l'utilisation des fonctions d'interpolation d'ordre élevé (polynôme de Legendre) pour une discrétisation des équations intégrales par une approche de type Galerkin. Les premiers essais de validation numérique de ce modèle ont été effectués sur un système MFL simplifié. Des premiers essais de validation expérimentale pour des données obtenues par EF ont été effectués en deux phases : La première a consisté à vérifier le distribution du champ magnétique pour un tube sain et en régime magnétostatique. La deuxième phase a consisté à calculer la réponse d'un défaut dans le tube ferromagnétique en tenant en compte les effets éventuels de la rotation du circuit magnétique par rapport au tube.

Au cours de ce travail, une implantation de la méthode d'analyse de sensibilité d'ordre élevé a été réalisée au sein de la méthode des éléments finis (MEF). Tout d'abord nous avons posé ta base théorique des développements: la MEF classique et surtout la méthode de dérivées d'ordre élevé qui engendre les résultats sous la forme de développements de Taylor par rapport aux paramètres de conception. En s'appuyant sur la structure typique d'un code de calculs EF, nous avons développé ensuite les modules d'un calcul "paramétré": la paramétrisation géométrique, la dérivation du maillage et la résolution paramétrée. L'originalité de cette implantation consiste dans une nouvelle organisation des calculs de dérivées qui sont basés notamment sur les opérations symboliques. Pour valider des étapes de calculs, plusieurs exemples ont été présentés. Enfin, les applications possibles des résultats paramétrés sont évoquées, entre autres, une procédure pour définir le modèle analytique équivalent de dispositifs électromagnétiques a été proposée. Notons que plusieurs problèmes de réalisation informatique ont été mise en évidence et résolus, ce qui représente un gros investissement en programmation. Des nombreuses perspectives sont ouvertes par cette approche, le travail mériterait donc d'être poursuivi.

Cette étude porte sur la modélisation et la simulation numérique d'un piano, en domaine temporel, par modèles phy- siques. Nous souhaitons rendre compte du comportement vibratoire et acoustique du piano, en prenant en compte les éléments principaux qui contribuent à la production du son. La table d'harmonie est modélisée par une équation bidimensionnelle de plaque épaisse, le système de Reissner Mindlin, pour un matériau orthotrope et hétérogène, dont l'amortissement dépend de la fréquence. Grâce aux équations de la vibroacoustique, la table rayonne dans l'air, dans lequel on souhaite calculer le champ acoustique complet autour de la ceinture du piano, que l'on suppose rigide. La table d'harmonie est d'autre part sollicitée par les cordes, à travers le chevalet où elles présentent un léger angle par rapport au plan horizontal. Chaque corde est modélisée par un système d'équations monodimensionnelles amorties dans lequel on prend en compte non seulement les ondes transversales excitées par le marteau, mais aussi la raideur à travers les ondes de cisaillement, ainsi que le couplage avec les ondes longi- tudinales provenant de la prise en compte des non linéarités géométriques. Le marteau est lancé avec une vitesse initiale vers un chœur de cordes, contre lequel il s'écrase avant d'être repoussé par les cordes. La force d'interaction dépend de façon non linéaire de l'écrasement du marteau.Le modèle complet de piano, que l'on souhaite résoudre numériquement, consiste donc en un système couplé d'équations aux dérivées partielles, dont chacune revêt des difficultés de nature différente : la corde est régie par un système d'équations non linéaires, la table d'harmonie est soumise à un amortissement dépendant de la fréquence, la propagation acoustique requiert un très grand nombre d'inconnues; auxquelles s'ajoute la difficulté inhérente aux couplages. D'une part, la stabilité numérique du schéma discret peut être compromise par la présence d'équations non linéaires et de nombreux couplages. Une méthode efficace pour garantir cette stabilité a priori est de construire un schéma qui conserve, ou dissipe, un équivalent discret de l'énergie physique d'un pas de temps au suivant. Une contribution majeure de ce travail a été de développer des schémas préservant une énergie discrète pour une classe de systèmes non linéaires dans laquelle s'inscrit le modèle de corde. D'autre part, afin d'augmenter l'efficacité de la méthode et de réduire le coût des calculs numériques, il est souhaitable de mettre à jour de façon découplée les inconnues liées aux différentes parties du problème, sur lesquelles la discrétisation en temps est faite de façon différente, afin de s'adapter aux spécificités de chacune. L'introduction de multiplicateurs de Lagrange nous permet de réaliser ce découplage artificiel grâce à des compléments de Schur adaptés. L'utilisation du code de calcul en situation réaliste montre le potentiel d'une telle modélisation d'un piano complet en domaine temporel. Au delà de très bien reproduire les mesures, il est possible d'étudier l'influence de certains phénomènes physiques (corde raide, non linéaire), de la géométrie ou encore des matériaux utilisés sur le comportement vibratoire général du piano, et sur le son en particulier. L'enrichissement spectral, ainsi que l'apparition des " partiels fantômes " et du précurseur non linéaire sont clairement mis en évidence pour les grandes amplitudes de jeu, soulignant l'intérêt de notre approche dans la compréhension du fonctionnement de l'instrument.

La pile à combustible (PAC) est une technologie maintenant maîtrisée qui permet de convertir efficacement le combustible hydrogène en énergie électrique et thermique avec un faible impact environnemental. L'extension de son utilisation dépend fortement de la qualité de l'interface électronique de puissance qui a pour l'objet d'adapter la pile à la charge. Le travail de recherche présenté dans ce mémoire de Doctorat s'intéresse à la conception d'une interface électronique de puissance pour la pile à combustible pour les trois domaines d'applications : télécommunication, transport et stationnaire. <br />Dans cette objective, l'approche modulaire, qui se base sur l'étude des convertisseurs unitaires et leurs modes de connexion, a été proposée afin de s'adapter à la modularité des stacks de PAC. Des modèles de pertes et de dimensionnement des composants passifs et semi conducteurs ont été construits. La technique de l'entrelacement est introduite afin de résoudre le problème de fort courant et permet d'optimiser le dimensionnement des inductances. Une nouvelle structure nommée double BOOST dual entrelacé (Interleaved Double Dual BOOST – IDD BOOST) a été proposée afin de résoudre les difficultés d'un convertisseur modulaire présentant une tension de sortie élevée et un rapport de tension important. Les stratégies de contrôle – commande multi sources ont été étudié en adaptant aux topologies de l'interface de puissance. Deux prototypes ont été réalisés afin de valider les résultats de prédiction de pertes ainsi que le contrôle commande multi sources.

La génération d'harmoniques d'ordre élevé (HHG), qui dans le domaine temporel se traduit par l'émission d'un train d'impulsion VUV attoseconde (1as =10-18s), a connu un grand intérêt scientifique depuis une dizaine d'années. Cette source constitue en effet un bon candidat pour la mise en oeuvre d'expériences pompe sonde visant à observer la dynamique électronique au coeur même des atomes et des molécules. <br />Au CELIA, nous avons implémenté une technique de post-compression qui nous a permi de comprimer nos impulsions laser IR de 40 fs à 9 fs (1fs=10-15s). Ces impulsions sont ensuite utilisée pour confiner la HHG. Étant donné que le processus de HHG est efficace uniquement si les impulsions IR génératrices sont polarisées linéairement, nous avons créé une porte dans le profil temporel de nos impulsions sub-10fs où la polarisation est linéaire pendant une durée inferieure à la durée de l'impulsion IR génératrice. Ceci nous permet de confiner la HHG en dessous d'un demi-cycle optique IR. Cette technique de porte d'ellipticité, complètement caractérisée dans cette thèse, nous a permis de confiner la HHG jusqu'à l'émission d'une à deux impulsions attosecondes. Afin de caractériser le profil temporel du train d'impulsions attosecondes, nous avons également implémenté un interféromètre à deux couleurs qui nous a permit de mesurer la phase harmonique et de reconstruire nos trains d'impulsions attosecondes.

Attendre un enfant est généralement un événement heureux pour un couple. Toutefois, lorsqu’une hospitalisation de la femme survient en raison d’une grossesse à risque élevé, l’expérience des futurs parents peut se transformer en véritable cauchemar. Jusqu’à maintenant, les savoirs disponibles suggèrent que cet événement imprévu ait de nombreux impacts sur les futurs parents. Néanmoins, ce sujet n’a été abordé que selon une perspective maternelle ou paternelle, sans s’intéresser au caractère systémique de la situation. Cette étude qualitative de cas multiples propose d’explorer l’expérience de couples vivant l’hospitalisation de la femme en raison d’une grossesse à risque élevé, dans une perspective systémique et constructiviste. Des entrevues semi-dirigées ont été réalisées auprès de 7 couples dont la femme était, au moment de l’entrevue, hospitalisée en raison d’une grossesse à risque élevé. L’analyse qualitative des données s’est inspirée de la thématisation et a tenu compte des cadres théoriques et épistémologique choisis, étant respectivement l’approche systémique familiale de Wright et Leahey (2013) et le constructivisme. Les résultats suggèrent que l’hospitalisation est marquée par l’intensité des émotions ressenties. Les futurs parents ressentent presque constamment des sentiments d’inquiétude, d’incertitude et de solitude. Il est toutefois intéressant de constater que le couple lui-même subit des changements au cours de l’hospitalisation. En ce sens, cet événement requiert une immense réorganisation, laquelle touche plusieurs systèmes gravitant autour du système conjugal. Enfin, afin de faire face à l’intensité et à l’imprévisibilité de la situation, les couples sont amenés à solliciter un soutien externe et à puiser dans leurs propres ressources internes. En lien avec les cadres théoriques et épistémologiques choisis, ces résultats conviennent de la nécessité d’adopter une vision systémique et contextuelle afin de développer des soins infirmiers concordants avec l’expérience des couples. Cela ramène à la nécessité de réaliser d’autres études afin de perfectionner la compréhension de l’expérience. Certaines pistes d’interventions infirmières sont néanmoins énoncées, afin d’amener les infirmières à considérer la situation du point de vue de l’expérience conjugale, dont les défis de réorganisations extérieurs et des relations conjugales. Enfin, il est proposé que l’infirmière puisse agir afin de renforcer la capacité des couples à composer avec les conséquences de l’hospitalisation anténatale. / Planning to have a child is most of the time a joyful event for a couple. However, the same event might becomes a real nightmare when the woman is hospitalized because of an at high risk pregnancy. To date, researches exploring this abrupt event have suggested that it may conduct to numerous impacts on future parents. Even though, the subject has only been explored using a maternal or paternal perspective, without considering the systemic aspect of the situation. This qualitative multiple case study suggests to explore the couples’ experience regarding the hospitalization of the woman because of an at high risk pregnancy, in a systemic and constructivist approach. Semi-directed interviews have been conducted with seven couples with a woman hospitalized because of an at high risk pregnancy, at the moment of the interview. Qualitative data analysis was inspired of thematization and took into account the theoretical and epistemological approach, respectively the Wright and Leahey (2013) family approach and constructivism. Results suggest the hospitalization is touched by intense emotions, such as feeling almost all the time worries, uncertainties and loneliness. It is however interesting to note that even the couple himself is touched by changes resulting of the hospitalization. According to that, this event requires a huge reorganization, which implies loads of systems gravitating around the couple’s system. Lastly, in order to face the intensity and the unpredictability of the situation, couples are asking for external support, while using their own internal resources. Concerning the theoretical and epistemological approach chosen to realize this study, results confirm the necessity of adopting a systemic and contextual vision in order to develop nursing cares appropriate with couples’ experience. Besides, it becomes essential to duplicate this type of study to perfection the actual comprehension of the event. However, nursing interventions relating to the external reorganization and the couple relation reorganization are proposed, in order to reinforce nurses’ consideration of the couple’s experience. Only then, it is suggested that nurses can influence couples’ capacity to face the consequences of the antenatal hospitalization.

Cette thèse porte sur la modélisation des ondes élastodynamiques dans deux situations particulières qui pénalisent les méthodes numériques utilisées pour simuler ces phénomènes. Dans la première partie, on se place dans le cas où les ondes de pression (ondes P) se propagent à une vitesse beaucoup plus grande que celle des ondes de cisaillement (ondes S). Les modèles numériques utilisés habituellement pour traiter cette configuration sont pénalisés par la plus petite vitesse qui dicte le choix du pas du schéma. Nous proposons ici un schéma qui découple numériquement, dans le volume, les ondes P et les ondes S, pour deux types de conditions de bord en utilisant la décomposition du déplacement en potentiels de Lamé, en deux dimensions. Les conditions aux limites de Dirichlet homogènes, qui sont des conditions essentielles pour la formulation classique en déplacement, deviennent des conditions naturelles, mais non standard, pour la formulation en potentiels qui se présente comme un système de deux équations d'ondes couplées par les conditions aux limites. Cette formulation préserve une énergie équivalente à l'énergie élastodynamique. Nous construisons un schéma éléments finis en espace et utilisons un thêta-schéma en temps sur les termes de bord afin de ne pas pénaliser la CFL et mener à une condition sur le pas de temps indépendante des termes de couplage au bord. Ce schéma préserve une énergie discrète. Le cas des conditions de surface libre mène à des instabilités. Nous les avons traitées comme des perturbations des conditions de Dirichlet, ce qui permet d'obtenir de bons résultats dans le domaine fréquentiel mais donne naissance à de sévères instabilités après discrétisation en temps. La seconde partie de la thèse est consacrée à la construction, l'analyse et la validation de conditions de transmission effectives (CTE) à travers une couche mince de matériau homogène et isotrope d'épaisseur constante h. Ici, la finesse de la couche affecte les schémas explicites usuels car le maillage de la couche avec des éléments suffisamment petits entraîne une diminution analogue du pas de temps critique via la condition CFL, tandis que l'on espère avec les CTE obtenir un pas de temps indépendant de l'épaisseur de la couche. Une analyse complète du cas de la bande mince rectiligne est donnée en deux et trois dimensions. Les conditions obtenues sont stables via la conservation d'une énergie et l'ordre de l'erreur d'approximation par rapport à l'épaisseur de la couche pour les conditions d'ordre 2 est de O(h^3). Des résultats numériques sont présentés pour les configurations bi et tridimensionnelles, ils valident les résultats de stabilité, d'estimation d'erreur et de conditions de stabilité de schémas en temps proposés, qui sont des modifications du schéma explicite utilisé en l'absence de couche mince. Enfin, le traitement d'une couche curviligne est effectué dans le cas bidimensionnel. Sa stabilité est à nouveau vérifiée par conservation d'énergie et des résultats numériques sont également présentés.