Dans cette thèse, nous présentons une nouvelle approche pour la simulation numérique des phénomènes électroconvectifs et électro-therrno-convectifs. La principale difficulté réside dans la détermination du champ électrique et de la distribution de densité volurnique de charges électriques. Dans cette approche, des schémas de type TVD (Total Variation Dirninishing) et IDC (ImprovedDeferred Correction) sont utilisées dans la discrétisation des flux convectifs et diffusifs par la méthode des volumes finis. La première partie de cette thèse présente certains aspects numériques liés à l'implémentation de ces schémas. Une approche unifiée pour les schémas convectifs TVD de type borné à haute résolution est présentée et diverses fonctions limiteur sont comparées. Dans une deuxième partie, l'électro-convection entre deux plaques parallèles est simulée. La méthodologie a étéévaluée et validée par la détermination des Critères de stabilité linéaire et non linéaire. Les différents scenarii d'évolution du développement de cette instabilité électroconvective vers l'état chaotique ont été définis. L'effet du mécanisme de diffusion la densité volumique de charge sur la boucle d'hystérésis et sur la structure de l'écoulement est étudié. L'influence du rapport d'aspect de la cavité est analysé. Enfin dans une dernière partie, nous étudions l'électro-thermo-convection lorsque le fluide est soumis simultanément à une injection unipolaire et à un gradient thermique. L'augmentation des transferts de chaleur a été caractérisée. / In this thesis, a numerical approach is presented to simulate the electro- and electro-thermo convection in dielectric liquids. The total variation diminishing (TVD) scheme and improved deferred correction (IDC) scheme are used to compute the convective and diffusive respectively. The aim of TVD scheme is to avoid non-physical oscillations and to capture high gradient of charge density. Some fundarnental aspects related to TVD and LDC schemes are investigated firstly. A unified approach for TVD schemes is explained and various limiter functions are compared. The connection among three methods for diffusive flux computation has been revealed. The original IDC scheme is improved by the application of 2nd order gradient evaluation method.The electro-convection between two parallel plates is then simulated. The methodology was assessed by the determination of the linear and nonlinear stability criterion. By continuously increasing the driving parameter, the successive instabilities and route to chaotic state has been defined. The effects of the diffusion mechanism for the charge density and vertical walls on the hysteresis 100p and the structure are also investigated. The last part is to simuiate electro-thermo-convection when injection and thermal gradient are simultaneously applied. Our solver was verified with a stationary and an overstable stability problem.The case that both heating and injection are from a bottom electrode has been analyzed in details. The neutral stabiliïy curve was reproduced. The existence of nonlinear phenornenon and the structure are highlighted.

