.

Chaque département français s’inscrit dans une histoire et une géographie, plus particulièrement ceux qui constituent ce qu'il est dorénavant appelé outre-mer. Au moment où les départements métropolitains viennent de connaître des élections selon des nouvelles modalités conduisant aux conseils départementaux sans qu'il ne soit plus question de leur suppression, Mayotte vient de s'affirmer, par la volonté référendaire du peuple comme le 101ème département français. Après une longue histoire née de la colonisation, les Comores n'ont pas choisi le même destin puisque ces îles de l'Océan Indien sont devenues La République des Comores faisant d'ailleurs l'objet d'un coup d'Etat et un nouveau territoire de la République française passant par les différents stades institutionnels des départementaux et territoires d'Outre-mer. C'est l'histoire institutionnelle et politique de Mayotte, les étapes qui ont conduits à sa départementalisation, les obstacles rencontrés autant que les conséquences à gérer qui suscitent l'attention, l’approche et le constat du chercheur dont l'éthique de conviction (il a vécu l’histoire se faisant comme mahorais) et de responsabilité (il a travaillé avec la même objectivité de tout chercheur métropolitain) a conduit à ce travail. Il se veut contribution à l'histoire constitutionnelle de la France, tentative d'explication des évolutions législatives de l'Outre-mer, exposition des problèmes (santé, éducation, migrants, droit du sol?) à régler. Il manquait au droit constitutionnel et politique d'Outre-mer une recherche permettant de mieux saisir la raison d'un nouveau département français, le 101ème , un des seuls dont Napoléon ne soit pas l'inventeur et qui néanmoins avec force, s'inscrit dans la grande tradition de 1789 et de la « République décentralisée des années 1982-2013 ». Si comme l’écrivait Gurvitch « la démocratie est le règne du droit » gageons que Mayotte, notamment grâce au droit public et malgré la gravité de la situation économique et sociale, politique et administrative ou même juridique qui a entraîné les réactions du défenseur des droits en 2013, s’inscrit dans la démocratie au même titre que les 100 autres départements français. / Every French department has its own history and geography, not least those which form the French overseas territories and departments. While the metropolitan departments have just voted for new departmental councils, Mayotte confirmed by referendum her wish to become the 101st department. After her long colonial history, the islands of the Comoros did not choose the same destiny. These Indian Ocean islands became a republic, and Mayotte a new territory of the French Republic. It is Mayotte’s institutional and political history, the various steps which led to its departmentalisation, the obstacles overcome as much as the consequences now needing to be managed, which have attracted the researcher’s attention and analysis. The work has been driven by an ethic of conviction (the researcher has experienced this history as a Mahorais) and of responsibility (he has applied the same objectivity as every metropolitan researcher). He wants to contribute to the constitutional history of France by trying to explain in this research the legislative evolution of the French overseas departments and by exposing the different problems (health, education, migration, right to nationality based on birthplace) which need to be addressed. Up until now there hasn’t been any research done which examines the reason for this new French department; the 101st, one of the only ones not founded by Napoleon, but which nevertheless joins the great tradition of 1789 and “the decentralised republic of 1982-2013

