Cette thèse porte sur une méthode cartésienne parallèle pour résoudre des problèmes elliptiques avec interfaces complexes et sur son application aux problèmes elliptiques en domaine irrégulier dans le cadre d'un modèle de croissance tumorale.La méthode est basée sur un schéma aux différences finies et sa précision est d’ordre deux sur tout le domaine. L'originalité de la méthode consiste en l'utilisation d'inconnues additionnelles situées sur l'interface et qui permettent d’exprimer les conditions de transmission à l'interface. La méthode est décrite et les détails sur la parallélisation, réalisé avec la bibliothèque PETSc, sont donnés. La méthode est validée et les résultats sont comparés avec ceux d'autres méthodes du même type disponibles dans la littérature. Une étude numérique de la méthode parallélisée est fournie.La méthode est appliquée aux problèmes elliptiques dans un domaine irrégulier apparaissant dans un modèle continue et tridimensionnel de croissance tumorale, le modèle à deux espèces du type Darcy . L'approche utilisée dans cette application est basée sur la pénalisation des conditions de transmission à l'interface, afin de imposer des conditions de Neumann homogènes sur le border d'un domaine irrégulier. Les simulations du modèle sont fournies et montrent la capacité de la méthode à imposer une bonne approximation de conditions au bord considérées. / This theses deals with a parallel Cartesian method to solve elliptic problems with complex interfaces and its application to elliptic irregular domain problems in the framework of a tumor growth model.This method is based on a finite differences scheme and is second order accurate in the whole domain. The originality of the method lies in the use of additional unknowns located on the interface, allowing to express the interface transmission conditions. The method is described and the details of its parallelization, performed with the PETSc library, are provided. Numerical validations of the method follow with comparisons to other related methods in literature. A numerical study of the parallelized method is also given.Then, the method is applied to solve elliptic irregular domain problems appearing in a three-dimensional continuous tumor growth model, the two-species Darcy model. The approach used in this application is based on the penalization of the interface transmission conditions, in order to impose homogeneous Neumann boundary conditions on the border of an irregular domain. The simulations of model are provided and they show the ability of the method to impose a good approximation of the considered boundary conditions. / Questa tesi introduce un metodo parallelo su griglia cartesiana per risolvere problemi ellittici con interfacce complesse e la sua applicazione ai problemi ellittici in dominio irregolare presenti in un modello di crescita tumorale.Il metodo è basato su uno schema alle differenze finite ed è accurato al secondo ordine su tutto il dominio di calcolo. L'originalità del metodo consiste nell'introduzione di nuove incognite sull'interfaccia, le quali permettono di esprimere le condizioni di trasmissione sull'interfaccia stessa. Il metodo viene descritto e i dettagli della sua parallelizzazione, realizzata con la libreria PETSc, sono forniti. Il metodo è validato e i risultati sono confrontati con quelli di metodi dello stesso tipo trovati in letteratura. Uno studio numerico del metodo parallelizzato è inoltre prodotto.Il metodo è applicato ai problemi ellittici in dominio irregolare che compaiono in un modello continuo e tridimensionale di crescita tumorale, il modello a due specie di tipo Darcy. L'approccio utilizzato è basato sulla penalizzazione delle condizioni di trasmissione sull'interfaccia, al fine di imporre condizioni di Neumann omogenee sul bordo di un dominio irregolare. Le simulazioni del modello sono presentate e mostrano la capacità del metodo di imporre una buona approssimazione delle condizioni al bordo considerate.

