L’effet des inclusions sur la turbulence de l’écoulement est un élément clé à comprendre afin de maîtriser le transport de milieux dispersés, dans le domaine du génie pétrolier, environnemental, agroalimentaire, génie de la réaction chimique ou transformation du solide. Les expériences de Matas et al. (PRL, 2003) ont mis en évidence un effet non monotone des particules isodenses (de densité égale à celle du fluide) sur la transition laminaire-turbulent, cet effet dépendant de la taille des particules et de leur concentration dans la suspension. Une petite quantité de particules de taille finie s’est avérée suffisante pour diminuer considérablement le seuil de transition laminaire turbulent. Nous avons utilisé des simulations numériques, basées sur une approche de type “Force Coupling Method” afin de comprendre cet effet. Les domaines de simulations étaient choisis pour accommoder le minimum de structures cohérentes suffisantes pour entretenir la turbulence. Nous avons particulièrement étudié la corrélation entre le comportement instationnaire de l’écoulement et la distribution instantanée de particules, en fonction de la configuration de l’écoulement (Couette plan ou écoulement en canal), de la forme des particules ainsi que leur inertie et concentration. Dans un écoulement de Couette plan turbulent, la contrainte pariétale est augmentée en présence des particules. Les profiles (dans la direction normale aux parois) de vitesse moyenne et des contraintes de Reynolds ne sont pas significativement modifiés en présence des particules, si la viscosité du fluide est remplacée par la viscosité effective de la suspension dans le calcul du nombre de Reynolds de l’écoulement. Par contre l’analyse temporelle et modale des fluctuations de l’écoulement suggère que les particules modifient légèrement le cycle de régénération de la turbulence, à travers l’augmentation d’énergie à petites échelles. En effet, la forme des streaks et le caractère intermittent de l’écoulement sont impactés par la présence des particules, surtout quand elles sont inertielles (de densité supérieure à celle du fluide). Ces résultats ont été publiés dans le journal Physical Review F (Wang et al., 2017). En outre, nous avons montré qu’à fraction volumique égale, les propriétés d’écoulement turbulent des suspensions de particules sphéroïdales de rapport de taille compris entre 0.5 et 2, sont similaires à celles des suspensions de particules sphériques. Le transfert de particules entre les différentes structures cohérentes de l’écoulement est analysé à la fin de la thèse. Néanmoins dans un écoulement en canal, les particules iso denses augmentent l’intensité des contraintes de Reynolds dans le plan transverse. Nous montrons que par leur concentration préférentielle dans les structures cohérentes à côté des parois (les éjections), elles influencent significativement le cycle de régénération en agissant sur tous les processus à la fois linéaires et non linéaires du cycle: la formation des streaks, puis leur rupture et la régénération des vortex alignés avec l’écoulement. La diminution du seuil de transition est la conséquence directe de cette modulation du cycle. / The effect of particles on turbulence is a key phenomenon in many practical industrial applications encountered in petroleum engineering, chemical reactors and food or solid processing (transport of slurries in pipes, reactive fluidized beds, and pneumatic transport of particles), environmental engineering (such as sand storm and Particulate Matter (PM) Pollution), and biological fluid mechanics (e.g. drug delivery in blood flow and inhaled particles through the respiratory system). The experiments of Matas et al. (PRL, 2003) have highlighted the non-monotonous effect of neutrally buoyant particles on the laminar-turbulent flow transition, depending on the particle-to-pipe size ratio and on the suspension volumetric concentration. A small amount of finite size particles allowed sustaining the turbulent state and decreasing the transition threshold significantly. The complex mechanisms related to particle flow interactions are often difficult to elucidate experimentally. During the last 4 decades, direct numerical simulations have proven to be a powerful tool for understanding the features of single-phase turbulent flows. Currently, it starts to play an important role in the investigation of suspension flows as well. Almost a decade after the experiments of Matas et al. (PRL, 2003), particle-resolved numerical simulations are able to evidence that at moderate concentration, particles have a significant impact on the unsteady nature of the flow, enhancing the transverse turbulent stress components and modifying the flow vortical structures (Loisel et al. Phys. Fluids, 2013; Yu et al. Phys. Fluids, 2013; Lashgari et al. PRL, 2015). In this work, we use particleresolved numerical simulations to understand the effect of finite sized particles on wall-bounded (pressure-driven or plane Couette) turbulent flows, slightly above the laminar-turbulent transition limit. We find that in turbulent Couette flow, wall-normal profiles of the flow velocity and Reynolds stress components reveal that there is no significant difference between single phase and two-phase flows at equivalent effective Reynolds number, except that the wall shear stress is higher for the two-phase flow. At concentration up to 10%, neutrally buoyant spherical particles have a negligible effect on both the intensity and intermittency of the Reynolds stress. However temporal and modal analysis of flow fluctuations, suggest that besides increasing small scale perturbation due to their rigidity, particles have an effect on the regeneration cycle of turbulence (streak formation, streak breakdown and streamwise vortex regeneration). Indeed, the shape of the streaks and the intermittent character of the flow (amplitude and period of oscillation of the modal fluctuation energy) are all altered by the particle presence, and especially by the inertial particles (Wang et al. Phys. Rev. Fluid, 2017). When the particle shape deviates from sphericity (spheroids with aspect ratios ranging between 0.5 and 2), the features of turbulent suspension flow are not significantly impacted. The transfer of particles between different coherent structures (along the regeneration cycle period) is analyzed at the end of the thesis. Nevertheless in channel flow, neutrally-buoyant spherical particles have a drastic impact on the regeneration cycle of turbulence, decreasing thereby the transition threshold. Particles enhance the intensity of the Reynolds stress although the frequency of burst events is decreased. Particles enhance the lift-up effect and act continuously within the buffer layer. Moreover, they increase the vorticity stretching, leading to smaller and more numerous wavy streaks for suspension flows compared to the single-phase configuration.