Plusieurs types d’applications industrielles complexes, relèvent du domaine multidisciplinaire de l’Electro-Hydro-Dynamique (EHD) où les interactions entre des particules chargées et des particules neutres sont étudiées dans le contexte couplé de la dynamique des fluides et de l’électrostatique. Dans cette thèse, nous avons étudié par voie de simulation numérique certains phénomènes Electro-Hydro-Dynamiques comme l’injection unipolaire, le phénomène de conduction dans les liquides peu conducteurs et le contrôle d’écoulement avec des actionneurs plasma à barrières diélectriques (DBD). La résolution de tels systèmes physiques complexes exige des ressources de calculs importantes ainsi que des solveurs CFD parallèles dans la mesure où ces modèles EHD sont mathématiquement raides et très consommateurs en temps de calculs en raison des gammes d’échelles de temps et d’espace impliquées. Cette thèse vise à accroitre les capacités de simulations numériques du groupe Electro-Fluido-Dynamique de l’Institut Pprime en développant un solveur parallèle haute performance basé sur des modèles EHD avancés. Dans une première partie de cette thèse, la parallélisation de notre solveur EHD a été réalisée avec des protocoles MPI avancés comme la topologie Cartésienne et les Inter-communicateurs. En particulier, une stratégie spécifique a été conçue pour prendre en compte la caractéristique multi-blocs structurés du code. La nouvelle version parallèle du code a été entièrement validée au travers de plusieurs benchmarks. Les tests de scalabilité menés sur notre cluster de 1200 cœurs ont montré d’excellentes performances. La deuxième partie de cette thèse est consacrée à la simulation numérique de plusieurs écoulements EHD typiques. Nous nous sommes intéressés entre autres à l’électroconvection induite par l'injection unipolaire entre deux électrodes plates parallèles, à l’étude des panaches électroconvectifs dans une configuration d'électrodes lame-plan, au mécanisme de conduction basé sur la dissociation de molécules neutres d'un liquide faiblement conducteur. Certains de ces nouveaux résultats ont été validés avec des simulations numériques entreprises avec le code commercial Comsol. Enfin, le contrôle d’écoulements grâce à un actionneur DBD a été simulé à l’aide du modèle Suzen-Huang dans diverses configurations. Les effets de l’épaisseur du diélectrique, de l’espacement inter-électrodes, de la fréquence de la tension appliquée et sa forme d’onde, sur la vitesse maximale du vent ionique induit ainsi que sur la force électrique moyenne ont été étudiés. / Several intricately coupled applications of modern industries fall under the multi-disciplinary domain of Electrohydrodynamics (EHD), where the interactions among charged and neutral particles are studied in context of both fluid dynamics and electrostatics together. The charge particles in fluids are generated with various physical mechanisms, and they move under the influence of external electric field and the fluid velocity. Generally, with sufficient electric force magnitudes, momentum transfer occurs from the charged species to the neutral particles also. This coupled system is solved with the Maxwell equations, charge transport equations and Navier-Stokes equations simulated sequentially in a common time loop. The charge transport is solved considering convection, diffusion, source terms and other relevant mechanisms for species. Then, the bulk fluid motion is simulated considering the induced electric force as a source term in the Navier-Stokes equations, thus, coupling the electrostatic system with the fluid. In this thesis, we numerically investigated some EHD phenomena like unipolar injection, conduction phenomenon in weakly conducting liquids and flow control with dielectric barrier discharge (DBD) plasma actuators.Solving such complex physical systems numerically requires high-end computing resources and parallel CFD solvers, as these large EHD models are mathematically stiff and highly time consuming due to the range of time and length scales involved. This thesis contributes towards advancing the capability of numerical simulations carried out within the EFD group at Institut Pprime by developing a high performance parallel solver with advanced EHD models. Being the most popular and specific technology, developed for the distributed memory platforms, Message Passing Interface (MPI) was used to parallelize our multi-block structured EHD solver. In the first part the parallelization of our numerical EHD solver with advanced MPI protocols such as Cartesian topology and Inter-Communicators is undertaken. In particular a specific strategy has been designed and detailed to account for the multi-block structured grids feature of the code. The parallel code has been fully validated through several benchmarks, and scalability tests carried out on up to 1200 cores on our local cluster showed excellent parallel speed-ups with our approach. A trustworthy database containing all these validation tests carried out on multiple cores is provided to assist in future developments. The second part of this thesis deals with the numerical simulations of several typical EHD flows. We have examined three-dimensional electroconvection induced by unipolar injection between two planar-parallel electrodes. Unsteady hexagonal cells were observed in our study. 3D flow phenomenon with electro-convective plumes was also studied in the blade-plane electrode configuration considering both autonomous and non-autonomous injection laws. Conduction mechanism based on the dissociation of neutral molecules of a weakly conductive liquid has been successfully simulated. Our results have been validated with some numerical computations undertaken with the commercial code Comsol. Physical implications of Robin boundary condition and Onsager effect on the charge species were highlighted in electro-conduction in a rectangular channel. Finally, flow control using Dielectric Barrier Discharge plasma actuator has been simulated using the Suzen-Huang model. Impacts of dielectric thickness, gap between the electrodes, frequency and waveform of applied voltage etc. were investigated in terms of their effect on the induced maximum ionic wind velocity and average body force. Flow control simulations with backward facing step showed that a laminar flow separation could be drastically controlled by placing the actuator at the tip of the step with both electrodes perpendicular to each other.