This dissertation addresses novel techniques in determining gyroless micro-satellite attitude and attitude rate. The main objective of this thesis is to explore the possibility of using commercially available low cost micro-light star sensor as a stand-alone sensor for micro-satellite attitude as well as attitude rate determination. The objective is achieved by developing accurate and computationally efficient algorithms for the realization of onboard operation of a low fidelity star sensor. All the algorithms developed here are tested with the measurement noise presented in the catalog of the sensor array STAR-1000. A novel accurate second order sliding mode observer (SOSMO) is designed for discrete time uncertain linear multi-output system. Our design procedure is effective for both matched and unmatched bounded uncertain ties and/or disturbances. The bound on uncertainties and/or disturbances is assumed to be unknown. This problem is addressed in this work using the second order multiple sliding modes approach. Second order sliding manifold and corresponding sliding condition for discrete time system is defined similar on the lines of continuous counterpart. Our design is not restricted to a particular class of uncertain (matched) discrete time system. Moreover, it can handle multiple outputs unlike single out-put systems. The observer design is achieved by driving the state observation error and its first order finite difference to the vicinity of the equilibrium point (0,0) in a finite steps and maintaining them in the neighborhood thereafter. The estimation synthesis is based on Quasi Sliding Mode (QSM) design. The problem of designing sliding mode observer for a linear system subjected to unknown inputs requires observer matching condition. This condition is needed to ensure that the state estimation error is a asymptotically stable and is independent of the unknown input during the sliding motion. In the absence of a matching condition, asymptotic stability of the reduced order error dynamics on the sliding surface is not guaranteed. However, unknown bounded inputs guarantee bounded error on state estimation. The QSM design guarantees an ultimate error bound by incorporating Boundary Layer (BL) in its design procedure. The observer achieves one order of magnitude improvement in estimation accuracy than the conventional sliding mode observer (SMO) design for an unknown input. The observer estimation errors, satisfying the given stability conditions, converge to an ultimate finite bound (with in the specified BL) of O(T2), where T Is the sampling period. A relation between sliding mode gain and boundary layer is established for the existence of second order discrete sliding motion. The robustness of the proposed observer with respect to measurement noise is also analyzed. The design algorithm is very simple to apply and is implemented for two examples with different classes of disturbances (matched and unmatched) to show the effectiveness of the design. Simulation results show the robustness with respect to the measurement noise for SOSMO. Second order sliding mode observer gain can be calculated off-line and the same gain can work for large band of disturbance as long as the disturbance acting on the continuous time system is bounded and smooth. The SOSMO is simpler to implement on board compared to the other traditional nonlinear filters like Pseudo-Linear-Kalman-filter(PLKF); Extended Kalman Filter(EKF). Moreover, SMO possesses an automatic adaptation property same as optimal state estimator(like Kalman filter) with respect to the intensity of the measurement noise. The SMO rejects the noisy measurements automatically, in response to the increased noise intensity. The dynamic performance of the observer on the sliding surface can be altered and no knowledge of noise statistics is required. It is shown that the SOSMO performs more accurately than the PLKF in application to micro-satellite angular rate estimation since PLKF is not an optimal filter. A new method for estimation of satellite angular rates through derivative approach is proposed. The method is based on optic flow of star image patterns formed on a star sensor. The satellite angular rates are derived directly from the 2D-coordinates of star images. Our algorithm is computationally efficient and requires less memory allocation compared to the existing vector derivative approaches, where there is also no need for star identification. The angular rates are computed using least square solution method, based on the measurement equation obtained by optic flow of star images. These estimates are then fed into discrete time second order sliding mode observer (SOSMO). The performance of angular rate estimation by SOSMO is compared with the discrete time First order SMO and PLKF. The SOSMO gives the best estimates as compared to the other two schemes in estimating micro-satellite angular rates in all three axes. The improvement in accuracy is one order of magnitude (around1.7984 x 10−5 rad/ sec,8.9987 x 10−6 rad/ sec and1.4222 x 10−5 rad/ sec in three body axes respectively) in terms of standard deviation in steady state estimation error. A new method and algorithm is presented to determine star camera parameters along with satellite attitude with high precision even if these parameters change during long on-orbit operation. Star camera parameters and attitude need to be determined independent of each other as they both can change. An efficient, closed form solution method is developed to estimate star camera parameters (like focal length, principal point offset), lens distortions (like radial distortion) and attitude. The method is based on a two step procedure. In the first step, all parameters (except lens distortion) are estimated using a distortion free camera model. In the second step, lens distortion coefficient is estimated by linear least squares (LS) method. Here the derived camera parameters in first step are used in the camera model that incorporates distortion. However, this method requires identification of observed stars with the catalogue stars. But, on-orbit star identification is difficult as it utilizes the values of camera calibrating parameters that can change in orbit(detector and optical element alignment get change in orbit due to solar pressure or sudden temperature change) from the ground calibrated value. This difficulty is overcome by employing a camera self-calibration technique which only requires four observed stars in three consecutive image frames. Star camera parameters along with lens (radial and decentering) distortion coefficients are determined by camera self calibration technique. Finally Kalman filter is used to refine the estimated data obtained from the LS based method to improve the level of accuracy. We consider the true values of camera parameters as (u0,v0) = (512.75,511.25) pixel, f = 50.5mm; The ground calibrated values of those parameters are (u0,v0) =( 512,512) pixel, f = 50mm; Worst case radial distortion coefficient affecting the star camera lens is considered to be k1 =5 x 10−3 .Our proposed method of attitude determination achieves accuracy of the order of magnitude around 6.2288 x 10−5 rad,3.3712 x 10−5 radand5.8205 x 10−5 rad in attitude angles φ,θ and ψ. Attitude estimation by existing methods in the literature diverges from the true value since they utilize the ground calibrated values of camera parameters instead of true values. To summarize, we developed a formal theory of discrete time Second Order Sliding Mode Observer for uncertain multi-output system. Our methods achieve the desired accuracy while estimating satellite attitude and attitude rate using low fidelity star sensor data. Our methods require lower on-board processing requirement and less memory allocation; thus are suitable for micro-satellite applications. Thus, the objective of using low fidelity star sensor as stand-alone sensor in micro-satellite application is achieved.

