My thesis addresses a systematic approach to stochastic models in population genetics; in particular, the Wright-Fisher models affected only by the random genetic drift. I used various mathematical methods such as Probability, PDE, and Geometry to answer an important question: \"How do genetic change factors (random genetic drift, selection, mutation, migration, random environment, etc.) affect the behavior of gene frequencies or genotype frequencies in generations?”. In a Hardy-Weinberg model, the Mendelian population model of a very large number of individuals without genetic change factors, the answer is simple by the Hardy-Weinberg principle: gene frequencies remain unchanged from generation to generation, and genotype frequencies from the second generation onward remain also unchanged from generation to generation. With directional genetic change factors (selection, mutation, migration), we will have a deterministic dynamics of gene frequencies, which has been studied rather in detail. With non-directional genetic change factors (random genetic drift, random environment), we will have a stochastic dynamics of gene frequencies, which has been studied with much more interests. A combination of these factors has also been considered. We consider a monoecious diploid population of fixed size N with n + 1 possible alleles at a given locus A, and assume that the evolution of population was only affected by the random genetic drift. The question is that what the behavior of the distribution of relative frequencies of alleles in time and its stochastic quantities are. When N is large enough, we can approximate this discrete Markov chain to a continuous Markov with the same characteristics. In 1931, Kolmogorov first introduced a nice relation between a continuous Markov process and diffusion equations. These equations called the (backward/forward) Kolmogorov equations which have been first applied in population genetics in 1945 by Wright. Note that these equations are singular parabolic equations (diffusion coefficients vanish on boundary). To solve them, we use generalized hypergeometric functions. To know more about what will happen after the first exit time, or more general, the behavior of whole process, in joint work with J. Hofrichter, we define the global solution by moment conditions; calculate the component solutions by boundary flux method and combinatorics method. One interesting property is that some statistical quantities of interest are solutions of a singular elliptic second order linear equation with discontinuous (or incomplete) boundary values. A lot of papers, textbooks have used this property to find those quantities. However, the uniqueness of these problems has not been proved. Littler, in his PhD thesis in 1975, took up the uniqueness problem but his proof, in my view, is not rigorous. In joint work with J. Hofrichter, we showed two different ways to prove the uniqueness rigorously. The first way is the approximation method. The second way is the blow-up method which is conducted by J. Hofrichter. By applying the Information Geometry, which was first introduced by Amari in 1985, we see that the local state space is an Einstein space, and also a dually flat manifold with the Fisher metric; the differential operator of the Kolmogorov equation is the affine Laplacian which can be represented in various coordinates and on various spaces. Dynamics on the whole state space explains some biological phenomena.

