Cette thèse est consacrée à l’exploration des enjeux philosophiques, religieux et littéraires d’un motif qui apparaît dans trois œuvres majeures de la latinité : celui de la révélation finale. Le dernier livre du De republica de Cicéron avec le « Songe de Scipion », le dernier livre des Métamorphoses d’Ovide avec le discours de Pythagore, et le dernier livre XI des Métamorphoses d’Apulée avec le dévoilement d’Isis, présentent en effet tous trois la particularité de mettre en scène une expérience de révélation. Il s’agit alors non seulement d’étudier la façon dont Cicéron, Ovide et Apulée élaborent une écriture de la révélation à partir de matériaux philosophiques et religieux qu’on aura préalablement replacés dans leur contexte historique ainsi que dans l’économie générale des œuvres et des interrogations respectives de ces trois auteurs, mais aussi de déterminer les effets produits par la place de ces révélations en ce lieu traditionnel de polarisation du sens qu’est la fin d’une œuvre : quel dialogue ces trois révélations entretiennent-elles avec la structure qui les précède ? dans quelle mesure ces révélations finales dans les textes sont-elles aussi des révélations finales des textes eux-mêmes, auxquels elles apportent une élucidation rétrospective ? et quelle Weltanschauung, quelle conception de la quête de la vérité cette dynamique traduit-elle dans chacun des cas ? Cette enquête sera précédée d’un chapitre préliminaire où l’on établira les outils conceptuels nécessaires à l’étude de ce motif et où l’on analysera les antécédents grecs de celui-ci, à savoir le deus ex machina dans la tragédie et le mythe eschatologique final chez Platon. / This doctoral thesis is devoted to the investigation of the philosophical, religious, and formal aspects of a motif which appears in three major works of Latin literature : that of final revelation. The last book of Cicero’s De republica with its “Dream of Scipio”, the last book of Ovid’s Metamorphoses with its discourse of Pythagoras, and the last book of Apuleius’ Metamorphoses with its unveiling of Isis, all share, indeed, the specificity of representing a revelatory experience. The aim of this research is then : on the one hand, to study the way Cicero, Ovid and Apuleius dream up a pattern of revelation on the basis of philosophical and religious material that we shall first replace in each case in its historical context and in the frame of the respective literary output and pervasive themes of thinking of our three writers ; on the other hand, to determine the effects given off by the setting of these revelations at the place that traditionally focuses meaning, the ending of the work : what kind of dialogue is set up between these three revelations and the preceding structure ? to what extent do those final revelations in the texts function at the same time as final revelations of the texts themselves, on which they cast retrospective lighting ? and what worldview, what conception of the search for truth does this dynamics convey in the three cases ? This investigation will be preceded by a preliminary chapter where we shall establish the conceptual tools required for the study of this motif and analyze its Greek antecedents, i.e. the deus ex machina in tragedy and the final eschatological myth in Plato.

Ces travaux portent sur la modélisation et la simulation numérique de la dispersion de matériaux nocifs (pulvérisations liquides ou gazeuses) en milieu urbain ou naturel (attentat ou explosion accidentelle survenant en zone peuplée, fuites de produits toxiques gazeux ou liquides, éclatement de réservoir..). Afin de prédire ces risques un outils de simulation tridimensionnel a été développé. Celui-ci est basé sur un modèle de milieu hétérogène afin de traiter des phénomènes dont la durée et les distances associées peuvent être très grandes. La topographie des milieux étudiées est prise en compte grâce à des données numériques d'´elévation ainsi que les conditions météo permettant l'utilisation de profils de température et de vent complexes. Les transferts de chaleur et de masse sont considérés, notamment au niveau des obstacles. Un schéma numérique d'ordre élevé en temps et en espace est utilisé pour calculer les concentrations massiques de polluants. Par ailleurs, un modèle d'écoulement gaz-particule a été développé et implémenté dans le code de calcul. L'instabilité d'une couche de fluide soumise à un important gradient de pression est également étudiée, ceci afin de mieux comprendre et de caractériser les conditions initiales à utiliser pour ce type d'écoulement, impliquant des couches de particules. / This work deals with the modeling and the numerical simulation of the dispersion of toxic cloud of dropplets or gas in uneven geometry such as urban environment, industrial plants and hilly environment. Examples of phenomena under study are the dispersion of chemical products from damaged vessels, gas diffusion in an urban environment under explosion conditions, shock wave propagation in urban environment etc. A 3D simulation code has been developed in this aim. To simplify the consideration of complex geometries, a heterogeneous discrete formulation has been developed. When dealing with large scale domains, such as hilly natural environment, the topography is reconstructed with the help of numerical elevation data. Meteorological conditions are also considered, concerning temperature and wind velocity profiles. Heat and mass transfers on subscale objects, such as buildings are studied. A high order numerical scheme in space and time is used to compute mass concentration of pollutant. A two-phase model for dilute gas-particles flow has been developed and implemented in the 3D simulation code. The instability of a fluid layer appearing under high pressure gradient is also studied. This analysis allows us a better understanding of initial conditions for similar problems involving particles layer.