This thesis presents the construction, the analysis and the verication of compact residual discretizations for the solution of conservation laws on unstructured meshes. <p>The schemes considered belong to the class of residual distribution (RD) or fluctuation splitting (FS) schemes. <p>The methodology presented relies on three main elements: design of compact linear first-order stable schemes for linear hyperbolic PDEs, a positivity preserving procedure mapping stable first-order linear schemes onto nonlinear second-order schemes with non-oscillatory shock capturing capabilities, and a conservative formulation enabling to extend the schemes to nonlinear CLs. These three design steps, and the underlying theoretical tools, are discussed in depth. The nonlinear RD schemes resulting from this construction are tested on a large set of problems involving the solution of scalar models, and systems of CLs. This extensive verification fills the gaps left open, where no theoretical analysis is possible. <p>Numerical results are presented on the Euler equations of a perfect gas, on a two-phase flow model with highly nonlinear thermodynamics, and on the shallow-water equations. <p>On irregular grids, the schemes proposed yield quite accurate and stable solutions even on very difficult computations. Direct comparisone show that these results are more accurate than the ones given by FV and WENO schemes. Moreover, our schemes have a compact nearest-neighbor stencil. This encourages to further develop our approach, toward the design of very high-order schemes, which would represent a very appealing alternative, both in terms of accuracy and efficiency, to now classical FV and ENO/WENO discretizations. These schemes might also be very competitive with respect to very high-order DG schemes. / Doctorat en sciences appliquées / info:eu-repo/semantics/nonPublished