Ziel dieser Arbeit ist die Untersuchung des Einflusses der stochastischen Interpretation der Langevin Gleichung mit zustandsabhängigen Diffusionskoeffizienten auf den Propagator des zugehörigen stochastischen Prozesses bzw. dessen Mittelwerte. Dies dient dem besseren Verständnis und der Interpretation von Messdaten von Diffusion in inhomogenen Systemen und geht einher mit der Frage der Form der Diffusionsgleichung in solchen Systemen. Zur Vereinfachung der Fragestellung werden in dieser Arbeit nur Systeme untersucht die vollständig durch einen ortsabhängigen Diffusionskoeffizienten und Angabe der stochastischen Interpretation beschrieben werden können. Dazu wird zunächst für mehrere experimentell relevante eindimensionale Systeme der jeweilige allgemeine Propagator bestimmt, der für jede denkbare stochastische Interpretation gültig ist. Der analytisch bestimmte Propagator wird dann für zwei exemplarisch ausgewählte stochastische Interpretationen, hier für die Itô und Klimontovich-Hänggi Interpretation, gegenübergestellt und die Unterschiede identifiziert. Für Mittelwert und Varianz der Prozesse werden die drei wesentlichen stochastischen Interpretationen verglichen, also Itô, Stratonovich und Klimontovich-Hänggi Interpretation. Diese systematische Untersuchung von inhomogenen Diffusionsprozessen kann zukünftig helfen diese Art von, in genau einer stochastischen Interpretation, driftfreien Systemen einfacher zu identifizieren. Ein weiterer wesentlicher Teil der Arbeit erweitert die Frage auf mehrdimensionale inhomogene anisotrope Systeme. Dies wird z.B. bei der Untersuchung von Diffusion in Flüssigkristallen mit inhomogenem Direktorfeld relevant. Obwohl hier, im Gegensatz zu eindimensionalen Systemen, der Propagator nicht allgemein berechnet werden kann, wird dennoch der Einfluss der Inhomogenität auf Messgrößen, wie die mittlere quadratische Verschiebung oder die Verteilung der Diffusivitäten, bestimmt. Anhand eines Beispiels wird auch der Einfluss der stochastischen Interpretation auf diese Messgrößen demonstriert. / The aim of this thesis is to investigate the influence of the stochastic interpretation of the Langevin equation with state-dependent diffusion coefficient on the propagator of the related stochastic process, or its averages, respectively. This helps to obtain a deeper understanding and to interpret measurement data of diffusion in inhomogeneous systems and is accompanied with the question of the proper form of the diffusion equation in such systems. To simplify the question, in this thesis only systems are considered which can be fully described by a spatially dependent diffusion coefficient and a given stochastic interpretation. Therefore, for several experimentally relevant one-dimensional systems, the respective general propagator is determined, which is valid for any possible stochastic interpretation. Then, the propagator for two exemplary stochastic interpretations, here the Itô and Klimontovich-Hänggi interpretation, are compared and the differences are identified. For mean and variance of the processes three major interpretations are compared, namely the Itô, the Stratonovich and the Klimontovich-Hänggi interpretation. This systematic research on inhomogeneous diffusion process may help in future to identify these kind of, in exactly one stochastic interpretation, drift-free systems more easily. Another important part of this thesis extends this question to multidimensional inhomogeneous anisotropic systems. This is of high relevance, for instance, for the research of diffusion in liquid crystalline systems with an inhomogeneous director field. Although, in contrast to one-dimensional systems, the propagator may not be calculated generally, the influence of the inhomogeneity on measurement data like the mean squared displacement or the distribution of diffusivities is determined. Based on one example, also the influence of the stochastic interpretation on these quantities is demonstrated.

Ce travail de thèse a pour objet l'étude de modèles markoviens qui résultent de la prise en compte d'incertitudes dans des systèmes possédant une dynamique hybride : entrées bruitées, dynamique mal connue, ou évènements aléatoires par exemple. De tels modèles, parfois qualifiés de Systèmes Hybrides Stochastiques (SHS), sont utilisés principalement en automatique et en recherche opérationnelle.<br /><br />Nous introduisons dans la première partie du mémoire la notion de processus diffusif par morceaux, qui fournit un cadre théorique général qui unifie les différentes classes de modèles "hybrides" connues dans la littérature. Différents aspects de ces modèles sont alors envisagés, depuis leur construction mathématique (traitée grâce au théorème de renaissance pour les processus de Markov) jusqu'à l'étude de leur générateur étendu, en passant par le phénomène de Zénon.<br /><br />La deuxième partie du mémoire s'intéresse plus particulièrement à la question de la "propagation de l'incertitude", c'est-à-dire à la manière dont évolue la loi marginale de l'état au cours du temps. L'équation de Fokker-Planck-Kolmogorov (FPK) usuelle est généralisée à diverses classes de processus diffusifs par morceaux, en particulier grâce aux notions d'intensité moyenne de sauts et de courant de probabilité. Ces résultats sont illustrés par deux exemples de modèles multidimensionnels, pour lesquels une résolution numérique de l'équation de FPK généralisée a été effectuée grâce à une discrétisation en volumes finis. La comparaison avec des méthodes de type Monte-Carlo est également discutée à partir de ces deux exemples.