L'objectif de ma thèse est d'étudier mathématiquement un indicateur de rupture de volatilité très utilisé par les praticiens en salle de marché. L'indicateur bandes de Bollinger appartient à la famille des méthodes dites d'analyse technique et donc repose exclusivement sur l'historique récente du cours considéré et un principe déduit des observations passées des marchés, indépendamment de tout modèle mathématique. Mon travail consiste à étudier les performances de cet indicateur dans un univers qui serait gouverné par des équations différentielles stochastiques (Black -Scholes) dont le coefficient de diffusion change sa valeur à un temps aléatoire inconnu et inobservable, pour un praticien désirant maximiser une fonction objectif (par exemple, une certaine utilité espérée de la valeur du portefeuille à une certaine maturité). Dans le cadre du modèle, l'indicateur de Bollinger peut s'interpréter comme un estimateur de l'instant de la prochaine rupture. On montre dans le cas des petites volatilités, que le comportement de la densité de l'indicateur dépend de la volatilité, ce qui permet pour un ratio de volatilité assez grand, de détecter via l'estimation de la distribution de l'indicateur dans quel régime de volatilité on se situe. Aussi, dans le cas des grandes volatilités, on montre par une approche via la transformée de Laplace, que le comportement asymptotique des queues de distribution de l'indicateur dépend de la volatilité. Ce qui permet de détecter le changement des grandes volatilités. Ensuite, on s'intéresse à une étude comparative entre l'indicateur de Bollinger et l'estimateur classique de la variation quadratique pour la détection de changement de la volatilité. Enfin, on étudie la gestion optimale de portefeuille qui est décrite par un problème stochastique non standard en ce sens que les contrôles admissibles sont contraints à être des fonctionnelles des prix observés. On résout ce problème de contrôle en s'inspirant de travaux de Pham and Jiao pour décomposer le problème initial d'allocation de portefeuille en un problème de gestion après la rupture et un problème avant la rupture, et chacun de ces problèmes est résolu par la méthode de la programmation dynamique . Ainsi, un théorème de verification est prouvé pour ce problème de contrôle stochastique. / The aim of my thesis is to study mathematically an indicator widely used by the practitioners in the trading market, and designed to detect changes in the volatility term . The Bollinger Bands indicator belongs to the family of methods known as technical analysis which consist in looking t the past price movement in order to predict its future price movements independently of any mathematical model. We study the performance of this indicator in a universe that is governed by a stochastic differential equations (Black-Scholes) such that the volatility changes at an unknown and unobservable random time, for a practitioner seeking to maximize an objective function (for instance, the expected utility of the wealth at a certain maturity). Within the framework of the model, Bollinger indicator can be interpreted as an estimator of the time at which the volatility changes its value. We show that in the case of small volatilities, the density behavior of the indicator depends on the value of the volatility, which allows that for large ratio of volatility, to detect via the distribution estimation in which regime of volatility we are. Also , for the case of large volatilities, we show by an approach via the Laplace transform that the asymptotic tails behavior of the indictor depends on the volatility value. This allows to detect a change for large volatilities. Next, we compare two indicators designed to detect a volatility change: the Bollinger bands and the quadratic variation indicators. Finally, we study the optimal portfolio allocation which is described by a non-standard stochastic problem in view of that the admissible controls need to be adapted to the filtration generated by the prices. We resolve this control problem by an approach used by Pham and Jiao to separate the initial allocation problem into an allocation problem after the rupture and an problem before the rupture, and each one of these problems is resolved by the dynamic programming method. Also, a verification theorem is proved for this stochastic control problem.

