Dans le domaine aéronautique, la détection des défauts sur structures composites se fait lors des phases de maintenance. Celles-ci sont prévues de manière à détecter les défauts avant que ceux-ci ne deviennent critiques pour la tenue des structures en service. L’optimisation des intervalles de maintenance représente donc un fort intérêt que ce soit au niveau financier ou opérationnel. Les dommages des structures composites sont aujourd’hui traités par la non-propagation des défauts : lorsque ceux-ci apparaissent et qu’ils sont détectés, les pièces sont immédiatement réparées ou remplacées. Afin d’améliorer ces intervalles de maintenance, le comportement des défauts soumis aux spectres de charges aéronautique nécessite d’être caractérisé. Cela permettrait une transition vers une philosophie de tolérance aux dommages, en définissant des phases de propagation non-critique pour la tenue mécanique des structures. Dans ce contexte, la propagation de délaminage est étudiée sous chargement de fatigue. En mode I de propagation, les effets de la fréquence de sollicitation et du rapport de charge sur la propagation de délaminage en fatigue sont déterminés. Ceux-ci amorcent l’étude de la propagation de délaminage sous spectre de charges à niveau variable. Les effets de l’historique de chargement sont mis en évidence. En mode II de propagation, les effets de la fréquence de sollicitation sur la propagation de délaminage sont abordés pour quatre matériaux à l’aide d’un montage d’essais vibratoires dédié. Une analyse thermique des essais est également conduite. De plus, les effets du rapport de charge sont déterminés pour ce mode de propagation.

Depuis de nombreuses années, le givrage a été identifié comme un danger dans le domaine de l’aéronautique. L’accrétion de givre se produit lorsque des gouttelettes d’eau surfondue se déposent sur une surface, en particulier le bord d’attaque d’une aile ou la lèvre d’entrée d’air moteur, et gèlent après l’impact. Ceci peut ensuite engendrer une dégradation des performances aérodynamiques, un dysfonctionnement des sondes ou encore un endommagement du moteur. C’est pourquoi cette problématique est étudiée avec attention. Les essais en vol et en soufflerie étant longs et coûteux, la simulation numérique de l’accrétion de givre est devenue un outil nécessaire dans le processus de conception et de certification des avions. Cette thèse s’inscrit dans le contexte de la modélisation 3D de l’accrétion de givre, et plus particulièrement des couches limites dynamique et thermique qui se développent autour du corps givré. Les outils numériques devant être rapides et robustes, l’approche proposée dans cette thèse pour le calcul aérodynamique est une méthode couplée Euler/couche limite intégrale. Ainsi, un modèle intégral est développé pour représenter le développement de la couche limite dynamique. La partie thermique est modélisée soit par une méthode simplifiée basée sur des approches algébriques, soit par une méthode intégrale. Cette modélisation des couches limites dynamique et thermique est valable sur paroi lisse ou rugueuse et permet de fournir notamment le coefficient de frottement et le coefficient d’échange thermique qui sont nécessaires pour un calcul d’accrétion de givre. Les équations intégrales de couche limite, associées à leurs relations de fermeture, sont ensuite résolues par une méthode Volumes-Finis sur maillage surfacique non structuré, qui est bien adaptée pour les géométries complexes. De plus, des traitements numériques spécifiques sont mis en œuvre pour améliorer la précision de la méthode au voisinage du point d’arrêt et pour rendre le code robuste au passage du décollement. Après la validation de la méthode de couche limite, le code est utilisé dans les chaînes de givrage 2D et 3D de l’ONERA pour des applications d’accrétion de givre. Ceci permet de montrer l’intérêt de la méthode en termes de robustesse et de précision par rapport aux codes de couche limite habituellement utilisés dans les codes de givrage actuels.