Le contexte de ce travail est le contrôle optimal géométrique appliqué à la mécanique céleste et au contrôle quantique. On s’est tout d’abord intéressé au problème de transfert orbital de satellite autour de la Terre à consommation minimale,qui amena à la réalisation du code HamPath, permettant tout d’abord la résolution de problème de contrôle optimal dont la loi de commande est lisse. Il se base sur le Principe du Maximum de Pontryagin (PMP) et sur la notion de point conjugué. Ce programme combine méthode de tir, méthodes homotopiques différentielles et calcul des conditions d’optimalité du deuxième ordre. Nous nous intéressons par la suite au contrôle quantique. On étudie tout d’abord le contrôle d’un système composé de deux types de particules de spin 1/2 ayant des temps de relaxation différents et dont la dynamique est gouvernée par les équations de Bloch. Ces deux sous-systèmes,correspondant aux deux types de particules, sont couplés par un même contrôle (un champ electromagnétique), le but étant alors d’amener la magnétisation des particules du premier type à zéro tout en maximisant celle du second (dans un système de coordonnées bien choisi). Ce modèle intervient en imagerie médicale par Résonance Magnétique Nucléaire et consiste à maximiser le contraste entre deux régions d’une même image. L’utilisation des outils géométriques et numériques aura permis de donner une très bonne synthèse sous-optimale pour deux cas particuliers (mélange sang oxygéné/désoxygéné et liquide cérébrospinal/eau). La dernière contribution de cette thèse porte sur l’étude d’un système quantique à deux niveaux d’énergie dontl a dynamique est régie par les équations de Lindblad. Le modèle est basé sur la minimisation d’énergie du transfert. On se restreint à un cas particulier pour lequelle Hamiltonien donné par le PMP est Liouville intégrable. On décrit alors les lieux conjugués et de coupure pour ce problème riemannien avec dérive / This work is about geometric optimal control applied to celestial and quantum mechanics. We first dealt with the minimum fuel consumption problem of transfering a satellite around the Earth. This brought to the creation of the code HamPath which permits first of all to solve optimal control problem for which the command law is smooth. It is based on the Pontryagin Maximum Principle (PMP) and on the notion of conjugate point. This program combines shooting method, differential homotopic methods and tools to compute second order optimality conditions. Then we are interested in quantum control. We study first a system which consists in two different particles of spin 1/2 having two different relaxation time. Both sub-systems are driven by the same control. The problem consists in bringing to zero the magnetization of one of the two system while maximizing the magnetization of the second one. This problem comes from constrast imaging in Nuclear Magnetic Resonance and consists in maximising the contrast between two areas of the image. The use of geometrical and numerical tools has given a very precise sub-optimal synthesis for two particular cases (deoxygenated/oxygenated blood and cerebrospinal fluid/water cases). The last contribution of this thesis is about the Lindblad equations in the two-level case. The model is based upon the minimisation of the transfer energy. We restrict the study to a particular case for which the Hamiltonian given by the PMP is Liouville integrable.We describe the conjugate and cut loci for this Riemannian with drift problem

Cette thèse est composée de deux parties indépendantes : la première partie traite des équations différentielles stochastiques dans le cadre de la G-espérance, tandis que la deuxième partie présente les résultats obtenus pour les équations différentielles stochastiques du seconde ordre. Dans un premier temps, on considère les intégrales stochastiques par rapport à un processus croissant, et on donne une extension de la formule d'Itô dans le cadre de la G-espérance. Ensuite, on étudie une classe d'équations différentielles stochastiques réfléchies unidimensionnelles dirigées par un G-mouvement brownien. Dans la suite, en utilisant une méthode de localisation, on prouve l'existence et l'unicité de solutions pour les équations différentielles stochastiques dirigées par un G-mouvement brownien, dont les coefficients sont localement lipschitziens. Enfin, dans le même cadre, on discute des problèmes de réflexion multidimensionnelle et on fournit quelques résultats de convergence. Dans un deuxième temps, on étudie une classe d'équations différentielles stochastiques rétrogrades du seconde ordre à croissance quadratique. Le but de ce travail est de généraliser le résultat obtenu par Possamaï et Zhou en 2012. On montre aussi l'existence et l'unicité des solutions pour ces équations, mais sous des hypothèses plus faibles. De plus, ce résultat théorique est appliqué aux problèmes de maximisation robuste de l'utilité du portefeuille en finance. / This thesis consists of two relatively independent parts : the first part concerns stochastic differential equations in the framework of the G-expectation, while the second part deals with a class of second order backward stochastic differential equations. In the first part, we first consider stochastic integrals with respect to an increasing process and give an extension of Itô's formula in the G-framework. Then, we study a class of scalar valued reflected stochastic differential equations driven by G-Brownian motion. Subsequently, we prove the existence and the uniqueness of solutions for some locally Lipschitz stochastic differential equations driven by G-Brownian motion. At the end of this part, we consider multidimensional reflected problems in the G-framework, and some convergence results are obtained. In the second part, we study the wellposedness of a class of second order backward stochastic differential equations (2BSDEs) under a quadratic growth condition on their coefficients. The aim of this part is to generalize a wellposedness result for quadratic 2BSDEs by Possamaï and Zhou in 2012. In this thesis, we work under some usual assumptions and deduce the existence and uniqueness theorem as well. Moreover, this theoretical result for quadratic 2BSDEs is applied to solve some robust utility maximization problems in finance.