This study focusses on modelling important aspects of the evolution of pulsar wind nebulae using two different approaches. The first uses a hydrodynamic model to simulate the morphological evolution of a spherically-symmetric composite supernova remnant that is expanding into a homogeneous interstellar medium. In order to extend this model, a magnetic field is included in a kinematic fashion, implying that the reaction of the fluid on the magnetic field is taken into account, while neglecting any counter-reaction of the field on the fluid. This approach is valid provided that the ratio of electromagnetic to particle energy in the nebula is small, or equivalently, for a large plasma β environment. This model therefore allows one to not only calculate the evolution of the convection velocity but also, for example, the evolution of the average magnetic field. The second part of this study focusses on calculating the evolution of the energy spectra of the particles in the nebula using a number of particle evolution models. The first of these is a spatially independent temporal evolution model, similar to the models that can be found in the literature. While spatially independent models are useful, a large part of this study is devoted to developing spatially dependent models based on the Fokker-Planck transport equation. Two such models are developed, the first being a spherically-symmetric model that includes the processes of convection, diffusion, adiabatic losses, as well as the non-thermal energy loss processes of synchrotron radiation and inverse Compton scattering. As the magnetic field geometry can lead to the additional transport process of drift, the previous model is extended to an axisymmetric geometry, thereby allowing one to also include this process. / PhD (Space Physics), North-West University, Potchefstroom Campus, 2014

Bereits vor zehn Jahren gab es den ersten großen E-Book-Boom. Während damals noch digitale Buchinhalte für spezielle Lesegeräte im Mittelpunkt der Diskussion standen, dominieren inzwischen webbasierte Angebote den stetig wachsenden Markt. Für Bibliotheken eröffnen sich durch diese E-Books der nächsten Generation völlig neue Möglichkeiten in der Informationsversorgung ihrer Benutzer. Sowohl der globale Wissensaustausch, als auch der Zugriff auf Fachinformationen können mit ihrer Hilfe beschleunigt werden. Dies spielt insbesondere für Spezialbibliotheken, wie die der Max-Planck-Gesellschaft, eine große Rolle. Denn gerade von ihnen wird erwartet, dass sie die Mitarbeiter ihrer Institute stets mit den aktuellsten Fachinformationen versorgen. Aufgrund der Vielzahl an Lizenzierungs- und Angebotsformen auf dem E-Book-Markt ist hier ein Anbietervergleich unerlässlich. Dieser erfolgt im Rahmen dieser Diplomarbeit für die Bibliotheken der Max-Planck- Gesellschaft. Hierfür werden zunächst die aktuelle E-Book-Marktsituation und die Strukturen der Max-Planck-Gesellschaft beschrieben. Anschließend untersucht die Verfasserin die bisherigen Entwicklungen in der zentralen und lokalen E-Book-Erwerbung der Max-Planck-Gesellschaft. Zu diesem Zweck werden u.a. eine Befragung unter den Institutsbibliothekaren und mehrere Experteninterviews durchgeführt. Auf diese Weise können geeignete E-Book-Anbieter für die Max-Planck-Gesellschaft ermittelt und Kriterien gefunden werden, um diese schließlich miteinander vergleichen zu können. Im Ergebnis des Vergleichs gibt die Verfasserin Empfehlungen für die weitere E-Book-Erwerbung in der Max-Planck-Gesellschaft ab.

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.

