Dans certaines simulations numériques exigeantes de mécanique des fluides, ilest nécessaire de simuler des écoulements multiphasiques impliquant de nombreuses contraintes simultanées : nombre de fluides important, évolutions compressibles à la fois isentropes et fortement choquées, équations d’états variables et contrastées, déformations importantes et transport surdes longues distances. Afin de remplir ces objectifs de manière robuste, il est nécessaire que la cohérence thermodynamique du schéma numérique soit vérifiée.Dans le premier chapitre, un schéma de type Lagrange plus projection est proposé pour la simulation d’écoulements diphasiques avec un modèle squelette à six équations et sans termes de dissipation. L’importance de la propriété de préservation des écoulements isentropiques est mise en évidence à l’aide d’une comparaison avec des résultats issus de la littérature pour le test deRansom. Ce chapitre souligne aussi certaines limitations de l’approche Lagrange plus projection pour simuler des modèles multiphasiques.Afin de pallier à ces limitations, une nouvelle procédure de dérivation est proposée afin de construire un schéma mimétique pour la simulation d’écoulements instationnaires compressibles dans un formalisme ALE direct (Arbitrary Lagrangian–Eulerian). La possibilité de choisir a prioriles degrés de liberté permet de s’inscrire dans une continuité avec les schémas historiques décalés, tout en imposant les conservations au niveau discret. L’équation de quantité de mouvement discrèteest obtenue par application d’un principe variationnel, assurant par construction la cohérence thermodynamique des efforts de pression. Cette approche est appliquée au cas d’écoulements monofluides comme preuve de concept au Chapitre 3, puis elle est étendue au cas d’écoulements à Nphasescompressibles au Chapitre 4. Des tests mono et multiphasiques montrent un comportement satisfaisant en terme de conservativité, versatilité aux mouvements de grilles et robustesse. / In some highly demanding fluid dynamics simulations, it appears necessary tosimulate multiphase flows involving numerous constraints at the same time : large numbers of fluids, both isentropic and strongly shocked compressible evolution, highly variable and contrasted equations of state, large deformations, and transport over large distances. Fulfilling such a challengein a robust and tractable way demands that thermodynamic consistency of the numerical scheme be carefully ensured.In the first chapter, a Lagrange plus remap scheme is proposed for the simulation of two-phase flows with a dissipation-free six-equation bakcbone model. The importance of the property of isentropic flow preservation is highlighted with a comparison with Ransom test results fromthe literature. This chapter also also point out certain limitations of the Lagrange plus remap approach for multiphase simulations.In order to overcome these limitations, a novel derivation procedure is proposed to construct a mimetic scheme for the simulation of unsteady and compressible flows in a direct ALE (ArbitraryLagrangian-Eulerian) formalism. The possibility to choose a priori the degrees of freedom allows to obtain a continuity with historical staggered scheme, while imposing conservativity at discretelevel. The discrete momentum evolution equation is obtained by application of a variational principle, thus natively ensuring the thermodynamic consistency of pressure efforts. This approach is applied to single-fluid flows as a proof of concept in Chapter 3, then it is extended to N-phasecompressible flows in Chapter 4. Single- and multi-phase tests show satisfactory behavior in terms on conservation, versatility to grid motions, and robustness.

This thesis deals with aspects of using the common, universal relational table editor as a simple information system, which would be fully independent on the underlying database system and partially independent on the underlying's database schema. One part of this thesis deals with potential of using such universal information system to create a framework to allow fast and easy development of small and medium information systems. The practical part of this thesis is the application, which implements the basics of simple relational table editor, which is fully independent on the underlying database provider and schema and servers as the demonstrative table editor.

Thèse diffusée initialement dans le cadre d'un projet pilote des Presses de l'Université de Montréal/Centre d'édition numérique UdeM (1997-2008) avec l'autorisation de l'auteur. / Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