L’amélioration des images sismiques peut aider à mieux contraindre l’exploration deshydrocarbures. Les ondes élastiques qui se propagent dans la Terre peuvent être classifiéescomme ondes de volume et ondes de surface. Si ces dernières sont les plus énergétiques,seules les ondes de volume sont couramment considérées comme des signaux utiles.Cependant, les ondes de surface sont utiles pour caractériser la proche sub-surface.Classiquement, les ondes de surface sont analysées dans des contextes de propriétésélastiques localement 1D.Nous proposons une modification de l’inversion des formes d’onde classique pourreconstruire des profils de propriétés 2D (la windowed-Amplitude Waveform Inversion, w-AWI). La w-AWI est spécialement robuste en ce qui concerne le choix du modèle initial.Nous appliquons la w-AWI aux données synthétiques ainsi qu’aux données réelles, montrantque cette approche permet de récupérer des propriétés 2D. / High-resolution seismic imaging is essential to improve results of hydrocarbon exploration.Elastic waves propagate in the Earth as body and surface waves, the latter being the mostenergetic ones. Body waves are preferred for exploration seismic imaging while surfacewaves are usually considered to be noise. However, it has been recognised that the nearsurface can be characterised by analysing surface waves and that such result may improvethe outcome of body-wave processing. Currently, surface waves analysis leads to retrievelocal 1D property profiles.We propose a waveform-based inversion procedure to derive 2D velocity models fromsurface waves. This method consists of a misfit functional modification of classical FullWaveform Inversion and we call it windowed-Amplitude Waveform Inversion (w-AWI). Weshow that w-AWI is robust regarding the choice of initial velocity model. We apply w-AWI tosynthetic and real data obtaining encouraging near-surface imaging results

Ce travail de thèse est constitué de l'étude de deux problèmes inverses associés à des équations paraboliques à coefficients périodiques. Dans la première partie, on a considéré une équation parabolique à coefficients et condition initiale périodiques. Notre travail a consisté à aborder le cas de coefficient à régularité faible et à minimiser les contraintes d'observations requises pour établir notre résultat de reconstruction du potentiel. On a commencé par établir un résultat d'existence et d'unicité de la solution dans un espace d'énergie adéquat. Ensuite, on a énoncé un principe du maximum adapté aux hypothèses du problème étudié et on a travaillé avec des coefficients mesurables et bornés. Enfin, on a reconstruit le potentiel en établissant une inégalité de Carleman. Le résultat d'identification a été obtenu via une inégalité de stabilité de type Lipschitz. Dans le second travail, on s'est intéressé à la détermination d'un coefficient périodique en espace du terme de réaction dans une équation de réaction-diffusion définie dans l'espace entier $mathbb{R}$. On établit un résultat d'unicité en utilisant un nouveau type d'observations. La nature du problème étudié, posé dans l'espace $mathbb{R}$, nous a permis d'utiliser la notion de vitesse asymptotique de propagation. On a prouvé l'existence de cette vitesse et on l'a caractérisé. On a surdéterminé le problème inverse en choisissant une famille de conditions initiales à décroi-ssance exponentielle. Notre principal résultat est que ce coefficient est déterminé de façon unique, à une symétrie près, par l'observation d'un continuum de vitesses asymptotiques de propagation. / This thesis consists in the study of two problems associated to inverse para-bolic equations with periodic coefficients. We are interested in identifying one coefficient by using two different methods. In the first part, we consider a parabolic equation with periodic coefficients and periodic initial condition. Our work consists to consider the case of coefficient with weak regularity and to minimize the constraints of observations which are required to establish our reconstruction result. We establish a result of existence and uniqueness of the solution in adequate energy space. Then we prove a maximum principle adapted to the hypothesis of the problem studied and we work with measurable and bounded coefficients. Finally, we reconstruct the potential by establishing a Carleman estimate. The identification result was achieved via an inequality of stability. In the second work, we are interested to determine a periodic coefficient of the reaction term defined in the whole space $mathbb{R}$. We establish a uniqueness result by using a new type of observations. The nature of the studied problem allowed us to use the notion of asymptotic speed of propagation. We prove the existence of this speed and we give its characterization. We overdetermin the inverse problem by choosing a family of initial conditions exponentially decaying. Our main result is that the coefficient is uniquely determined up to a symmetry, by the observation of a continuum of asymptotic speed of propagation.