La demande pour les ressources énergétiques est en hausse alors que leur disponibilité est en baisse. Dans ce contexte, l’industrie du transport et de l’énergie est à la recherche de petits moteurs efficaces et puissants et le Rim Rotor Rotary Ramjet Engine (R4E) pourrait correspondre à ces critères. Or, en ce moment, le potentiel de ce moteur est limité, car son compresseur supersonique entraîne des pertes d’efficacité lorsque le rotor tourne à son nombre de Mach tangentiel optimal qui est de 2. Le présent mémoire compile toutes les notions requises pour comprendre le fonctionnement d’un compresseur supersonique lors de son démarrage et de concevoir le compresseur le plus approprié pour le R4E, tant en démarrage qu’en régime permanent. Pour se faire, des concepts de cascades inspirés des compresseurs et des méthodes de démarrage des moteurs ramjet actuels ont été générés et validés à l’aide de modèles analytiques. Les concepts sont par la suite essayés expérimentalement sous la forme de cascades à l’aide d’une soufflerie supersonique. Bien que le modèle analytique montre que les cascades munies de canaux de purge soient plus performantes et plus robustes en conditions off-design, ces dernières n’ont jamais démarré lors des expérimentations même si les canaux ont été agrandis et multipliés. Ainsi, parmi tous les concepts essayés, celui qui démarre par survitesse et qui comporte des canaux de succion de couche limite à son col a donné les meilleurs résultats. Il est très stable et permet d’obtenir un ratio de pression statique de 4.25 et un recouvrement de pression totale de 89 %, pour une efficacité isentropique de 92 % à un nombre de Mach tangentiel de 2. Par contre, il est à noter qu’il n’a pas été possible de mesurer la pression totale. Elle a plutôt été estimée à partir des images de strioscopie tirées lors des essais. Comme on ne dispose pas d’une structure permettant d’essayer le compresseur rotatif à Mach 2, il a fallu approximer l’influence de l’accélération centrifuge sur l’écoulement de la cascade et trouver un moyen d’intégrer le nouvel aubage à la roue. Un modèle permettant d’estimer les paramètres d’une couche limite se développant sur une plaque plane en rotation a permis de déduire que l’accélération transverse n’aurait qu’un effet légèrement favorable, puisqu’il permet d’amincir l’épaisseur de déplacement, réduisant ainsi les risques d’interaction en la couche limite et les chocs. Finalement, les canaux de succion de couche limite du compresseur pourraient permettre d’alimenter un système de refroidissement qui limiterait la température à la jante à 820 K. Le R4E pourrait devenir l’avenir des systèmes de régénération électrique pour les véhicules hybrides. Il serait aussi intéressant pour une utilisation dans les petites centrales thermiques des régions éloignées. Ce grand potentiel d’utilisation provient de la grande densité de puissance du moteur, de sa simplicité et de son très faible coût de fabrication et de maintenance.

Un modèle de cavitation à poche partielle a été développé et implémenté dans un code potentiel bidimensionnel stationnaire puis dans un code potentiel tridimensionnel instationnaire. L'originalité de la méthode est d'utiliser la technique des vitesses de transpiration pour simuler la présence de la poche de cavitation. Cette méthode permet une fermeture naturelle de la poche sans avoir à imposer un quelconque artifice dans le modèle ou dans la méthode de calcul. Adaptable à tout solveur de modèle d'écoulement, le modèle couplé aux codes potentiels permet, grâce à cette technique, une simulation raisonnablement rapide de l'écoulement cavitant sur un hydrofoil ou une hélice en régime d'écoulement stationnaire ou instationnaire. Le modèle a été validé par comparaison avec des essais en régime stationnaire sur un profil bidimensionnel. On retrouve numériquement les longueurs de poche observées lors des essais ainsi que les pressions mesurées. Quant à l'écoulement tridimensionnel, des résultats disponibles dans la littérature n'ont permis que la validation qualitative au niveau de la longueur de poche. Une campagne d'essais en collaboration avec le BSHC pour mesurer la portance et la traînée d'un hydrofoil a alors été effectuée avec pour objectif la validation tridimensionnelle quantitative. Les résultats obtenus montrent une bonne concordance entre les mesures expérimentales et les résultats numériques, ce qui permet de conclure de façon satisfaisante quant à la validation du modèle. Finalement, des résultats de simulations instationnaires sur des hélices en régime cavitant sont présentés. Les résultats numériques montrent que la cavitation à poche partielle affecte très peu les performances hydrodynamiques de l'hélice.

.

Cette étude est consacrée à l'analyse par simulations numériques instationnaires tridimensionnelles LES de la transition laminaire-turbulent dans un canal à injection pariétale représentatif de l'écoulement dans les chambres de moteurs à propergol solide. Pour déclencher le phénomène de transition vers la turbulence avec la méthode LES, il a été nécessaire d'introduire des fluctuations dans l'écoulement. Cette étude se focalise sur la recherche d'une méthodologie de calcul permettant de générer des perturbations dans l'écoulement. Une première méthode, basée sur le comportement du schéma numérique du code de calcul CEDRE de l'ONERA, a été utilisée pour générer des fluctuations dans l'écoulement. Une deuxième méthode prenant en compte les données expérimentales relevées sur le montage en gaz froid VECLA de l'ONERA a ensuite 't' développée. Cette dernière utilise une modélisation stochastique dans le cadre de l'equation de Langevin. Après implantation de cette méthode de calcul stochastique dans le code CEDRE, une étude spécifique a été menée pour analyser le comportement du modèle en relation avec la problématique de l'écoulement transitionnel dans le canal avec injection pariétale. Les simulations numériques instationnaires LES de l'écoulement en transition dans le montage VECLA ont ensuite été réalisées avec le code de calcul CEDRE. Les résultats ont été confrontés aux données expérimentales ainsi qu'a une simulation stationnaire RANS de l'écoulement. Les résultats mettent en évidence l'intérêt de la modélisation stochastique de Langevin sur la description du champ fluctuant dans le canal avec injection pari'etale.