L’imagerie du sous-sol par équations d’onde est une application de l’ingénierie pétrolière qui mobilise des ressources de calcul très importantes. On dispose aujourd’hui de calculateurs puissants qui rendent accessible l’imagerie de régions complexes mais des progrès sont encore nécessaires pour réduire les coûts de calcul et améliorer la qualité des simulations. Les méthodes utilisées aujourd’hui ne permettent toujours pas d’imager correctement des régions très hétérogènes 3D parce qu’elles sont trop coûteuses et /ou pas assez précises. Les méthodes d’éléments finis sont reconnues pour leur efficacité à produire des simulations de qualité dans des milieux hétérogènes. Dans cette thèse, on a fait le choix d’utiliser une méthode de Galerkine discontinue (DG) d’ordre élevé à flux centrés pour résoudre l’équation des ondes acoustiques et on développe un schéma d’ordre élevé pour l’intégration en temps qui peut se coupler avec la technique de discrétisation en espace, sans générer des coûts de calcul plus élevés qu’avec le schéma d’ordre deux Leap-Frog qui est le plus couramment employé. Le nouveau schéma est comparé au schéma d’ordre élevé ADER qui s’avère plus coûteux car il requiert un plus grand nombre d’opérations pour un niveau de précision fixé. De plus, le schéma ADER utilise plus de mémoire, ce qui joue aussi en faveur du nouveau schéma car la production d’images du sous-sol consomme beaucoup de mémoire et justifie de développer des méthodes numériques qui utilisent la mémoire au minimum. On analyse également la précision des deux schémas intégrés dans un code industriel et appliqués à des cas test réalistes. On met en évidence des phénomènes de pollution numériques liés à la mise en oeuvre d'une source ponctuelle dans le schéma DG et on montre qu'on peut éliminer ces ondes parasites en introduisant un terme de pénalisation non dissipatif dans la formulation DG. On finit cette thèse en discutant les difficultés engendrées par l'utilisation de schémas numériques dans un contexte industriel, et en particulier l'effet des calculs en simple précision. / Oil engineering uses a wide variety of technologies including imaging wave equation which involves very large computing resources. Very powerful computers are now available that make imaging of complex areas possible, but further progress is needed both to reduce the computational cost and improve the simulation accuracy. The current methods still do not allow to image properly heterogeneous 3D regions because they are too expensive and / or not accurate enough. Finite element methods turn out to be efficient for producing good simulations in heterogeneous media. In this thesis, we thus chose to use a high order Discontinuous Galerkin (DG) method based upon centered fluxes to solve the acoustic wave equation and developed a high-order scheme for time integration which can be coupled with the space discretization technique, without generating higher computational cost than the second-order Leap Frog scheme which is the most widely used . The new scheme is compared to the high order ADER scheme which is more expensive because it requires a larger number of computations for a fixed level of accuracy. In addition, the ADER scheme uses more memory, which also works in favor of the new scheme since producing subsurface images consumes lots of memory and justifies the development of low-memory numerical methods. The accuracy of both schemes is then analyzed when they are included in an industrial code and applied to realistic problems. The comparison highlights the phenomena of numerical pollution that occur when injecting a point source in the DG scheme and shows that spurious waves can be eliminated by introducing a non-dissipative penalty term in the DG formulation. This work ends by discussing the difficulties induced by using numerical methods in an industrial framework, and in particular the effect of single precision calculations.

Lors de la formation du point chaud dans une expérience de Fusion par Confinement Inertiel, le plasma au centre de la sphère de deutérium-tritium peut être loin de l'équilibre thermodynamique local. Dans la première partie, on décrit donc un modèle cinétique ionique de type Vlasov-Fokker-Planck susceptible de prendre en compte ces déséquilibres. Après avoir rappelé les grandes étapes pour résoudre numériquement le système obtenu, on introduit la notion de moyenne entropique pour définir un nouveau schéma numérique traitant les collisions ion-électron homogènes en espace. Ce schéma est conservatif, stable et entropique sous un critère de type CFL dans sa version explicite. Dans sa version semi-implicite, on établit que ce schéma conserve l'équilibre thermodynamique. Le temps de calcul pour résoudre les équations cinétiques étant très important, il est nécessaire d'étudier la possibilité de ne résoudre ces équations que là où c'est nécessaire c'est à dire principalement au centre de la sphère de deutérium-tritium. Dans la seconde partie, on propose donc une technique de couplage cinétique-fluide, la formation du point chaud étant traitée avec le modèle cinétique, le reste avec les équations d'Euler à deux températures (températures ionique et électronique). Les ions deutérium et tritium pouvant ne pas être à l'équilibre thermodynamique, on s'est ensuite posé la question de la validité des formules analytiques donnant le taux de réaction nucléaire, formules établies en supposant que le plasma est à l'équilibre thermodynamique. Dans la troisième partie, on propose donc une méthode de type Monte-Carlo pour résoudre numériquement les équations cinétiques de type Boltzmann qui décrivent les réactions de fusion thermonucléaire et on montre qu'effectivement, les déséquilibres thermodynamiques rencontrés lors de la formation du point chaud peuvent invalider les formules usuelles.