L'imagerie terrestre urbaine se développe grâce à la diffusion grand public de visualisateurs immersifs au niveau de la rue. L'objectif de cette thèse est de proposer une méthode physique d'estimation de la réflectance des matériaux d'une scène urbaine à partir des images d'un système mobile d'acquisition, comme par exemple le véhicule d'acquisition Stereopolis, développé au laboratoire MATIS de l'Institut National de l'Information Géographique et Forestière (IGN).Pour cela, on considère que l'on dispose d'un modèle 3D de la scène, segmenté en zones de réflectance homogène (décrite par un modèle paramétrique), d'un jeu d'images géoréférencées corrigées radiométriquement. On propose une méthode d'estimation des réflectances par minimisation de l'écart entre les images réelles acquises par le capteur et des images simulées depuis le même point de vue. Après une modélisation phénoménologique des différentes composantes de la luminance arrivant au niveau d'un capteur imageur dans le domaine visible, une méthode utilisant le lancer de rayons sert à simuler cette luminance. Cela constitue ici le problème direct. Cet outil de simulation nécessite la connaissance de l'illumination de la scène, i.e. la répartition et la puissance des sources de lumière. Pour une acquisition avec Stereopolis, on ne réalise généralement pas de mesures atmosphériques et/ou radiométriques permettant de déterminer l'illumination avec des codes de transfert radiatif ; c'est pourquoi on propose une méthode d'estimation de la luminance du ciel à partir des images, en utilisant les pixels qui voient le ciel. L'éclairement solaire direct, non accessible directement dans les images, est estimé par une méthode ombre-soleil grâce à une plaque de référence placée sur le toit du véhicule. L'algorithme d'inversion du système se fait par minimisation d'une fonction coût constituée par la différence pixel à pixel entre les images simulées avec certains paramètres de réflectance et les images réelles. Cela nécessite de nombreuses simulations par lancer de rayons, car l'algorithme est itératif en raison des réflexions multiples entre les objets qui doivent être calculées avec des paramètres de réflectances initiaux. Afin d'éviter ces très coûteux lancers de rayons, on propose un algorithme de lancer de rayons formel qui stocke la luminance simulée comme des fonctions des paramètres de réflectances au lieu d'une valeur numérique. Cela permet de mettre à jour les images de luminance simulées par simple évaluation de ces fonctions avec le jeu de paramètres courant. La minimisation elle-même est effectuée par la méthode du gradient conjugué. Des résultats sur des scènes synthétiques permettent de faire une première validation de la méthode. Cependant, l'application sur un jeu d'images issues de Stereopolis pose plusieurs difficultés, notamment liées l'étalonnage radiométrique et à la segmentation du modèle 3D utilisé en entrée / Urban terrestrial imagery is widely used through online viewers of street-level images. The MATIS of the French National Geographical and Forester Data Institute (IGN) has developed its own mobile-mapping vehicle, Stereopolis, dedicated to research purposes. In this thesis, we develop a physically-based method to retrieve the reflectance of urban materials from a set of images shot by a mobile-mapping vehicle. This method uses a 3D model of the scene (segmented in areas of homogeneous reflectance, modeled by a parametric formula) and a set of georeferenced and radiometrically corrected images. We present a method for the estimation of the reflectances of the materials by minimizing the difference between real acquired images and simulated images from the same point of view. A modelisation of the physical phenomena leading to the formation of optical images is presented. A code using ray-tracing algortihm is used to compute the radiance at the sensor level. This is the direct problem for our estimation of reflectances. The illumination of the scene must be an input of this simulation tool; now, in a typical urban mobile-mapping acquisition neither atmospheric nor radiometric measurements are performed, that could be used to determine the illumination with a radiative transfer code. We propose an estimation of the illumination using directly the sky pixels in the acquired images. The direct solar irradiance cannot be estimated from the images because of overexposure, but we can use a shadows casted on a dedicated reference plate placed on top of the vehicle. The reflectances estimation is performed by minimizing a cost function; this cost function is the pixel-wise difference between simulated images (with current reflectances parameters) and acquired images. This leads to numerous ray-tracing simulations as the algorithm is iterative due to interreflections between the objects of the scene, that are computed using initial parameters. In order to prevent these costly ray-tracing simulations, we propose a symbolic ray-tracing algorithm that computes the radiance as a symbolic function of the reflectances parameters instead of a numerical value. Then, each iteration of the minimization algorithm is only an evaluation of a symbolic function. Results are shown on synthetic scenes to perform a first validation the estimation method. Using this method on real Stereopolis images remains difficult, mainly due to the radiometric calibration of cameras and the segmentation on the 3D model