L'objectif de cette thèse est l'étude de la représentation probabiliste des différentes classes d'EDPSs non-linéaires(semi-linéaires, complètement non-linéaires, réfléchies dans un domaine) en utilisant les équations différentielles doublement stochastiques rétrogrades (EDDSRs). Cette thèse contient quatre parties différentes. Nous traitons dans la première partie les EDDSRs du second ordre (2EDDSRs). Nous montrons l'existence et l'unicité des solutions des EDDSRs en utilisant des techniques de contrôle stochastique quasi- sure. La motivation principale de cette étude est la représentation probabiliste des EDPSs complètement non-linéaires. Dans la deuxième partie, nous étudions les solutions faibles de type Sobolev du problème d'obstacle pour les équations à dérivées partielles inteégro-différentielles (EDPIDs). Plus précisément, nous montrons la formule de Feynman-Kac pour l'EDPIDs par l'intermédiaire des équations différentielles stochastiques rétrogrades réfléchies avec sauts (EDSRRs). Plus précisément, nous établissons l'existence et l'unicité de la solution du problème d'obstacle, qui est considérée comme un couple constitué de la solution et de la mesure de réflexion. L'approche utilisée est basée sur les techniques de flots stochastiques développées dans Bally et Matoussi (2001) mais les preuves sont beaucoup plus techniques. Dans la troisième partie, nous traitons l'existence et l'unicité pour les EDDSRRs dans un domaine convexe D sans aucune condition de régularité sur la frontière. De plus, en utilisant l'approche basée sur les techniques du flot stochastiques nous démontrons l'interprétation probabiliste de la solution faible de type Sobolev d'une classe d'EDPSs réfléchies dans un domaine convexe via les EDDSRRs. Enfin, nous nous intéressons à la résolution numérique des EDDSRs à temps terminal aléatoire. La motivation principale est de donner une représentation probabiliste des solutions de Sobolev d'EDPSs semi-linéaires avec condition de Dirichlet nul au bord. Dans cette partie, nous étudions l'approximation forte de cette classe d'EDDSRs quand le temps terminal aléatoire est le premier temps de sortie d'une EDS d'un domaine cylindrique. Ainsi, nous donnons les bornes pour l'erreur d'approximation en temps discret. Cette partie se conclut par des tests numériques qui démontrent que cette approche est effective. / The objective of this thesis is to study the probabilistic representation (Feynman-Kac for- mula) of different classes ofStochastic Nonlinear PDEs (semilinear, fully nonlinear, reflected in a domain) by means of backward doubly stochastic differential equations (BDSDEs). This thesis contains four different parts. We deal in the first part with the second order BDS- DEs (2BDSDEs). We show the existence and uniqueness of solutions of 2BDSDEs using quasi sure stochastic control technics. The main motivation of this study is the probabilistic representation for solution of fully nonlinear SPDEs. First, under regularity assumptions on the coefficients, we give a Feynman-Kac formula for classical solution of fully nonlinear SPDEs and we generalize the work of Soner, Touzi and Zhang (2010-2012) for deterministic fully nonlinear PDE. Then, under weaker assumptions on the coefficients, we prove the probabilistic representation for stochastic viscosity solution of fully nonlinear SPDEs. In the second part, we study the Sobolev solution of obstacle problem for partial integro-differentialequations (PIDEs). Specifically, we show the Feynman-Kac formula for PIDEs via reflected backward stochastic differentialequations with jumps (BSDEs). Specifically, we establish the existence and uniqueness of the solution of the obstacle problem, which is regarded as a pair consisting of the solution and the measure of reflection. The approach is based on stochastic flow technics developed in Bally and Matoussi (2001) but the proofs are more technical. In the third part, we discuss the existence and uniqueness for RBDSDEs in a convex domain D without any regularity condition on the boundary. In addition, using the approach based on the technics of stochastic flow we provide the probabilistic interpretation of Sobolev solution of a class of reflected SPDEs in a convex domain via RBDSDEs. Finally, we are interested in the numerical solution of BDSDEs with random terminal time. The main motivation is to give a probabilistic representation of Sobolev solution of semilinear SPDEs with Dirichlet null condition. In this part, we study the strong approximation of this class of BDSDEs when the random terminal time is the first exit time of an SDE from a cylindrical domain. Thus, we give bounds for the discrete-time approximation error.. We conclude this part with numerical tests showing that this approach is effective.