La sûreté de fonctionnement occupe une place prépondérante dans la conception de systèmes critiques, puisqu'un dysfonctionnement peut être dangereux pour les utilisateurs ou l'environnement. Les concepteurs doivent également démontrer aux autorités de certification que les risques encourus sont acceptables. Pour cela, le concepteurs définissent une architecture contenant un ensemble de mécanismes de sûreté permettant de mitiger ou tout du moins limiter la probabilité d’occurrence des risques identifiés. L'objectif de ce travail est de développer une méthode automatique et générique de synthèse d’architecture assurant formellement le respect d’exigences de sûreté. Cette activité de synthèse est formalisée comme un problème d'exploration de l'espace des architectures c'est-à-dire trouver un candidat appartenant à un espace de recherche fini, respectant les exigences de sûreté. Ainsi nous proposons un processus de résolution complet et correct des problèmes d'exploration basé sur l'utilisation des solveurs SMT. Les contributions principales sont: 1- La formalisation de la synthèse comme un problème de Satisfiabilité Modulo Théorie (SMT) afin d’utiliser les solveurs existants pour générer automatiquement une solution assurant formellement le respect des exigences; 2- Le développement de méthodes d’analyse spécialement conçues pour évaluer efficacement la conformité d’une architecture vis-à-vis d’un ensemble d’exigences; 3- La définition d'un langage KCR permettant de formuler les problèmes d'exploration et l'implantation des méthodes de résolution au sein de l'outil KCR Analyser.

Dans cette thèse, le problème du contrôle d’attitude de satellites agiles à l’aide de grappes d’actionneurs gyroscopiques (AGs) est considéré et plus particulièrement son application au contrôle de micro/nanosatellites (10-100 kg). Des outils d’analyse topologique sont tout d'abord développés. La comparaison de différentes configurations de grappes justifie le choix d'une géométrie pyramidale à six actionneurs gyroscopiques. Des analyses plus approfondies de cette grappe, avec et sans cas de panne, permettent de déduire des contraintes que la loi de pilotage doit vérifier pour être adaptée à ce système. Pour y répondre, après analyse de la littérature, une nouvelle structure de loi de pilotage ainsi qu’une formulation différente des équations cinématiques sont développées. Cette structure est basée sur l’algorithme du filtre de Kalman étendu. Elle a pour avantages de répondre aux exigences en termes de calcul temps réel au bord des satellites, de flexibilité sur la gestion des contraintes et de facilité d’adaptation en cas de pannes. En outre, une procédure de génération de boucle de commande, englobant la loi de pilotage et un contrôleur robuste du système, est proposée. La généralisation de cette boucle de commande est illustrée sur des bras manipulateurs à base fixe et spatiaux. En parallèle, l’étude du passage des singularités internes intraversables mène à une nouvelle stratégie d’évitement de ces singularités. Des simulations sur des modèles de satellites représentatifs illustrent les résultats. La grappe d’actionneurs et la boucle de commande développées seront testées dans le cadre d’une expérimentation en microgravité.

L'objectif de cette thèse est de contribuer au développement et à la validation d'outils numériques permettant la Simulation aux Grandes Echelles (SGE) de l'allumage d'un écoulement turbulent diphasique dans une chambre de combustion. Pour ce faire, une méthode de dépôt d'énergie modélisant l'apport d'énergie lié au claquage de la bougie d'allumage a été implémentée dans la chaîne de calcul CEDRE. Cette méthode a été validée sur une simulation de l'allumage d'un écoulement laminaire purement gazeux d'air et de propane. Une SGE de l'écoulement d'air du MERCATO a été effectué à l'aide du solveur Navier-Stokes CHARME de CEDRE. Cette simulation reproduit fidèlement l'écoulement turbulent non-réactif dans la chambre de combustion. Une méthode d'injection simplifiée FIMUR a été ajoutée au solveur lagrangien SPARTE de CEDRE. Dans cette méthode, des gouttes sont injectées directement au nez de l'injecteur avec une distribution de vitesse et de taille imposée. Une SGE de l'écoulement turbulent diphasique dispersé non-réactif dans la chambre MERCATO a ensuite été réalisée avec cette méthode. La comparaison des champs particulaires moyens de vitesse et de taille obtenus par simulation numérique avec les données expérimentales est satisfaisante. Enfin, des SGE de l'allumage de la chambre MERCATO ont été effectuées à partir du champ diphasique non-réactif simulé et de la méthode de dépôt d'énergie développée. Selon l'instant du dépôt d'énergie, les simulations conduisent à des allumages réussis ou ratés. La propagation de la flamme dans la chambre pour un allumage réussi a fait l'objet d'une analyse détaillée pour tenter de déterminer les principaux facteurs l'influençant.