In this thesis, we study the existence of random periodic solutions for both nonlinear dissipative stochastic functional differential equations (SFDEs) and semilinear nondissipative SFDEs in C([-r,0],R^d). Under some sufficient conditions for the existence of global semiflows for SFDEs, by using pullback-convergence technique to SFDE, we obtain a general theorem about the existence of random periodic solutions. By applying coupled forward-backward infinite horizon integral equations method, we perform the argument of the relative compactness of Wiener-Sobolev spaces in C([0,τ],C([-r,0]L²(Ω))) and the generalized Schauder's fixed point theorem to show the existence of random periodic solutions.

In this thesis, we study the existence of stationary solutions for two cases. One is for random difference equations. For this, we prove the existence and uniqueness of the stationary solutions in a finite-dimensional Euclidean space Rd by applying the coupling method. The other one is for semi linear stochastic evolution equations. For this case, we follows Mohammed, Zhang and Zhao [25]'s work. In an infinite-dimensional Hilbert space H, we release the Lipschitz constant restriction by using Arzela-Ascoli compactness argument. And we also weaken the globally bounded condition for F by applying forward and backward Gronwall inequality and coupling method.

Cryptography is the study of encryptying and decrypting messages and deciphering encrypted messages when the code is unknown. We consider Λπ(Δx, Δy) which is a count of how many ways a permutation satisfies a certain property. According to Hawkes and O'Connor, the distribution of Λπ(Δx, Δy) tends to a Poisson distribution with parameter ½ as m → ∞ for all Δx,Δy ∈ (Z/qZ)m - 0. We give a proof of this theorem using the Stein-Chen method: As qm approaches infinity, the distribution of Λπ(Δx, Δy) is approximately Poisson with parameter ½. Error bounds for this approximation are provided.

In this thesis, we study the existence of pathwise random periodic solutions to both the semilinear stochastic differential equations with linear multiplicative noise and the semilinear stochastic partial differential equations with linear multiplicative noise in a Hilbert space. We identify them as the solutions of coupled forward-backward infinite horizon stochastic integral equations in general cases, and then perform the argument of the relative compactness of Wiener-Sobolev spaces in C([0, T],L2Ω,Rd)) or C([0, T],L2(Ω x O)) and Schauder's fixed point theorem to show the existence of a solution of the coupled stochastic forward-backward infinite horizon integral equations.

Une méthode de couplage continu-discret a été développée pour simuler les mécanismes complexes d'endommagement de la silice soumise à un choc laser de haute puissance. Dans un premier temps, une classification des méthodes numériques existantes a été faite pour choisir celles les mieux adaptées à la simulation du comportement sous choc de la silice. Comme résultat de cette classification, deux méthodes ont été retenues: la méthode des éléments discrets (DEM) et la méthode des éléments naturels contraints (CNEM). Ces méthodes sont alors couplées en se basant sur la technique dite "Arlequin". Puis, un modèle numérique permettant de tenir compte des différents phénomènes qui caractérise le comportement de la silice sous haute pression a été développé. Pour bien caractériser les mécanismes de fissuration de la silice à l’échelle microscopique, un nouveau modèle de rupture a été développé dans ce travail. Finalement, ces deux modèles, modèle de comportement et modèle de rupture, ont été intégrés dans la méthode du couplage pour simuler d'un point de vue mécanique le choc laser sur un échantillon en silice. / A discrete-continuum coupling approach has been developed to simulate the laser-induced damage in silica glass. First, a classification of the different numerical methods has been performed to select the ones that best meet the objectives of this work. Acting upon this classification, the Discrete Element Method (DEM) and the Constrained Natural Element Method (CNEM) have been retained. Subsequently, a coupling approach between these methods has been proposed. This approach is based on the Arlequin technique. In the second part, a numerical model of the silica glass mechanical behavior has been developed to better characterize the silica glass response under highly dynamic loadings and particularly loading generated by a laser beam. To correctly characterize the silica glass cracking mechanisms, a new fracture model has been proposed in this work. Finally, all these developments have been used to simulate the laser-induced damage in silica glass.