Cette thèse s'intéresse aux liens entre l'évolution de la coopération et la sélection de second ordre. Dans une première partie, nous montrons comment des organismes digitaux adaptent leurs génomes pour encoder les gènes liées à la coopération d'une manière plus contrainte (suppression d'évolvabilité), notamment à l'aide d'opérons et d'overlaps impliquant aussi des gènes essentiels. Dans une deuxième partie, nous testons expérimentalement cette vision des overlaps de gènes comme "contrainte évolutive" grâce à des outils d'algorithmique et de biologie synthétique que nous avons développés. Dans une troisième partie, nous utilisons des simulations par agents pour montrer comment une forme de division du travail peut être interprétée comme un système coopératif à la lumière de la théorie évolutive moderne. Dans une dernière partie, nous montrons que la dispersion spatiale des allèles coopératives obtenue par des phénomènes de "genetic hitchiking" joue un rôle important dans l'évolution de la coopération, quand bien même ce mécanisme de dispersion s'applique aussi à des allèles non coopératives, grâce à la "relatedness" (aux loci codant pour la coopération) crée par l'invasion locale de mutations bénéfiques (à des loci non liés à la coopération) et par l'équilibre complexe entre ces mutations bénéfiques et la robustesse mutationnelle. L'ensemble de ces résultats appelle à une prise en compte plus importante des pressions sélectives de second ordre dans l'étude de l'évolution sociale, et au développement de modèles plus réalistes qui permettraient d'intégrer de telles forces évolutives. Nous insistons également sur l'importance du paysage mutationnel dans l'étude des populations bactériennes, et montrons le potentiel croissant de la biologie synthétique comme outil d'étude de ce paysage et de l'évolution microbienne en général. / In the first part, I show how digital organisms adapt their genomes to encode cooperation-related genes in a more constrained way (evolvability suppression), especially using operons and overlaps also involving essential genes. In the second part, we experimentally test this view of gene overlaps as an evolutionary constraint, using both algorithmic and synthetic biology tools that we have developed. In the third part, I use agent-based simulations to show how a form of division of labour can be interpreted as a cooperative system in the light of modern evolutionary theory. In the final part, I show that the patterns of dispersal of cooperative alleles due to hitchhiking phenomena play an important role in the evolution of cooperation. The last result holds even though the hitchhiking mechanisms also applies to non-cooperative alleles, thanks to the relatedness (at cooperation-related loci) created by the local invasion of beneficial mutations (at loci not related to cooperation). The beneficial mutations form a complex and interesting equilibrium with mutational robustness, which I investigate using in silico evolution. On the whole, these results call for a more careful consideration of the second-order selection pressures in the study of social evolution, and show the necessity for more realistic models allowing to integrate such evolutionary forces. My thesis research specifically highlights the importance of the mutational landscape in the study of microbial populations and shows the increasing potential of synthetic biology as a tool to study such landscape and microbial evolution in general.