La détection d'impacts est une préoccupation majeure pour suivre l'intégrité d'une structure évoluant dans un environnement comportant des risques de collision. Ce travail consiste à développer une méthode d’identification d'impact applicable à une structure composite de grandes dimensions qui soit à la fois robuste, peu coûteuse en temps de calcul, et qui nécessite peu de capteurs. Dans un premier temps, la démarche a été de décrire l'image du point d'impact dans la réponse vibratoire par un vecteur de participations modales. L'idée a ensuite été d'introduire l'hypothèse d'existence d'une famille de modes discriminants permettant d'établir le lien bijectif entre les participations modales et le point d'impact. Une procédure d'estimation de ces participations modales à partir d'un unique point de mesure a été proposée. Dans un second temps, la démarche a été d'étendre la procédure à l'identification d'une loi décrivant les paramètres principaux d'une force d'impact. Des conditions portant sur les paramètres de mesure et les propriétés modales de la structure sont déterminées pour garantir la précision et la robustesse de l'identification. Ce travail a permis de développer une approche en rupture par rapport à l'état de l'art, en ce sens qu'elle ne nécessite, en théorie, qu'un seul point de mesures vibratoires pour identifier un impact. Des essais expérimentaux sur un pavillon d'A350 ont permis de valider cette approche sur une structure composite de grandes dimensions. / Many industrial structures operate in an environment with a high risk of collision. The detection of impacts and the assessment of their severity is a major preoccupation in Structural Health Monitoring. This work deals with the development of an impact identification technique that is applicable to a large composite structure, numerically robust and time efficient, and that requires a low number of sensors. The first step was to describe the image of the impact point in the vibration response as a modal participation vector. The idea was to introduce the existence of a discriminating modes familly to make a bijective link between the modal participations and the impact point. A least-squares procedure is developed to estimate those modal participations by measuring the vibration response of a single point on the structure. The second step was to extend the procedure to the identification of a parametric law representing the impact load history. Some conditions on measurement parameters and modal properties of the structure are identified to guarantee both the accuracy and the robustness of the procedure. The approach developed in this work is new regarding the state of the art: only one measurement point is considered for identifying an impact event. An experimental validation on an A350 crown panel indicates that the methodology is valid for large composite structures.