Dans ce travail de thèse, nous nous sommes intéressés à la modélisation des phénomènes de diffraction d’ondes électromagnétiques par des objets à symétrie de révolution complexes et fortement hétérogènes. La méthode que nous développons ici consiste en un couplage entre équations aux dérivées partielles (EDP) et équations intégrales (EI). Cette idée est essentiellement connue pour avoir deux avantages. Le premier est que les hétérogénéités de l’objet sont prises en compte naturellement dans la formulation du problème. Le deuxième est dû à l’utilisation des équations intégrales qui donnent une représentation exacte des solutions dans le milieu extérieur en fonction des courants surfaciques. Le domaine de simulation peut ainsi être ramené à l’objet lui-même. L’utilisation de développements en séries de Fourier combinés à la propriété d’invariance par rotation de l’objet permet alors la réduction du problème global 3D à un ensemble dénombrable de problème 2D.L’étude de ces problèmes nous a conduit à décomposer notre analyse en plusieurs parties,chacune ayant à traiter une partie du problème complet ou les méthodes d’intégrations numériques. Ces dernières étant difficiles à réaliser dans le cas des équations intégrales.Nous avons tout d’abord étudié un problème de Maxwell intérieur pour lequel nous avons développé une nouvelle méthode d’éléments finis d’ordre élevé dont nous avons montré l’efficacité et la précision sur de multiples exemples. Puis, nous avons étudié le problème de diffraction d’ondes planes pour des objets parfaitement conducteurs. La méthode d’éléments finis de frontière employée est alors construite par extension de la méthode précédente via l’opérateur de trace tangentielle. En combinant ces deux études, nous avons résolu le problème couplé en introduisant la propriété de symétrie de révolution dans une formulation variationnelle bien choisie. Par construction, les éléments finis qui y sont utilisés sont alors naturellement adaptées. L’algorithme de parallélisation de la méthode de couplage est finalement présentée et des comparaisons entre notre code AxiMax et un code 3D sont illustrées. Dans tous les cas, nous montrons que la méthode d’éléments finis d’ordre élevé permet d’obtenir des résultats d’une grande précision en fonction de la qualité des paramètres de simulation. / In this thesis, we are interested in modeling diffraction of electromagnetic waves by axisymmetric and highly heterogeneous objects. Our method consists in a coupling between partial differential equations and integral equations. This idea is mainly interesting for two reasons : heterogeneities are handled naturally in the formulation and integral equations give an analytical representation of solutions outside the object based on surface currents.These advantages allow us to limit the domain of simulation to the object itself. In addition,using Fourier series combined with the rotational invariance property of the object, the 3D problem is reduced to a countable set of 2D problems. The study of these problems is split into several parts. Each part has to deal with aspecific problem as for example the numerical integration of singular integrals which is difficult to achieve. As a first step, we study time-harmonic Maxwell’s equations in a bounded domain for which we develop a new high-order finite element method and present its efficiency and accuracy on many examples. Secondly, we consider the diffraction of plane waves by perfect electric conductors to analyse integral equations for these kind of object.The boundary finite element method applied is defined by extension of the previous one via tangential trace operator. Then, we solve the coupled problem using a well chosen formulation based on the previous studies for which our finite element method is naturally adapted by construction. In order to evaluate its efficiency, a comparison is performed between our program « AxiMax » and one based on a purely 3D model. To conclude, in the last two chapters, we present the numerical integration method and the multi-processing algorithm developed in AxiMax. In all cases, we put forward the fact that our finite element method provides accurate results depending on the quality of the simulation parameters.