Dans un contexte de densification des villes et de leurs réseaux de transport, les gens sont de plus en plus exposés au bruit. Ainsi, le résultat de l'étude d'impact vibro-acoustique joue un rôle primordial dans l'expansion du réseau ferroviaire. L'une des principales sources est le bruit de roulement : La rugosité de la surface de la roue et du rail produit un déplacement imposé sur ces derniers. Ce déplacement entraine une réponse vibratoire des roues et de la voie ferrée et leurs rayonnements acoustiques. Cette thèse propose une amélioration de la modélisation vibro-acoustique de la voie ferrée.Pour la réponse vibratoire, le coté infini de la voie et sa déformation dans les 3 dimensions rendent les modèles analytiques et les éléments finis non-optimales dans la gamme de fréquence de l’audible. La méthode élément fini semi-périodique (SAFEM) est utilisée dans cette thèse pour modéliser une voie à support continue. Elle est ensuite couplée au théorème de Floquet pour modéliser une voie à support périodique. Cependant, cette technique génère des problèmes numériques qui ont imposé un algorithme adapté. La méthode d'Arnoldi du second ordre (SOAR) est utilisée avant de résoudre l'équation SAFEM permet de résoudre ces problèmes ainsi qu’apporter la stabilité requise. Des comparaisons avec d’autres modèles et des données expérimentales permettent de valider la méthode.Pour le rayonnement acoustique, la simulation de grand domaine en haute fréquence rendent inadapté l'utilisation de techniques conventionnelles (FEM, BEM, ...). La méthode proposée ici : la théorie variationnelle du rayon complexe est particulièrement bien adaptée à ce cas. Les principales caractéristiques de l'approche VTCR sont l'utilisation d'une formulation faible du problème acoustique, qui permet de considérer automatiquement les conditions limites entre sous-domaines. Ensuite, l'utilisation d'une répartition intégrale des ondes planes dans toutes les directions permet de simuler le champ acoustique. Les inconnues du problème sont leurs amplitudes. Cette méthode qui a déjà montré son efficacité pour les domaines fermés a été étendue au domaine ouvert et couplée à la réponse vibratoire. Des comparaisons avec des solutions analytiques et des simulations FEM à basse fréquence permettent de valider la méthode. / In a context of urban and transport network densification, people are increasingly exposed to noise. Consequently, the result of vibro-acoustic impact assessment has a pivotal role in rail network expansion. One of the main sources is the rolling noise: Roughness on the wheel and rail surface produce an imposed displacement one the both. This last, generates vibrational response of wheels and the railway track and their acoustic radiation. This PhD thesis presents some improvements of the vibro-acoustic railway track modelling.Concerning vibrational response, the infinite dimension in the longitudinal direction of the track and its deformation in the 3 dimensions, make the analytical models and finite elements non-optimal. The Semi-analytical finite element method (SAFEM), used in this thesis, is particularly well adapted in this case. Firstly, it is used to model railway track on a continuous support. Then, it is coupled with Floquet theorem to model tracks with a periodic support. However, this technique suffers from numerical problems that imposed an adapted algorithm. The second-order Arnoldi method (SOAR) is used to tackle them. This reduction allows to eliminate critical values improving the robustness of the method. Comparison with existing techniques and experimental results validate this model.Concerning acoustic radiation, big domains simulations at high frequency are almost unfeasible when using conventional techniques (FEM, BEM,…). The method used in this thesis, the Variational theory of complex ray (VTCR) is particularly well adapted to these cases. The principal features of VTCR approach are the use of a weak formulation of the acoustic problem, which allows to consider automatically boundary conditions between sub-domains. Then, the use of an integral repartition of plane waves in all the direction allow to simulate the acoustic field. The unknowns of the problem are their amplitudes. This method well assessed for closed domain, has been extended to open domain and coupled to vibrational response of the rail. Comparison with analytic solution and FEM simulation at low frequency allow to validate the method.Coupling these both methods allowed to simulate complex real life vibro-acoustic scenarios. Result of different railway tracks are presented and validated