Dans les systèmes informatiques temps réel, les tâches logicielles sont contraintes par le temps. Pour garantir la sûreté du système critique controlé par le système temps réel, il est primordial d'estimer de manière sûre le pire temps d'exécution de chaque tâche. Les performances des processeurs actuels du commerce permettent de réduire en moyenne le temps d'exécution des tâches, mais la complexité des composants d'optimisation de la plateforme rendent difficile l'estimation du pire temps d'exécution. Il existe différentes approches d'estimation du pire temps d'exécution, souvent ségréguées et difficilement généralisables ou au prix de modèles coûteux. Les approches probabilistes basées mesures existantes sont vues comme étant rapides et simples à mettre en œuvre, mais souffrent d'un manque de systématisme et de confiance dans les estimations qu'elles fournissent. Les travaux de cette thèse étudient les conditions d'application de la théorie des valeurs extrêmes à une suite de mesures de temps d'exécution pour l'estimation du pire temps d'exécution probabiliste, et ont été implémentées dans l'outil diagxtrm. Les capacités et les limites de l'outil ont été étudiées grâce à diverses suites de mesures issues de systèmes temps réel différents. Enfin, des méthodes sont proposées pour déterminer les conditions de mesure propices à l'application de la théorie des valeurs extrêmes et donner davantage de confiance dans les estimations.

L’estimation de la charge d’un satellite et du risque de décharge nécessite dans certains cas la prise en compte dans les modèles numériques d’échelles spatiales très différentes. En particulier, les interconnecteurs présents à la surface des générateurs solaires d’un satellite sont susceptibles de modifier son équilibre électrostatique lors de missions spatiales rencontrant un environnement plasma dense. Une modélisation classique de cet effet nécessiterait le maillage d’éléments à des échelles submillimétriques, sur un satellite de plusieurs dizaines de mètres d’envergure, ce qui rendrait la simulation extrêmement onéreuse en temps de calcul. De plus, ces interconnecteurs sont parfois fortement chargés positivement par rapport à l’environnement, ce qui empêche l’application du modèle de Maxwell-Boltzmann classiquement utilisé pour les populations d’électrons. Dans une première partie, nous avons développé une méthode itérative de type Patch adaptée à la résolution du problème non-linéaire de Poisson-Boltzmann pour la simulation du plasma spatial. Cette méthode numérique multigrille permet la simulation de l’impact d’éléments de petite taille à la surface d’un satellite complet. Dans une seconde partie, nous avons développé un schéma correctif permettant d’utiliser le modèle de Maxwell-Boltzmann pour la population d’électrons, malgré la présence de surfaces satellites chargées positivement, en y ajoutant un terme de correction calculé à l’aide de la méthode Particle-in-Cell. Nous avons montré que ce schéma permet, tout en limitant le coût en calculs, de déterminer avec précision les courants collectés par les surfaces du satellites, qu’elles soient chargées négativement ou positivement.

La réduction des émissions polluantes et l'amélioration des performances des turboréacteurs nécessitent une connaissance détaillée des phénomènes physiques mis en jeu dans une chambre de combustion. L'atomisation du carburant résulte du cisaillement engendré par un fort écoulement d'air généré dans l'injecteur. La simulation numérique directe d'écoulements avec interface permet de simuler l'ensemble du processus d'atomisation. L'utilisation de maillages Cartésiens permet la réalisation de calculs HPC efficaces et précis. Mais, une des complexités de l'atomisation vient d'une interaction forte entre le comportement de la nappe liquide et l'écoulement gazeux dans les conduites de l'injecteur, rendant impératif la simulation de l'injecteur complet. Ceci étant impossible avec des maillages Cartésiens structurés, l'objectif de cette thèse est de développer une méthode d'interface immergée permettant l'inclusion d'objets solides dans un domaine de calcul, indépendamment du maillage, afin de réaliser des DNS du système d'injection complet. Les équations de Navier-Stokes incompressibles diphasiques sont résolues à l'aide d'un algorithme de projection, l'interface liquide-gaz étant transportée avec une méthode CLSVOF conservative en masse et quantité de mouvement. La présence du solide est prise en compte grâce à la méthode d'interface immergée. Cette méthode a été appliquée à la simulation numérique de nappes liquides cisaillées pour une configuration d'injecteur utilisée en essais à l'ONERA et a permis une meilleure prédiction de la fréquence de battement de la nappe.