Le diagnostic du SAGE Loire au niveau de la plaine du Forez (42) relate une dégradation du potentiel hydrogéologique et une pollution des nappes superficielles par les nitrates et les pesticides. Ces problèmes nécessitent une gestion durable et équilibrée de la ressource en eau. Une partie de la demande en eau est assurée par la ressource souterraine. Les eaux souterraines exploitées sont : les eaux des formations alluviales superficielles pour l'irrigation et l'AEP ; et les eaux minérales des niveaux aquifères profonds du bassin sédimentaire sous-jacents aux alluvions pour le thermalisme et de mise en bouteille.La synthèse des connaissances disponibles a révélé des insuffisances dans la description des aquifères superficiels. En effet, elles ne permettent pas d'appréhender comment fluctue la ressource dans le temps et ne prennent pas en compte l'origine de la recharge des aquifères superficiels. Ainsi, le travail de thèse se focalise sur ces deux points particuliers :- l'identification du fonctionnement transitoire des aquifères est abordée par caractérisation physique des écoulements s'appuyant sur l'analyse de 4 points de suivi piézométrique permanent et la modélisation numérique des écoulements permettant de reproduire les fluctuations piézométriques observées.- l'origine de la recharge des aquifères alluviaux est abordée par caractérisation géochimique des écoulements (analyses en ions majeurs et isotopes δ180, δ2H, δ13C et A14C). Cette approche a mis en évidence des apports d'eaux profondes alimentant localement les aquifères superficiels.La méthodologie proposée, couplée physique et géochimie, permet de mieux préciser la vulnérabilité des aquifères investigués, en vue d'une meilleure gestion future.

Les calculs Monte-Carlo de criticité permettent d'estimer le facteur de multiplication effectif ("k-effectif") d'un système fissile au cours d'itérations simulant la propagation d'une population de neutrons, formant une chaîne de Markov. L'initialisation arbitraire de la population des neutrons simulés peut biaiser fortement l'estimation du k-effectif du système, défini comme la moyenne de la séquence des k-effectifs estimés à chaque itération. Un modèle simplifié de cette séquence de k-effectifs d'étapes est établi à partir du contexte technique d'exploitation industrielle des calculs Monte-Carlo de criticité. Des tests statistiques, inspirés des propriétés du pont brownien, sont construits pour discriminer la stationnarité de la séquence des k-effectifs d'étapes. Le régime transitoire initial éventuellement détecté est alors supprimé pour améliorer l'estimation du k-effectif du système. Les différentes déclinaisons de cette méthodologie sont détaillées puis comparées, d'une part sur un plan d'expériences représentatif des calculs Monte-Carlo de criticité, et d'autre part sur des calculs réels de configurations de criticité. Finalement, les performances observées sur ces applications permettent d'envisager une exploitation pertinente dans les calculs Monte-Carlo de criticité industriels.