The escape of trajectories is a ubiquitous phenomenon in open dynamical systems and stochastic processes. If escape occurs repetitively for a statistical ensemble of trajectories, the population of remaining trajectories often undergoes an exponential decay characterised by the so-called escape rate. Its inverse defines the lifetime of the decaying state, which represents an intrinsic property of the system. This paradigm is fundamental to nucleation theory and reaction-rate theory in chemistry, physics, and biology.<p><p>In many circumstances, escape is activated by the presence of noise, which may be of internal or external origin. This is the case for thermally activated escape over a potential energy barrier and, more generally, for noise-induced escape in continuous-time or discrete-time dynamics. <p><p>In the weak-noise limit, the escape rate is often observed to decrease exponentially with the inverse of the noise amplitude, a behaviour which is given by the van't Hoff-Arrhenius law of chemical kinetics. In particular, the two important quantities to determine in this case are the exponential dependence (the ``activation energy') and its prefactor.<p><p>The purpose of the present thesis is to develop an analytical method to determine these two quantities. We consider in particular one-dimensional continuous and discrete-time systems perturbed by Gaussian white noise and we focus on the escape from the basin of attraction of an attracting fixed point.<p><p>In both classes of systems, using path-integral methods, a formula is deduced for the noise-induced escape rate from the attracting fixed point across an unstable fixed point, which forms the boundary of the basin of attraction. The calculation starts from the trace formula for the eigenvalues of the operator ruling the time evolution of the probability density in noisy maps. The escape rate is determined by the loop formed by two heteroclinic orbits connecting back and forth the two fixed points in a two-dimensional auxiliary deterministic dynamical system. The escape rate is obtained, including the expression of the prefactor to van't Hoff-Arrhenius exponential factor./L'échappement des trajectoires est un phénomène omniprésent dans les systèmes dynamiques ouverts et les processus stochastiques. Si l'échappement se produit de façon répétitive pour un ensemble statistique de trajectoires, la population des trajectoires restantes subit souvent une décroissance exponentielle caractérisée par le taux d'échappement. L'inverse du taux d'échappement définit alors la durée de vie de l'état transitoire associé, ce qui représente une propriété intrinsèque du système. Ce paradigme est fondamental pour la théorie de la nucléation et, de manière générale, pour la théorie des taux de transitions en chimie, en physique et en biologie.<p><p>Dans de nombreux cas, l'échappement est induit par la présence de bruit, qui peut être d'origine interne ou externe. Ceci concerne en particulier l'échappement activé thermiquement à travers une barrière d'énergie potentielle, et plus généralement, l'échappement dû au bruit dans les systèmes dynamiques à temps continu ou à temps discret.<p><p>Dans la limite de faible bruit, on observe souvent une décroissance exponentielle du taux d'échappement en fonction de l'inverse de l'amplitude du bruit, un comportement qui est régi par la loi de van't Hoff-Arrhenius de la cinétique chimique. En particulier, les deux quantités importantes de cette loi sont le coefficient de la dépendance exponentielle (c'est-à-dire ``l'énergie d'activation') et son préfacteur.<p><p>L'objectif de cette thèse est de développer une théorie analytique pour déterminer ces deux quantités. La théorie que nous présentons concerne les systèmes unidimensionnels à temps continu ou discret perturbés par un bruit blanc gaussien et nous considérons le problème de l'échappement du bassin d'attraction d'un point fixe attractif. Pour s'échapper, les trajectoires du système bruité initialement contenues dans ce bassin d'attraction doivent alors traverser un point fixe instable qui forme la limite du bassin.<p><p>Dans le présent travail, et pour les deux types de systèmes, une formule est dérivée pour le taux d'échappement du point fixe attractif en utilisant des méthodes d'intégrales de chemin. Le calcul utilise la formule de trace pour les valeurs propres de l'opérateur gouvernant l'évolution temporelle de la densité de probabilité dans le système bruité. Le taux d'échappement est déterminé en considérant la boucle formée par deux orbites hétéroclines liant dans les deux sens les deux points fixes dans un système dynamique auxiliaire symplectique et bidimensionnel. On obtient alors le taux d'échappement, comprenant l'expression du préfacteur de l'exponentielle de la loi de van't Hoff-Arrhenius. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Dans le cadre des recherches de la qualité des ressources en eau, l’étude du processus de transfert de masse du sol vers la nappe phréatique constitue un élément primordial pour la compréhension de la pollution de cette dernière. En effet, les éléments polluants solubles à la surface (produits liés aux activités humaines tels engrais, pesticides...) peuvent transiter vers la nappe à travers le milieu poreux qu’est le sol. Ce scénario de transfert de pollution repose sur deux phénomènes : la pluie qui génère la masse d’eau à la surface et la dispersion de celle-ci à travers le milieu poreux. La dispersion de masse dans un milieu poreux naturel comme le sol forme un sujet de recherche vaste et difficile aussi bien au plan expérimental que théorique. Sa modélisation constitue une préoccupation du laboratoire EMMAH, en particulier dans le cadre du projet Sol Virtuel dans lequel un modèle de transfert (modèle PASTIS) a été développé. Le couplage de ce modèle de transfert avec en entrée un modèle décrivant la dynamique aléatoire de la pluie est un des objectifs de la présente thèse. Ce travail de thèse aborde cet objectif en s’appuyant d’une part sur des résultats d’observations expérimentaux et d’autre part sur de la modélisation inspirée par l’analyse des données d’observation. La première partie du travail est consacrée à l’élaboration d’un modèle stochastique de pluie. Le choix et la nature du modèle sont basés sur les caractéristiques obtenus à partir de l’analyse de données de hauteur de pluie recueillies sur 40 ans (1968-2008) sur le Centre de Recherche de l’INRA d’Avignon. Pour cela, la représentation cumulée des précipitations sera assimilée à une marche aléatoire dans laquelle les sauts et les temps d’attente entre les sauts sont respectivement les amplitudes et les durées aléatoires entre deux occurrences d’événements de pluie. Ainsi, la loi de probabilité des sauts (loi log-normale) et celle des temps d’attente entre les sauts (loi alpha-stable) sont obtenus en analysant les lois de probabilité des amplitudes et des occurrences des événements de pluie. Nous montrons alors que ce modèle de marche aléatoire tend vers un mouvement brownien géométrique subordonné en temps (quand les pas d’espace et de temps de la marche tendent simultanément vers zéro tout en gardant un rapport constant) dont la loi de densité de probabilité est régie par une équation de Fokker Planck fractionnaire (FFPE). Deux approches sont ensuite utilisées pour la mise en œuvre du modèle. La première approche est de type stochastique et repose sur le lien existant entre le processus stochastique issu de l’équation différentielle d’Itô et la FFPE. La deuxième approche utilise une résolution numérique directe par discrétisation de la FFPE. Conformément à l’objectif principal de la thèse, la seconde partie du travail est consacrée à l’analyse de la contribution de la pluie aux fluctuations de la nappe phréatique. Cette analyse est faite sur la base de deux relevés simultanées d’observations de hauteurs de pluie et de la nappe phréatique sur 14 mois (février 2005-mars 2006). Une étude statistique des liens entre les signaux de pluie et de fluctuations de la nappe est menée comme suit : Les données de variations de hauteur de nappe sont analysées et traitées pour isoler les fluctuations cohérentes avec les événements de pluie. Par ailleurs, afin de tenir compte de la dispersion de masse dans le sol, le transport de la masse d’eau pluviale dans le sol sera modélisé par un code de calcul de transfert (modèle PASTIS) auquel nous appliquons en entrée les données de hauteurs de pluie mesurées. Les résultats du modèle permettent entre autre d’estimer l’état hydrique du sol à une profondeur donnée (ici fixée à 1.6m). Une étude de la corrélation entre cet état hydrique et les fluctuations de la nappe sera ensuite effectuée en complément à celle décrite ci-dessus pour illustrer la possibilité de modéliser l’impact de la pluie sur les fluctuations de la nappe / Within the research quality of water resources, the study of the process of mass transfer from soil to groundwater is a key element for understanding the pollution of the latter. Indeed, soluble contaminants to the surface (related to human activities such fertilizers, pesticides products ...) can transit to the web through the porous medium that is the ground. This scenario transfer pollution based on two phenomena: the rain that generates the body of water to the dispersion and the surface thereof through the porous medium. The dispersion of mass in a natural porous medium such as soil forms a subject of extensive research and difficult both experimental and theoretical grounds. Its modeling is a concern EMMAH laboratory, particularly in the context of Virtual Sol project in which a transfer model (PASTIS model) was developed. The coupling of this transfer model with input a model describing the dynamics of random rain is one of the objectives of this thesis. This thesis addresses this goal by relying in part on the results of experimental observations and also on modeling inspired by the analysis of observational data. The first part of the work is devoted to the development of a stochastic model of rain. The choice and nature of the model are based on the features obtained from the analysis of data collected rainfall over 40 years (1968-2008) on the Research Centre INRA Avignon. For this, the cumulative rainfall representation will be treated as a random walk in which the jumps and waiting times between jumps are the amplitudes and durations between two random occurrences of rain events. Thus, the probability jumps (log-normal distribution) and that of waiting between jumps (Law alpha-stable) time is obtained by analyzing the laws of probability amplitudes and occurrences of rain events. We show that the random walk model tends towards a subordinate in time geometric Brownian motion (when space step and time step walking simultaneously tend to zero while maintaining a constant ratio), the law of probability density is governed by a Fokker Planck fractional (FFPE). Two approaches are then used to implement the model. The first approach is based on stochastic type and the relationship between the stochastic process derived from the differential equation of Itô and FFPE. The second approach uses a direct numerical solution by discretization of the FFPE. Accordance with the main objective of the thesis, the second part of the work is devoted to the analysis of the contribution of rain to fluctuations in groundwater. We approach this analysis on the basis of two simultaneous records of observations of rainfall amounts and groundwater over 14 months (February 2005-March 2006). A statistical study of the relationship between the signals of rain and fluctuating water will be conducted. Data sheet height variations are analyzed and processed to isolate coherent fluctuations with rain events. In addition, to take account of the mass dispersion in the soil, the mass transport of storm water in the soil layer is modeled by a calculation code transfer (PASTIS model) which we apply input data measured heights of rain. The model results allow between another estimate soil water status at a given depth (here set at 1.6m). A study of the correlation between the water status and fluctuating water will then be performed in addition to that described above to illustrate the ability to model the impact of rain on the water table fluctuations