Dans cette thèse, nous proposons de construire de meilleurs modèles Processus Gaussiens (GPs) en intégrant les connaissances antérieures avec des données expérimentales. En raison du coût élevé de l’exécution d’expériences sur les systèmes physiques, les modèles numériques deviennent un moyen évident de concevoir des systèmes physiques. Traditionnellement, ces modèles ont été construits expérimentalement et itérativement ; une méthode plus rentable de construction de modèles consiste à utiliser des algorithmes d’apprentissage automatique. Nous démontrons comment créer des modèles en intégrant une connaissance antérieure en modifiant les fonctions de covariance. Nous proposons des modèles GP pour différents phénomènes physiques en mécanique des fluides. De même, les lois physiques entre plusieurs sorties peuvent être appliquées en manipulant les fonctions de covariance. Pour chaque application, nous comparons le modèle proposé avec le modèle de l’état de l’art et démontrons les gains de coût ou de performance obtenus.

Dans le cadre des systèmes de positionnement par satellite GNSS (« Global Navigation Satellite Systems »), l’intégrité de la navigation d’un utilisateur est gérée en réception par la détection, l’identification voire l’exclusion de mesures de pseudo-distance jugées erronées. Généralement basés sur le concept a posteriori RAIM (« Receiver Autonomous Integrity Monitoring »), les algorithmes de contrôle autonome d’intégrité fournissent de hautes performances pour l’aviation civile, dont le contexte de navigation est caractérisé par une forte visibilité des satellites et peu de signaux parasites captés par l’antenne réceptrice. L’algorithme WLSR RAIM est communément utilisé dans ce cadre. Néanmoins, les techniques RAIM ne sont pas compatibles avec la navigation terrestre en milieu contraint. En effet, le contexte urbain est notamment caractérisé par un masquage récurrent des signaux satellitaires directs ainsi que la réception de multi-trajets générés par l’environnement proche du récepteur. RAIM ne prend pas en compte l’ensemble des données disponibles en réception, dégradant ainsi fortement ses performances. Il est donc nécessaire de développer des méthodes de contrôle d’intégrité compatibles avec un tel contexte de navigation. Pour cela, la thèse propose d’étudier l’apport d’informations GNSS a priori non utilisées par les techniques RAIM. Deux paramètres principaux ont été exploités : le signal GNSS brut reçu et les estimations de directions d’arrivée des signaux satellitaires DOA (« Direction Of Arrival »). La première étape a consisté à implémenter une méthode a priori qui évalue la cohérence du positionnement estimé par rapport au signal brut directement reçu. Cette méthode a été nommée Direct-RAIM (D-RAIM) et a démontré une forte sensibilité de détection, permettant d’anticiper d’éventuels risques sur la navigation et de caractériser plus finement la qualité de l’environnement proche du récepteur. Toutefois, le caractère a priori de l’approche engendre de potentielles non détection d’erreurs en cas de modèle de signal défectueux. Afin de contourner cette limitation, un couplage WLSR RAIM – D-RAIM a été développé, nommé Hybrid-RAIM (H-RAIM). Une telle approche permet de combiner robustesse et sensibilité apportées par ces techniques respectives. Le second axe de recherche a mis en évidence la contribution de l’information des DOA dans un contrôle autonome d’intégrité. L’intégration d’un réseau d’antennes en réception permet d’obtenir l’estimation des DOA pour l’ensemble de la constellation visible. Théoriquement, l’évolution jointe des DOA est directement liée à l’attitude du réseau. Cet aspect permet donc de détecter toute incohérence sur une ou plusieurs voies en cas d’estimation(s) de DOA biaisée(s), par rapport à l’ensemble de la constellation. L’algorithme RANSAC (« RANdom SAmple Consensus») a été utilisé afin de détecter tout comportement aberrant dans l’estimation des DOA, et ainsi mesurer la confiance que l’utilisateur peut placer dans chaque voie. L’algorithme WLSR RAIM RANSAC a ainsi été implémenté. L’intégration de la composante DOA permet d’ajouter un degré de liberté dans le contrôle autonome d’intégrité côté récepteur et ainsi d’affiner la détection voire l’exclusion d’erreurs. Au cours de cette thèse, un récepteur logiciel a été implémenté, permettant de traiter des signaux Galileo, de la génération du signal jusqu’au positionnement puis au contrôle d’intégrité. Ce récepteur a pu être évalué à partir de données simulées en environnement urbain.