Les émissions de méthane (CH[indice inférieur 4]), gaz à effet de serre provoquant le réchauffement climatique doivent être contrôlées. Les biofiltres peuvent être utilisés pour atteindre cet objectif. Les émissions de CH[indice inférieur 4] issues des industries agroalimentaires ou du traitement des eaux peuvent être accompagnées de vapeurs d’alcool. La présence simultanée de CH[indice inférieur 4], polluant à limitation par transfert de masse et d’alcool, polluant à limitation cinétique dans un mélange gazeux peut induire des limitations dans le biofiltre. L’objectif principal de cette recherche est l’évaluation des limitations dans un biofiltre traitant le CH[indice inférieur 4] en présence ou en absence de vapeur d’alcool en régime permanent ou transitoire. Dans un premier temps, une revue de littérature s’est penchée sur les limitations basées sur le transfert de masse et la cinétique lors de l’enlèvement de polluants organiques dans un biofiltre. Par la suite, l’élimination du CH[indice inférieur 4] a été effectuée dans un biofiltre afin d’évaluer l’influence de la concentration à l’entrée du biofiltre sur la performance du biofiltre. Une capacité d’élimination maximale de 45 g [indice supérieur -3] h[indice supérieur -1] a été obtenue pour une charge à l’entrée de 87 g [indice supérieur -3] h[indice supérieur -1] du biofiltre. Le biofiltre a toléré des charges par à-coups de CH[indice inférieur 4] de même que des privations de CH[indice inférieur 4] et de nutriments. Par conséquent, les comportements en régimes permanent et transitoire d’élimination du CH[indice inférieur 4] en présence de vapeurs d’éthanol ont été étudiés dans un biofiltre ayant un lit filtrant inorganique sous des temps de résidence en fût vide (EBRT) de 6, 3 et 1.5 minutes. L’ajout d’éthanol sur 3 cycles a été effectué en fonction des 3 EBRTs. Un EBRT de 6 min correspondant à des charges à l’entrée de CH4 et d’éthanol de 4.5 et de 132 g [indice supérieur -3] h[indice supérieur -1] a induit des limitations mineures en ce qui a trait à l’enlèvement du CH[indice inférieur 4] et de l’éthanol. En régime transitoire, la période de récupération après les 3 cycles a nécessité 10 à 25 jours. Ce délai est relié à la présence d’éthanol dans le lixiviat. Dans un dernier temps, deux biofiltres ayant un garnissage de pierres et un garnissage mixte ont été comparés pour l’enlèvement du CH[indice inférieur 4] et de l’éthanol présents dans un mélange gazeux en régime permanent. La section inférieure du biofiltre a permis l’élimination totale de l’éthanol. De plus, lors de l’élimination totale de l’éthanol dans la section inférieure du biofiltre, la production de dioxyde de carbone (CO[indice inférieur 2]) dépasse 16 g[indice supérieur -3] h[indice supérieur -1], pour des charges à l’entrée de CH[indice inférieur 4] et d’éthanol de 11 et 13 g m[indice supérieur -3] h[indice supérieur -1] respectivement. Par ailleurs, une concentration en éthanol dans le lixiviat excédant 2500 géthanol m[indice supérieur -3] lixiviat a été obtenue. Les biofiltres ont démontré une flexibilité pour des charges par à-coups d’éthanol suivies de périodes de carence. Le principal inconvénient du biofiltre à lit de pierres par rapport au biofiltre mixte est une perte de charge élevée dans la section inférieure du biofiltre. Une période de carence est un excellent moyen de contrer la perte de charge. / Abstract : Since methane (CH[subscript 4]) is a greenhouse gas with hazardous effects for global warming, every effort should be made to reduced methane emissions. Biofilters are potential candidates for CH[subscript 4] removal. In food and beverage industries as well as ethanol refineries, the feed of the biofilter might be a mixture of CH[subscript 4] emissions from wastewater treatment unit and ethanol emissions from other units. The presence of CH[subscript 4] as a mass transfer limited and ethanol vapor as a kinetic limited pollutant in a mixture can produce several limitations in a biofilter. The main objective of this research is to evaluate the limitations of CH[subscript 4] biofiltration or in the presence of ethanol vapors under steady and transient state conditions. First, a literature review was provided on mass transfer and kinetic limited organic pollutants removal in biofilters and the related limitations. Subsequently, the CH[subscript 4] elimination was assessed in a biofilter in order to evaluate the effect of CH[subscript 4] inlet concentration in the range of 1000 to 13000 ppmv and a gas flow rate of 3 L min[supercript -1] on the biofilter performance. A maximum CH[subscript 4] elimination capacity (ECmax) of 45 g m[superscript -3] h[superscript -1] was obtained for a CH[subscript 4] inlet load (IL) of 87 g m[superscript -3] h[superscript -1]. The biofilter tolerated CH[subscript 4] shock loads as well as different types of CH[subscript 4] and nutrient starvations. Subsequently, the steady state and dynamic behaviors of CH[subscript 4] elimination in the presence of ethanol vapor was studied in an inorganic bed biofilter with empty bed residence times (EBRTs) of 6, 3 and 1.5 min. Ethanol addition was performed in 3 cycles based on the EBRTs. An EBRT of 6 min with corresponding CH[subscript 4] and ethanol inlet loads of 132 and 4.5 gpollutant m[superscript -3] h[superscript -1] respectively, caused the least limitations for the simultaneous removal of CH[subscript 4] and ethanol in the biofilter. According to dynamic behavior of the biofilter, the recovery time after the three cycles took from 10 to 25 days. The delayed biofilter recovery was linked to the presence of ethanol in the liquid effluent. Finally, a stone-based bed and a hybrid packing biofilter were compared for CH[subscript 4] and ethanol removal in a mixture under steady and transient state conditions. Ethanol was completely removed in the bottom sections of both biofilters. A large carbon dioxide (CO[subscript 2]) production rate exceeding 18 g m[superscript -3] h[superscript -1] occurred in the bottom sections for CH[subscript 4] and ethanol inlet loads of 11 and 13 g m[superscript -3] h[superscript -1] respectively. In addition, an ethanol concentration in the leachate exceeding 2500 gethanol m[superscript -3] leachate was obtained for both biofilters. The biofilters were flexible to an ethanol shock load followed by a starvation period. The main drawback of the stone based bed biofilter compared to the hybrid packing biofilter was an excess pressure drop in the bottom section. Starvation was found an effective strategy for reducing the pressure drop.