Cette thèse porte sur les équations gyro-cinétiques et traite un développement rigoureux deslimites de l'équation de Vlasov avec différents opérateurs de collision dans un champ magnétiquefort, ainsi que du développement de méthodes numériques.On commence par une étude de l'opérateur de moyenne. L'opérateur de moyenne a été développé parM. Bostan dans le cadre général d'une équation pour laquelle une partie du transport estfortement pénalisée. Puis, on applique ces résultats généraux aux deux régimes limites que nousétudions : le régime du rayon de Larmor fini et le régime centre-guide.On s'intéresse au calcul précis et explicite de la moyenne de l'opérateur de Fokker-Planck-Landau. On se place pour cela dans le cas du régime du rayon de Larmor fini. Avant de réaliserles calculs sur l'opérateur de Fokker-Planck-Landau, qui contient des convolutions et des termesde diffusion, il semble raisonnable de calculer la moyenne de l'opérateur de relaxation deBoltzmann, dont l'expression est plus simple.On se place ensuite dans le cas du régime centre-guide et on présente un schéma numérique basésur une décomposition micro-macro de la fonction de distribution des particules qui provientd'un travail en collaboration avec N. Crouseilles et M. Lemou. On obtient un schéma uniformémentconsistant avec le modèle continu, pour tout ordre du champ magnétique. Des simulationsnumériques, basées sur cette approche, ont été réalisées à l'aide d'un code de calcul 2D quel'on a développé durant cette thèse.On présente ensuite un projet réalisé dans le cadre du Cemracs 2012, consacré à la modélisationdes écoulements sanguins dans le réseau veineux cérébral. / The main subject of this thesis is the gyro-kinetic equation. We present a rigourousdeveloppement of the Vlasov equation limits with different collision operator in a strongmagnetic field and numerical methods.We start with a study of the gyro-average operator. The average operator has been introduced byM. Bostan in the case of an equation where part of the transport is highly penalised. Then weapply our results at the two approximation we study : the finite Larmor radius approximation andthe guiding-center approximation.We first focus on the precise and explicit computation of the Fokker-Planck-Landau operatoraverage in the finite Larmor radius approximation. The Fokker-Planck-Landau operator containsconvolution and diffusion terms, it is then reasonable to first compute the average of theBoltzmann relaxation operator.We then focus on the guiding-center approximation and present a numerical scheme based on amicro-macro decomposition of the particles distribution fonction which comes from a joint workwith N. Crouseilles and M. Lemou. We obtain a scheme which is uniformly consistant with thecontinuous model for any order of the magnetic field. Numerical simulation based on thisapproach are presented.The last chapter of this thesis presents a project which was realised during the Cemracs 2012concerning the modelisation of blood flow in cerebral veins.