Dans la thèse présente, nous nous intéressons à l'étude de la relative hyperbolicité des produits semi-direct des produits libres, ainsi que le problème de conjugaison pour certains automorphismes de ces produits libres.Plus précisement, pour un produit libre $$G=G_1astdotsast G_past F_k$$ un automorphisme $phi$ est intitulé atoroidal s'il ne fixe pas (ni aucune de ses puissances) la classe de conjugaison d'un élément hyperbolique de $G$. Cet automorphisme est appelé completement irréductible si le système de facteurs libres est le plus grand qui est fixé par toutes les puissances de cet automorphisme. Il est appelé toral si pour tous les $i$, il existe $g_iin G$ tel que ${rm ad}_{g_i}circ phi|_{G_i}$ est identité sur le facteur libre $G_i$. Nous disons qu'il a la condition centrale si pour chaque $i$, il existe $g_iin G$ conjugue $phi(G_i)$ à $G_i$, et s'il existe un élément non trivial de $G_irtimes_{{rm ad}_{g_i} circ phi|_{G_i}} mathbb{Z}$ qui est central dans $G_irtimes_{{rm ad}_{g_i} circ phi|_{G_i}} mathbb{Z}$.Nous prouvons, dans le Théorème 4.28, que si $phi$ est atoroidal et completement irréductible, et si le produit libre est non-elementaire ($kgeq 2$ ou $ p+k geq 3$), le groupe $Grtimes_phi mathbb{Z}$ est relativement hyperbolique (relativement a des suspensions de chaque $G_i$). Après, dans le Théorème 6.10, nous prouvons le même résultat si $phi$ est atoroidal avec la condition centrale. Nous prouvons aussi dans le Théorème 7.21 que si tous les $G_i$ sont abelien, le problème de conjugaison est solvable pour les automorphismes atoroidaux, toraux. Ces sont des analogues du résultat de Brinkmann [7] (celui qui a donné le résultat d'hyperbolicité pour les groupes libres), et du résultat de Dahmani [12] (celui qui a résolu le problème de conjugaison des automorphismes hyperboliques). / In this thesis, we are interested in the study of the relative hyperbolicity of the suspensions of free products, as well as the conjugacy problem of certain automorphisms of free products.To be more precise, given a free product $$G=G_1astdotsast G_past F_k$$ an automorphism $phi$ is said atoroidal if no power fixes the conjugacy class of an hyperbolic element. It is called fully irreducible if the given free factor system $[G_1],dots,[G_p]$ is the largest one that is fixed by every power of the automorphism. It is said toral if for all $i$, there exists $g_iin G$ such that ${rm ad}_{g_i}circ phi|_{G_i}$ is the identity on the free factor $G_i$. It is said to have central condition if for each $i$, there exists $g_iin G$ conjugating $phi(G_i)$ to $G_i$, and if there exists a non-trivial element of $G_irtimes_{{rm ad}_{g_i} circ phi|_{G_i}} mathbb{Z}$ that is central in $G_irtimes_{{rm ad}_{g_i} circ phi|_{G_i}} mathbb{Z}$.We prove, in Theorem 4.28, that if $phi$ is atoroidal and fully irreducible, and if the free product is non-elementary ($kgeq 2$ or $ p+k geq 3$), the group $Grtimes_phi mathbb{Z}$ is relatively hyperbolic (relative to the mapping torus of each $G_i$). Then in Theorem 6.10 we prove the same result holds if $phi$ is atoroidal with central condition. We also prove in Theorem 7.21 that if all $G_i$ are abelian, the conjugacy problem is solvable for toral atoroidal automorphisms. These are analogue of the result of Brinkmann [7] (which gave the hyperbolicity result for free groups) and the result of Dahmani [12] (which solved the conjugacy problem of hyperbolic automorphisms).