The diploma thesis is focused on design and optimalization of digital octave and fraction-octave band filters. This thesis describe the behavior of filters in systems with fixed point arithmetics and investigate the impact of quantization coefficients for frequency response of filter. Filters, whitch has been designed, are implemented into simple software in C. Designed filters are in accordance with standard IEC 61260.

Cette thèse s'inscrit dans l'étude mathématique de la percolation, qui regroupe une famille de modèles présentant une transition de phase. Des avancées majeures au cours des quinze dernières années, notamment l'invention du SLE et la preuve de l'invariance conforme de la percolation de Bernoulli critique, nous permettent aujourd'hui d'avoir une image très complète de la percolation de Bernoulli sur le réseau triangulaire. Cependant, de nombreuses questions demeurent ouvertes, et ont motivé notre travail.La première d'entre elle est l'universalité de la percolation plane, qui affirme que les propriétés macroscopiques de la percolation plane critique ne devraient pas dépendre du réseau sous-jacent à sa définition. Nous montrons, dans le cadre de la percolation Divide and Color, un résultat qui va dans le sens de cette universalité et identifions, dans ce contexte, des phénomènes macroscopiques indépendants du réseau microscopique. Une version plus faible d'universalité est donnée par la théorie de Russo-Seymour-Welsh (RSW), et sa validité est connue pour la percolation de Bernoulli (sans dépendance) sur les réseaux plans suffisamment symétriques. Nous étudions de nouveaux arguments de type RSW pour des modèles de percolation avec dépendance. La deuxième question que nous avons abordée est celle de l'absence d'une composante connexe ouverte infinie au point critique, une question importante du point de vue physique, puisqu'elle traduit la continuité de la transition de phase. Dans deux travaux en collaboration avec Hugo Duminil-Copin et Vladas Sidoravicius, nous montrons que la transition de phase est continue pour la percolation de Bernoulli sur le graphe Z^2x{0,...,k}, et pour la percolation FK sur le réseau carré avec paramètre q inférieur ou égal à 4. Enfin, la dernière question qui nous a guidés est la localité du point critique : la donnée des boules de grands rayons d'un graphe suffit-elle à identifier avec une bonne précision la valeur du point critique? Dans un travail en collaboration avec Sébastien Martineau, nous répondons de manière affirmative à cette question dans le cadre des graphes de Cayley de groupes abéliens. / This thesis is part of the mathematical study of percolation theory, which includes a family of models with a phase transition. Major advances in the 2000s, including the invention of SLE and the proof of conformal invariance of critical Bernoulli percolation, provide us with a very complete picture of the Bernoulli percolation process on the triangular lattice. Fortunately, many questions remain open, and motivated our work.The first of these is the universality of planar percolation, which states that the macroscopic properties of critical planar percolation should not depend on the underlying graph. We study this question in the framework of Divide and Color percolation, and prove in this context a result that goes in the direction of universality. A weaker universality statement is given by the theory of Russo-Seymour-Welsh (RSW), which is known to hold for planar Bernoulli percolation (without dependence) on sufficiently symmetric graphs. We study new RSW-type arguments for percolation models with dependence.The second question is the absence of an infinite cluster at the critical point, an important question from a physical point of view, equivalent to the continuity of the phase transition. In two different joint works with Hugo Duminil-Copin and Vladas Sidoravicius, we show that the phase transition is continuous for Bernoulli percolation on the graph Z^2 x {0,...,k} and for FK percolation on the square lattice with parameter q smaller than or equal to 4.Finally, the last question that guided us is the locality of the critical point: is it possible to determine with good accuracy the critical value for Bernoulli percolation on a graph if we know only the balls with large radii? Jointly with Sébastien Martineau, we answer positively to this question in the framework of Cayley graphs of abelian groups.