La simulation des phénomènes électromagnétiques transitoires (EMT) est capitale dans plusieurs spécialisations de l'électrotechnique pour analyser, comprendre et prévoir divers comportements de réseaux électriques ou de dispositifs d'électronique de puissance commutée. Ces simulations sont typiquement réalisées sans contraintes temporelles : la simulation d'un événement peut être plus rapide ou plus lente que l'événement lui-même. Toutefois, dans un contexte de validation de systèmes de contrôle réel, les simulations doivent être effectuées en temps réel afin de soumettre l'équipement réel à des signaux réalistes avec le bon contenu spectral sur la plage fréquentielle requise. Une simulation opérant avec un pas de temps Tₛ est dite en temps réel si elle est en mesure de produire une nouvelle solution complète et d'effectuer tous les services nécessaires (entrées/sorties, acquisition, monitorage, etc.) en un laps de temps de Tₛ, et ce de manière maintenue et périodique. Les besoins de simulation temps réel ne cessent de croître pour étudier l'intégration des différents dispositifs d'électronique de puissance avec système de contrôle et protection qui prolifèrent sur les réseaux électriques et pour évaluer leurs interactions potentielles. Les réseaux à simuler sont de plus en plus vastes parce que ces dispositifs sont de plus en plus distribués spatialement. De plus, les modélisations sont de plus en plus complexes afin de représenter le plus fidèlement possible les phénomènes réels sur des plages de fréquence de plus en plus larges. Il existe donc un besoin criant d'augmenter la capacité de simulation temps réel. Pour y arriver, plusieurs groupes de recherche explorent l'utilisation d'unités de traitement alternatives telles les FPGA et GPU pour complémenter ou remplacer les processeurs standards tandis que d'autres explorent des modélisations et techniques pour augmenter la capacité de simulation des processeurs standards tels les équivalents fréquentiels, la simulation multitaux et la simulation hybride. Toutefois, la majorité de ces efforts vise l'augmentation des capacités en temps différé et non en temps réel. Afin de combler ce décalage, le présent ouvrage vise l'augmentation des capacités de simulation en temps réel des processeurs standards pour les grands et très grands réseaux électriques modernes c'est-à-dire où l'électronique de puissance foisonne et impose des pas de temps de plus en plus faibles. Deux avenues sont explorées : le parallélisme massif et l'augmentation de capacité par l'utilisation de simulation à taux multiples, qui mènera à la simulation hybride. En premier lieu, le parallélisme massif de la simulation temps réel est présenté. La puissance de calcul de superordinateurs n'est pas triviale à canaliser pour la simulation EMT temps réel à cause des difficultés de synchronisation. Une fois maîtrisé, le Superdome Flex de HPE a effectué la simulation en temps réel d'un réseau de 16432 nœuds électriques avec 372 cœurs de calcul opérant avec un pas de temps de 17 µs. Ces travaux sur les plateformes matérielles a également permis de découvrir les excellentes performances temps réel d'ordinateurs personnels ordinaires bas de gamme (RTPC) et très accessibles, ouvrant la voie à une démocratisation de la simulation temps réel. En second lieu, le cadre de simulation multitaux temps réel est présenté ainsi que les améliorations nécessaires à son bon fonctionnement. En simulant avec des taux plus lents que le taux de base EMT, il est possible d'augmenter substantiellement la capacité de calcul d'une plateforme donnée. Des simulations avec des pas de temps habituellement réservés aux FPGA (5 et 10 µs) sont effectuées sur de modestes RTPC En exploitant le cadre de simulation multitaux, il est ensuite montré comment la simulation hybride EMT/stabilité transitoire permet d'augmenter drastiquement (d'un facteur 23 à 10 µs et 1 ms) les capacités de simulation en temps réel d'un cœur de calcul. Toutefois, ce gain se fait au prix d'une perte de l'interaction harmonique entre le système EMT et celui de stabilité transitoire. L'implémentation temps réel de cette simulation hybride en temps réel est facilitée par l'utilisation d'une interface basée sur une modélisation hybride des lignes de transport et d'un protocole d'interaction parallèle, ce qui diffère de la méthode classique des images et des protocoles itératifs sophistiqués. Finalement, ces trois approches sont comparées : selon la situation, l'une de ces approches sera à préconiser selon l'ampleur de la simulation, la plage fréquentielle désirée et la quantité de ressources temps réel disponible. / The simulation of transient electromagnetic (EMT) phenomena is essential in several specializations of electrical engineering to analyze, understand and predict various behaviors of power systems and switched power electronic devices. These simulations are typically performed without time constraint: the simulation of an event can be faster or slower than the real event itself. However, in the context of validating real control systems, simulations must be performed in real-time in order to subject real equipment to realistic signals with the right spectral content over the required frequency range. A simulation operating with a time step Tₛ is said to be in real-time if it is able to produce a complete new solution and perform all the necessary services (inputs/outputs, acquisition, monitoring, etc.) in a period of time of Tₛ in a maintained and periodic manner. The need for real-time simulation continues to grow to study the integration of the various power electronic devices with control and protection systems that proliferate in power systems and to assess their potential interactions. The power systems to be simulated are increasingly vast because the spatial distribution of these devices is increasing. Moreover, model complexity is ever increasing in order to represent as faithfully as possible the real phenomena over increasingly wide frequency ranges. There is therefore a crying need to increase the capacity for real-time simulation. Several research groups are exploring the use of alternative processing units such as FPGAs and GPUs to complement or replace standard processors, while others are exploring models and techniques to increase the simulation capacity of standard processors such as frequency equivalents, multi-rate simulation and hybrid simulation. However, the majority of these efforts are aimed at increasing capacity in offline simulations, not for real time. In order to bridge this gap, this work aims to increase the real-time simulation capabilities of standard processors for large and very large modern power systems, i.e. where power electronics abound and impose smaller and smaller time steps. Two research avenues are explored: massive parallelism and increasing capacity through the use of multi-rate simulation, which will lead to hybrid simulation. First, the massive parallelism of real-time simulation is presented. The computing power of supercomputers is not trivial to harness for real-time EMT simulations because of synchronization difficulties. Once mastered, the HPE Superdome Flex performed real-time simulation of a network of 16,432 power nodes with 372 computer cores operating at a time step of 17 µs. This work on hardware platforms has also made it possible to discover the excellent real-time performance of ordinary low-end, but very accessible, personal computers (RTPC), paving the way for a democratization of real-time simulation. Secondly, the real-time multi-rate simulation framework is presented as well as the improvements necessary for its proper functioning. By simulating with rates slower than the base EMT rate, it is possible to substantially increase the computational capacity of a given platform. Simulations with time steps usually reserved for FPGAs (5 and 10 µs) are carried out on modest RTPCs. By exploiting the multi-rate simulation framework, it is then shown how the EMT / transient stability hybrid simulation can drastically increase (by a factor of 23 with 10 µs and 1 ms) the real-time simulation capacities of a computer core. However, this gain comes at the cost of a loss of the harmonic interaction between the EMT and transient stability systems. The real-time implementation of this real-time hybrid simulation is facilitated by the use of an interface based on hybrid transmission line modeling and a simple parallel interaction protocol, which differs from the classic equivalent image method and sophisticated iterative interaction protocols. Finally, these three approaches are compared: depending on the situation, one of these approaches will be more suited depending on the desired size of the simulation, the desired frequency range and the amount of available real-time resources.