Ce travail traite des méthodes de calcul numérique pour les simulations hydrodynamiques appliquées principalement sur des produits développés par ANDRITZ HYDRO. Il s’agit ici de mettre en place des schémas d’ordre élevé pour des simulations CFD en utilisant le code de calcul ASPHODEL développé et utilisé par ANDRITZ HYDRO. Les principales motivations sont l’augmentation de la fiabilité des résultats de calculs numériques avec un coût de calcul raisonnable. Cette fiabilité s’exprime à travers l’augmentation de la précision et de la robustesse des schémas numériques. Le code de calcul ASPHODEL est basé sur la méthode sans maillage SPH-ALE. Mélange entre les volumes finis et la méthode SPH (Smoothed Particle Hydrodynamics), la méthode SPH-ALE emploie un ensemble de points appelés particules servant à la discrétisation du domaine fluide. Elle permet en particulier de par son caractère sans maillage, un suivi des surfaces libres sans effort de calcul supplémentaire. Cet aspect est véritablement attrayant pour bon nombre d’applications industrielles notamment la simulation des écoulements à surface libre se produisant dans une turbine Pelton, mais également le remplissage d’une turbine Francis. Cependant, le bémol à cette méthode est son manque de précision spatiale. En effet les points de calcul étant mobiles, les opérateurs spatiaux doivent être en mesure de conserver leur précision et leur robustesse au cours du temps. La qualité des résultats en est du coup impactée, en particulier le champ de pression souvent excessivement bruité. La montée en ordre et l’amélioration de la consistance des opérateurs pour un vaste panel de configurations géométriques sont donc les enjeux de ce travail. En utilisant des outils inspirés par les volumes finis non-structurés, il est possible d’améliorer les opérateurs spatiaux. En effet, la montée en ordre ou p-raffinement peut notamment se faire avec des reconstructions d’ordres élevés pour évaluer les états aux interfaces des problèmes de Riemann. La sommation des flux numériques résolus par un solveur de Riemann est ensuite retravaillée pour obtenir un schéma numérique d’ordre global cohérent. Le même soucis de cohérence avec les schémas en temps doit d’ailleurs être pensé. Le gain de précision apporté par les schémas numériques d’ordre élevé est comparé avec un raffinement spatial, c’est à dire une augmentation du nombre des particules de taille plus petite, aussi appelé h-raffinement. La méthode SPH-ALE améliorée est ensuite testée sur des cas représentatifs des applications visées. En conclusion, les développements effectués dans cette étude ont été guidés par l’application en turbine Pelton principalement mais il va de soi qu’ils sont applicables à des écoulements sans surface libre dans les turbines Francis par exemple. Ce travail montre les possibilités d’une méthode sans maillage pour des cas d’écoulements complexes autour de géométrie tournantes. / This work deals with numerical methods for hydrodynamic testing applied mainly on products developed by ANDRITZ HYDRO. This is to put in place high order schemes for CFD simulations using the ASPHODEL calculation code developed and used by ANDRITZ HYDRO. The main reasons are the increased reliability of the results of numerical calculations with a reasonable computational cost. This reliability is expressed through increasing the accuracy and robustness of numerical schemes. The ASPHODEL computer code is based on the meshfree method SPH-ALE. Mix between finite volume method and SPH (Smoothed Particle Hydrodynamics), the SPH-ALE method uses a set of points called particles serving as the fluid domain discretization. It allows track free surfaces without additional computational effort. This is truly attractive for many industrial applications including the simulation of free surface flows occurring in a Pelton turbine, but also filling a Francis turbine. However, the downside of this method is its lack of spatial accuracy. Indeed calculation points are mobile, space operators must be able to keep their accuracy and robustness over time. The quality of results is impacted especially the pressure field is often excessively noisy. The rise in order and improving the consistency of the operators for a wide range of geometric configurations are the challenges of this work. Using tools inspired by the unstructured finite volume, it is possible to improve the spatial operators. Indeed, the increasing order or p-refinement particular can be done with reconstructions of high order to assess the conditions at the interfaces of Riemann problems. The summation of discret fluxes solved by Riemann solver is then reworked to obtain a coherent global order scheme. The same concern for consistency with temporal schemes should also be considered. The precision gain provided by numerical schemes of higher orders is compared with a spatial refinement ie an increase in the number of smaller particles ; also called h -refinement . Improved SPH -ALE method is then tested on representative cases of intended applications. In conclusion, the developments made in this study were guided in accordance mainly with the Pelton turbine but it goes without saying that they are applicable to non- free surface flows in Francis turbines for example. This work shows the possibilities of a free mesh method for cases of complex flow around rotating geometry.