Prévoir et estimer le risque de canicule est un enjeu politique majeur. Évaluer la probabilité d'apparition des canicules et leurs sévérités serait possible en connaissant la température en temps continu. Cependant, les extrêmes journaliers (maxima et minima) sont parfois les seules données disponibles. Pour modéliser la dynamique des températures, il est courant d'utiliser un processus d'Ornstein-Uhlenbeck. Une estimation des paramètres de ce processus n'utilisant que les suprema journaliers observés est proposée. Cette nouvelle approche se base sur une minimisation des moindres carrés faisant intervenir la fonction de répartition du supremum. Les mesures de risque liées aux canicules sont ensuite obtenues numériquement. Pour exprimer explicitement ces mesures de risque, il peut être utile d'avoir la loi jointe du processus d'Ornstein-Uhlenbeck et de son supremum. L'étude se limite tout d'abord à la fonction de répartition / densité jointe du point final du processus et de son supremum. Cette probabilité admet une densité, solution de l'équation de Fokker-Planck, obtenue explicitement et utilisant les fonctions spéciales paraboliques cylindriques. La preuve de l'expression de la densité repose sur une décomposition sur une base hilbertienne de l'espace via une méthode spectrale. On étudie également le processus d'Ornstein-Uhlenbeck oscillant, dont le paramètre de drift est constant par morceaux selon le signe du processus. La transformée de Laplace du temps d'atteinte de ce processus est déterminée et la probabilité que le processus soit positif en un temps donné est calculée. / Forecasting and assessing the risk of heat waves is a crucial public policy stake. Evaluate the probability of heat waves and their severity can be possible by knowing the temperature in continuous time. However, daily extremes (maxima and minima) might be the only available data. The Ornstein-Uhlenbeck process is commonly used to model temperature dynamic. An estimation of the process parameters using only daily observed suprema of temperatures is proposed here. This new approach is based on a least square minimization using the cumulative distribution function of the supremum. Risk measures related to heat waves are then obtained numerically. In order to calculate explicitly those risk measures, it can be useful to have the joint law of the Ornstein-Uhlenbeck process and its supremum. The study is _rst limited to the joint density / distribution of the endpoint and supremum of the Ornstein-Uhlenbeck process. This probability admits a density, solution of the Fokker-Planck equation and explicitly obtained as an expansion involving parabolic cylinder functions. The proof of the density expression relies on a decomposition on a Hilbert basis of the space via a spectral method. We also study the oscillating Ornstein-Uhlenbeck process, which drift parameter is piecewise constant depending on the sign of the process. The Laplace transform of this process hitting time is determined and we also calculate the probability for the process to be positive on a fixed time.