Ma thèse de doctorat se concentre principalement sur le comportement en temps long des processus de Markov, les inégalités fonctionnelles et les techniques relatives. Plus spécifiquement, Je vais présenter les taux de convergence sous-exponentielle explicites des processus de Markov dans deux approches : la méthode Meyn-Tweedie et l’hypocoercivité (faible). Le document se divise en trois parties. Dans la première partie, Je vais présenter quelques résultats importants et des connaissances connexes. D’abord, un aperçu de mon domaine de recherche sera donné. La convergence exponentielle (ou sous-exponentielle) des chaînes de Markov et des processus de Markov (à temps continu) est un sujet d’actualité dans la théorie des probabilité. La méthode traditionnelle développée et popularisée par Meyn-Tweedie est largement utilisée pour ce problème. Dans la plupart des résultats, le taux de convergence n’est pas explicite, et certains d’entre eux seront brièvement présentés. De plus, la fonction de Lyapunov est cruciale dans l’approche Meyn-Tweedie, et elle est aussi liée à certaines inégalités fonctionnelles (par exemple, inégalité de Poincaré). Cette relation entre fonction de Lyapounov et inégalités fonctionnelles sera donnée avec les résultats au sens L2. En outre, pour l’exemple de l’équation cinétique de Fokker-Planck, un résultat de convergence exponentielle explicite de la solution sera introduite à la manière de Villani : l’hypocoercivité. Ces contenus sont les fondements de mon travail, et mon but est d’étudier la décroissance sous-exponentielle. La deuxième partie, fait l’objet d’un article écrit en coopération avec d’autres sur les taux de convergence sous-exponentielle explicites des processus de Markov à temps continu. Comme nous le savons, les résultats sur les taux de convergence explicites ont été donnés pour le cas exponentiel. Nous les étendons au cas sous-exponentielle par l’approche Meyn-Tweedie. La clé de la preuve est l’estimation du temps de passage dans un ensemble ”petite”, obtenue par Douc, Fort et Guillin, mais pour laquelle nous donnons une preuve plus simple. Nous utilisons aussi la construction du couplage et donnons une ergodicité sous exponentielle explicite. Enfin, nous donnons quelques applications numériques. Dans la dernière partie, mon second article traite de l’équation cinétique de Fokker-Planck. Je prolonge l’hypocoercivité à l’hypocoercivité faible qui correspond à inégalité de Poincaré faible. Grâce à cette extension, on peut obtenir le taux de convergence explicite de la solution, dans des cas sous-exponentiels. La convergence est au sens H1 et au sens L2. A la fin de ce document, j’étudie le cas de l’entropie relative comme Villani, et j’obtiens la convergence au sens de l’entropie. Enfin, Je donne deux exemples pour les potentiels qui impliquent l’inégalité de Poincaré faible ou l’inégalité de Sobolev logarithmique faible pour la mesure invariante. / My Ph.D dissertation mainly focuses on long time behavior of Markov processes, functional inequalities and related techniques. More specifically, I will present the computable sub-exponential convergence rate of the Markov process in two approaches : Meyn-Tweedie’s method and (weak) hypocoercivity. The paper consists of three parts. In the first part, I will introduce some important results and related knowledge. Firstly, overviews of my research field are given. Exponential (or subexponential) convergence of Markov chains and (continuous time) Markov processes is a hot issue in probability. The traditional method - Meyn-Tweedie’s approach is widely applied for this problem. Most of the results about convergence rate is not explicit, and some of them will be introduced briefly. In addition,Lyapunov function is crucial in Meyn-Tweendie’s aproach, and it is also related to some functional inequalities (for example, Poincar´e inequality). The relationship of them will be given with results in L2 sense. Furthermore, as a example of kinetic Fokker-Planck equation, a computable result of exponential convergence of the solution of it will be introduced in Villani’ way - hypocoercivity. These contents are foundations of my work, and my destination is to study the sub-exponential decay. In the second part, it is my article cooperated with others about subexponential convergence rate of continuous time Markov processes. As we all know, the explicit results of convergence rate is about the exponential case. We extend them to sub-exponential case in Meyn-Tweedie’s approach. The key of the proof is the estimation of the hitting time to small set which was got by Douc, Fort and Guillin, for which we also propose an alternative simpler proof. We also use coupling construction as others and give a quantitative sub-exponential ergodicity. At last, we give some calculations for examples. In the last part, my second article deal with the kinetic Fokker-Planck equation. I extend the hypocoercivity to weak hypocoercivity which correspond to weak Poincar´e inequality. Through the extension, one can get the computable rate of convergence of the solution, which is also sub-exponential case. The convergence is in H1 sense and in L2 sense. In the end of this paper, I study the relative entropy case as C.Villani, and get convergence in entropy. Finally, I give two examples for potentials that implies weak Poincar´e inequality or weak logarithmic Sobolve inequality for invarient measure.