Le travail de ma thèse est basé sur ces 4 points :i) Transport laplacien d'une cellule absorbante :Soit un certain espèce (cellule) de concentration C(x), qui diffuse dans un milieu homogène et isotrope à partir d'une lointaine source localisée sur la frontière fermée $partial Omega_{0}$ vers une interface compact semi-perméable $partial Omega$ (membrane de la "cellule") à laquelle elle disparaisse àun taux d'absorption donné : W>=0. La concentration C (transport laplacien avec un coefficient de diffusion D) satisfaite le problème (P1) (voir la thèse). On s'intéresse à résoudre le problème (P1) en dimension dim = 2; 3 et à calculer les courants local et total à travers les frontières des $partial Omega$ et $partial Omega_{0}$ qui seront utiles pour résoudre le problèmeinverse de localisation. Pour faciliter les calculs et les rendre explicites, on prend $partial Omega$ et $partial Omega_{0}$ avec des formes géométriquement régulières, précisément des boules, en distinguant les deux cas : $Omega$ et $Omega_{0}$ sont concentriques ou non-concentriques. Pour le cas non-concentriques , on utilise la technique de transformation conforme et le développement orthogonal en série de Fourier pour résoudre le problème (P1) en cas bidimensionnel. Tandis que en cas tridimensionnel, on résout le problème (P1) en utilisant le développement orthogonal suivant les fonctions sphériques harmoniques.ii) Problème inverse de localisationOn s'intéresse dans cette partie à résoudre le problème inverse de localisation associé au problème (P1) où les domaines $Omega$ et $Omega_{0}$ sont considérés avec des formes géométriques régulières (précisément des boules) . Ce problème consiste à trouver les conditions de Dirichlet-Neumann sur $partial Omega_{0}$ (courant local, courant total) suffisantes pour déterminer la position de la cellule $partial$ (par rapport à $Omega_{0}$), dont ces conditions sont disponibles par une suite des mesures expérimentales.iii) Problème invesre géomètrique :Dans cette partie on traite un autre type de problème inverse qui consiste à trouver la forme géométrique de la cellule en sachant les conditions de Dirichlet-Neumann au bord extérieur(partial Omega_{0}) qui sont mésurables par une suite d'expérience. Ce type du problème, on l'appelle le problème inverse géométrique. On résout ce problème en utilisant des techniques concernant les fonctions harmoniques et les transformations conformes.iv) Opérateur de Dirichlet-NeumannOn étudie l'opérateur de Dirichlet-Neumann relatif au problème (P1) dans les dimension deux et trois en distinguant les deux cas concentriques et non-concentriques. Ensuite, on montre que cet opérateur de Dirichlet-Neumann engendre certain semi-groupe qu'on l'appelle semi-groupe de Lax. Enfin, on construit ce semi-groupe de Lax associé à cet opérateur en cas tridimensionnel concentriques afin de vérifier que ce semi-groupe admet les mêmes propriétés que celui dans le cas général. / The outline of my thesisi) Let some "species" of concentration C(p), x 2 Rd, diuse stationary in the isotropic bulk from a (distant) source localised on the closed boundary $partial Omega_{0}$ towards a semipermeable compact interface $partial Omega$ of the cell $Omega in Omega_{0}$ where they disappear at a given rate $W >= 0$. Then the steady field of concentrations C satisfy the problem $(P1)$. (see the Thesis). We interest to solve (P1) in Twodimensional and Tridimensional cases and to calculate the local and total flux in order to solving the localisation inverse problem. In order to make easy the calculations, we take $Omega$ and $Omega_{0}$ with a regularly geometricals forms by distinguishing the two cases : Concentrics and non-concentrics case. For the non-cncentrics case, we use the conformal mapping technique for resolving the problem (P1) in the twodimensional case. whereas in the tridimensional case, we use the development according to the spherical harmonics functions.ii) Localisation inverse problemThe aim of the localisation inverse problem is to find the necessary Dirichlet-to-Neumann conditions in order to determine the position of thecell $Omega$, where these conditions are measurable.iii) Geometrical inverse problemOur main results concerns a formal solution of the geometrical inverse problem for the form of absorbing domains. We restrict this study to two dimensions and we study it by the conformal mapping technique and harmonic functions.iv) Dirichlet-to-Neumann operatorWe study the Dirichlet-to-Neumann operatot relative to problem (P1) in the twodimensional and tridimensionnal cases by distinguishing the two cases : Concentrics and non-concentrics case. We prove that the Dirichlet-to-Neumann operator generates some semi-group, we call it the Lax semi-group. Finally we construct this semi group and verify that this demi-group satisfies the generals properties of a operator.