Ce travail présente quelques résultats concernant le comportement en temps fini et en temps long de méthodes numériques pour des équations stochastiques. On s'intéresse d'abord aux équations différentielles stochastiques de Langevin et de Langevin amorti. On montre un résultat concernant l'analyse d'erreur faible rétrograde de ses équations par des schémas numériques implicites. En particulier, on montre que l'erreur entre le générateur associé au schéma numérique et la solution d'une équation de Kolmogorov modifiée est d'ordre élevé par rapport au pas de discrétisation. On montre aussi que la dynamique associée au schéma numérique est exponentiellement mélangeante. Dans un deuxième temps, on étudie le comportement en temps long d'une discrétisation en temps et en espace d'une EDPS semi-linéaire avec un bruit blanc additif, qui possède une unique mesure invariante . On considère une discrétisation en temps par un schéma d'Euler et en espace par une méthode des éléments finis. On montre que la moyenne, par rapport aux lois invariantes (qui n'est pas forcément unique) associées à l'approximation, par des fonctions tests suffisamment régulières est proche de la quantité correspondante pour . Plus précisément, on étudie la vitesse de convergence par rapport aux différents paramètres de discrétisation. Enfin, on s'intéresse à une EDPS semi-linéaire avec un bruit blanc additif dont le terme non-linéaire est un polynôme. On étudie la convergence au sens faible d'une approximation en temps par un schéma de splitting implicite.

Ce mémoire résume les travaux que nous avons réalisés de 1995 à 2010. Ces travaux ont eu pour thème l'étude et la discrétisation, d'une part, de modèles cinétiques de type Fokker-Planck ou de type Boltzmann semi-classiques et, d'autre part, de modèles fluides de type Euler ou de type Navier-Stokes à bas nombre de Mach. L'équation de Fokker-Planck étudiée modélise les collisions entre ions et électrons dans un plasma chaud, et concerne ici la fusion par confinement inertiel. Les équations de Boltzmann semi-classiques étudiées sont de deux types. Le premier type modélise la réaction de fusion thermonucléaire entre un ion deuterium et un ion tritium donnant une particule alpha et un neutron, et concerne également ici la fusion par confinement inertiel. Le deuxième type - connu sous le nom d'équations de Wang-Chang & Uhlenbeck - modélise ici les transitions d'énergie quantique dans les couches électroniques d'atomes d'uranium et de fer provoquées par les collisions entre ces mêmes atomes au sein du procédé SILVA de Séparation Isotopique par Laser sur Vapeur Atomique. Nous avons étudié les propriétés de base de ces deux types d'équations de Boltzmann semi-classiques, et, dans le cas des équations de Wang-Chang & Uhlenbeck, nous avons proposé un algorithme de couplage cinétique-fluide. L'étude de cet algorithme nous a incité à étudier la notion de relaxation dans un mélange binaire de gaz et de fluides non-miscibles, et à souligner les points communs de cette approche avec la théorie cinétique standard. L'étude de modèles moyennés pour des mélanges de fluides non-miscibles nous a amené à proposer et à discrétiser un modèle sans ondes acoustiques modélisant la déformation d'une interface entre deux fluides non-miscibles provoquée par de forts gradients thermiques à bas nombre de Mach. Puis, afin d'améliorer la précision des calculs tout en en maîtrisant le coût, nous avons également étudié la possibilité de résoudre sur un maillage dynamique de type AMR un modèle simplifié de déformation d'interface. Ces études à bas nombre de Mach nous ont également incités à analyser sur maillage cartésien le mauvais comportement à bas nombre de Mach des schémas de type Godunov appliqués au système d'Euler compressible. Enfin, nous avons justifié l'algorithme LBM dans le cas de l'équation de la chaleur.

Cette thèse présente la construction de schémas distribuant le résidu (RD) d'ordre très élevés, pour la discrétisation d'équations d'advection-diffusion multidimensionnelles et stationnaires sur maillages non structurés. Des schémas linéaires ainsi que des schémas non linéaires sont considérés. Une approximation de la solution polynomiale par morceaux et continue sur chaque élément est adoptée, de plus une procédure de reconstruction du gradient que celle de la solution numérique est utilisée afin d'avoir une représentation continue de la solution numérique et de son gradient. Il est montré que le gradient doit être reconstruit avec la même précision de la solution, sans quoi la précision formel du schéma numérique est perdue dans les cas où les effets de diffusion prévalent sur les effets d'advection, et aussi quand l'advection et la diffusion sont également importants. Ensuite, la méthode est étendue à des systèmes d'équations, en particulier aux équations de Navier-Stokes et aux équations RANS. La précision, l'efficacité et la robustesse du solveur RD implicite sont démontrées sur plusieurs cas tests. / The construction of compact high-order Residual Distribution schemes for the discretizationof steady multidimensional advection-diffusion problems on unstructuredgrids is presented. Linear and non-linear scheme are considered. A piecewise continuouspolynomial approximation of the solution is adopted and a gradient reconstructionprocedure is used in order to have a continuous representation of both thenumerical solution and its gradient. It is shown that the gradient must be reconstructedwith the same accuracy of the solution, otherwise the formal accuracy ofthe numerical scheme is lost in applications in which diffusive effects prevail overthe advective ones, and when advection and diffusion are equally important. Thenthe method is extended to systems of equations, with particular emphasis on theNavier-Stokes and RANS equations. The accuracy, efficiency, and robustness of theimplicit RD solver is demonstrated using a variety of challenging aerodynamic testproblems.

Les travaux de cette thèse concerne l’étude d’une méthode élément finis d’ordre spatial élevé que l’on peut assimilé à une extension du schéma de Yee. On parle alors de méthode différences finies d’ordre élevé. Après avoir donné, dans un premier chapitre, un historique non exhaustif des principales méthodes utilisées pour résoudre les équations de Maxwell dans le cadre de problèmes de CEM et montré l’ intérêt de disposer d’un solveur de type "différences finies d’ ordre élevé", nous présentons dans un deuxième chapitre le principe de la méthode. Nous donnons pour cela les caractéristiques du schéma spatial et temporel en précisant les conditions de stabilité de la méthode. En outre, dans une étude purement numérique, nous étudions la convergence du schéma. On se focalise ensuite sur la possibilité d’utiliser des ordres spatiaux variable par cellules dans chaque direction de l’espace. Des comparaisons avec le schéma de Yee et un schéma de Galerkin Discontinu particulier sont ensuite effectuées pour montrer les gains en coûts calcul et mémoire et donc l’intérêt de notre approche par rapport aux deux autres. Dans un troisième chapitre, nous nous intéressons à l’étude de modèles physiques indispensable au traitement d’un problème de CEM. Pour cela, nous nous focalisons particulièrement sur un modèle de fil mince oblique, des modèles de matériaux volumiques et minces et enfin sur la prise en compte de sol parfaitement métallique dans une agression de type onde plane. Chaque modèle est détaillé et validé par comparaison avec des solutions analytiques ou résultant de la littérature, sur des exemples canoniques. Le quatrième chapitre est dédié à une technique d’hybridation entre notre méthode et une approche Galerkin Discontinu en vue de traiter des géométries possédant des courbures. Nous donnons pour cela une stratégie d’hybridation basée sur l’échange de flux qui garantie au niveau continue la conservation d’une énergie. Nous présentons ensuite quelques exemples montrant la validité de notre approche dans une stratégie multi-domaines/multi-méthodes que nous précisons. Enfin le dernier chapitre de cette thèse concerne l’exploitation de notre méthode sur des cas industriels en comparaisons avec d’autres méthodes ou des résultats expérimentaux. / This thesis is about the study of a high spatial finite element method whichcan be assimilated at an extension of the Yee schema. In the next, this method is also called high order finite difference method. In the first chapter, we give a non exhaustive recall of the major methods used to treat EMC problems and we show the necessity to have this kind of schema to simulate efficiently some EMC configurations. In the second chapter, the principle of the numerical method is presented and a stability condition is given. A numerical study analysis of the schema convergence is also done. Next, we show the interest to have the possibility to use local spatial order by cell in each direction of the computational domain. Some canonic examples are given to show the advantages interms of CPU time and memory storage of the method by comparison with Yee’s scheme and DG approach. In the third chapter, we define and validate on several examples,some physical models as thin wire, materials and perfectly metallic ground in presence of a plane wave, to have the possibility to treat EMC problems. The fourth chapter is about a hybridization strategy between our high order FDTD method and a DG schema.We focalize our study on a hybrid method which provides an energy conservation of the continuous problem. A numerical example is given to validate the method. Finally, in the last chapter, we present some simulations on industrial problems to show the possibility of the method to treat realistic EMC problems.

The aim of this thesis is the development project of the firm TREFAL, Spol. s.r.o. The object SOO2 - Administrative building - that is a constituent of the area of TREFAL, Spol. s.r.o. is discussed topic of this thesis. The main focus is the construction development of the monolithic supporting construction objects. The thesis further depicts and investigates accompanying and overall technical report, an implementation study of the main technological stages, situational designs, a solution of the broader traffic relations including a proposal for the supplementation of the construction development, BOZP plan, budget and timeline, a proposal for the new local communication and following textual and design documentation according to the thesis assignments.