On s’intéresse dans cette thèse à montrer que la solution approchée, par la méthode des volumes finis, converge vers la solution renormalisée de problèmes elliptiques ou paraboliques à donnée L1. Dans la première partie nous étudions une équation de convection-diffusion ellliptique à donnée L1. En adaptant la stratégie développée pour les solutions renormaliséesà la méthode des volumes finis, nous montrons que la solution approchée converge vers l’unique solution renormalisée.Dans la deuxième partie nous nous intéressons à un problème parabolique nonlinéaire à donnée L1. En utilisant une version discrète de résultats de compacité classiques, nous montrons que les résultats obtenues dans le cas elliptique restentvrais dans le cas parabolique. Dans la troisième partie nous montrons des résultats similaires pour une équationparabolique doublement non-linéaire à donnée L1. Le caractère doublement nonlinéaire de l’équation crée des difficultés supplémentaires par rapport à la partie précédente, notamment car la règle de dérivation en chaîne ne s’applique pas dansle cas discret. Enfin, dans la quatrième partie, nous utilisons les résultats établis précédemment pour étudier un système de type thermoviscoélasticité. Nous montrons que la solution approchée, obtenue par un schéma éléments finis-volumes finis, converge vers une solution faible-renormalisée du système. / In this thesis we are interested in proving that the approximate solution, obtained by the finite volume method, converges to the unique renormalized solution of elliptic and parabolic equations with L1 data. In the first part we study an elliptic convection-diffusion equation with L1 data. Mixing the strategy developed for renormalized solution and the finite volume method,we prove that the approximate solution converges to the unique renormalized solution. In the second part we investigate a nonlinear parabolic equation with L1 data. Using a discrete version of classical compactness results, we show that the results obtaines previously in the elliptic case hold true in the parabolic case. In the third part we prove similar results for a doubly nonlinear parabolic equation with L1 data. The doubly nonlinear character of the equation makes new difficulties with respect to the previous part, especially since the chain rule formula does not apply in the discrete case. Finaly, in the fourth part we use the results established previously to investigate a system of thermoviscoelasticity kind. We show that the approximate solution,obtaines by finite element-finite volume scheme, converges to a weak-renormalized solution of the system.

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 fi nies 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ée 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 a l'interface, afin de imposer des conditions de Neumann homogènes sur le bord 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.

L'objectif de ce travail est de prendre en compte les instationnarités naturelles à grande échelle dans les écoulements décollés et à un coût plus faible que la LES, tout en s'intéressant à la modélisation des effets de paroi par des modèles statistiques au second ordre. S'inspirant des approches de Durbin, le modèle à pondération elliptique EB-RSM reproduit l'effet non-local de blocage, en résolvant une équation différentielle sur le terme de pression. La limite à deux composantes de la turbulence est bien prédite en canal. Ce modèle est appliqué à la marche descendante, dans une approche URANS. Nous avons montré que les erreurs numériques peuvent être suffisantes pour exciter le mode le plus instable de la couche cisaillée, et aboutir à une solution instationnaire. La solution est stationnaire quand on raffine le maillage, rendant l'URANS peu fiable. Récemment, Schiestel \& Dejoan ont proposé le modèle hybride non-zonal PITM. Le coefficient $C_{\e_2}$ de l'équation de la dissipation devient fonction de la coupure dans le spectre, et la valeur $C_{\e_1}=3/2$ est déduite par ces auteurs. Nous avons donné une formulation plus générale où la valeur de $C_{\e_1}$ est quelconque. Pour offrir un formalisme plus cohérent aux modèles hybrides non-zonaux dans les écoulements de paroi, une approche basée sur un filtrage temporel est proposée. Enfin, l'adaptation du modèle EB-RSM dans un cadre hybride a été réalisée. Les résultats en canal sont encourageants : la transition continue d'un modèle RANS en proche paroi à une LES au centre du canal est mise en évidence. Le transfert d'énergie des échelles modélisées vers celles résolues est bien reproduit quand on raffine le maillage.

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.

L'idée d'utiliser des bases d'ondelettes dans l'analyse numérique (résolution des équations elliptiques, aux dérivées partielles, intégrales) s'est imposée depuis que ces bases ont fait preuve de leur efficacité dans le traitement du signal. Deux problèmes se posent quant au calcul de la solution dans une base d'ondelettes : - problème 1 : l'étude de la structure de la matrice associée à un noyau K d'un opérateur intégral T dans une base d'ondelettes ; - problème 2 : l'adaptation des techniques de discrétisation de Galerkin aux bases d'ondelettes. Cette thèse contribue à l'étude de ces problèmes par l'introduction d'une nouvelle classe d'opérateurs définis par leur matrice représentative dans une base d'ondelettes et caractérisés par les dérivées fractionnaires de leurs noyaux.

La presente recherche doctorale en analyse numerique (et calcul scientifique) aborde le probleme du controle de la dissipation d'entropie numerique associee a une discretisation donnee. Cette problematique constitue un veritable challenge numerique non encore completement resolu a ce jour.<br /><br />Le travail se decompose en deux parties principales dont les caracteristiques sont reellement differentes.<br /><br />La premiere partie concerne l'approximation numerique des solutions (instationnaires en 1D et stationnaires en 2D) du systeme des equations de Navier-Stokes a plusieurs pressions independantes. Ce systeme est hyperbolique et possede des champs vraiment non lineaires sous des hypotheses classiques, mais s'ecrit naturellement sous forme non conservative.<br /><br />La deuxieme partie est dediee a l'approximation numerique des solutions instationnaires en 1D de quelques systemes de lois de conservation de type soit hyperbolique mais dont les champs possedent un defaut de vraiment non linearite, ou soit mixte hyperbolique-elliptique.<br /><br />Dans toutes ces situations motivees par des applications physiques concretes, le controle de la dissipation d'entropie joue un role déterminant dans la caracterisation des solutions recherchees. Les schemas numeriques proposes dans ce manuscrit sont obtenus par une analyse fine des bilans d'entropie associes.

Ce travail se compose de deux parties indépendantes. La première est consacrée à l'étude de la méthode du décalage, dite aussi méthode du Shift, pour le calcul d'espérances mathématiques en dimension grande ou infinie. Pour l'essentiel, la méthode du décalage est la mise en oeuvre informatique du théorème ergodique ponctuel de Birkhoff pour l'opérateur de décalage (à gauche ou, à défaut, à droite). La deuxième partie s'attache au problème de l'approximation des mesures invariantes pour les chaînes de Markov.

Cette thèse porte sur l'étude, la conception matérielle, la validation théorique et pratique, et enfin la comparaison de différents opérateurs arithmétiques pour des cryptosystèmes basés sur les courbes elliptiques (ECC). Les solutions proposées doivent être robustes contre certaines attaques par canaux cachés tout en étant performantes en matériel, tant au niveau de la vitesse d'exécution que de la surface utilisée. Dans ECC, nous cherchons à protéger la clé secrète, un grand entier, utilisé lors de la multiplication scalaire. Pour nous protéger contre des attaques par observation, nous avons utilisé certaines représentations des nombres et des algorithmes de calcul pour rendre difficiles certaines attaques ; comme par exemple rendre aléatoires certaines représentations des nombres manipulés, en recodant certaines valeurs internes, tout en garantissant que les valeurs calculées soient correctes. Ainsi, l'utilisation de la représentation en chiffres signés, du système de base double (DBNS) et multiple (MBNS) ont été étudiés. Toutes les techniques de recodage ont été validées théoriquement, simulées intensivement en logiciel, et enfin implantées en matériel (FPGA et ASIC). Une attaque par canaux cachés de type template a de plus été réalisée pour évaluer la robustesse d'un cryptosystème utilisant certaines de nos solutions. Enfin, une étude au niveau matériel a été menée dans le but de fournir à un cryptosystème ECC un comportement régulier des opérations effectuées lors de la multiplication scalaire afin de se protéger contre certaines attaques par observation. / This PhD thesis focuses on the study, the hardware design, the theoretical and practical validation, and eventually the comparison of different arithmetic operators for cryptosystems based on elliptic curves (ECC). Provided solutions must be robust against some side-channel attacks, and efficient at a hardware level (execution speed and area). In the case of ECC, we want to protect the secret key, a large integer, used in the scalar multiplication. Our protection methods use representations of numbers, and behaviour of algorithms to make more difficult some attacks. For instance, we randomly change some representations of manipulated numbers while ensuring that computed values are correct. Redundant representations like signed-digit representation, the double- (DBNS) and multi-base number system (MBNS) have been studied. A proposed method provides an on-the-fly MBNS recoding which operates in parallel to curve-level operations and at very high speed. All recoding techniques have been theoretically validated, simulated extensively in software, and finally implemented in hardware (FPGA and ASIC). A side-channel attack called template attack is also carried out to evaluate the robustness of a cryptosystem using a redundant number representation. Eventually, a study is conducted at the hardware level to provide an ECC cryptosystem with a regular behaviour of computed operations during the scalar multiplication so as to protect against some side-channel attacks.

Cette thèse porte principalement sur l'étude des solutions de certaines équations aux dérivées partielles elliptiques via l'indice de Morse, y compris des solutions stables, i.e. quand l'indice de Morse est égal à zéro. Elle comporte deux parties indépendantes.Dans la première partie, sous des hypothèses sur-linéaires et sous-critiques sur f, on établit d'abord une estimation explicite de la norme L [infini] des solutions de -Δu = f(u) avec u = 0 sur le bord, via leurs indices de Morse. On propose une approche plus transparente et plus souple que le travail de Yang [1998], ce qui nous permet de traiter des non linéarités très proches de la croissance critique. Les résultats obtenus nous ont motivé de travailler sur des équations polyharmoniques (-Δ)ku = f(x; u) avec notamment k = 2 et 3. Avec des hypothèses semblables à Yang [1998] sur f et des conditions au bord convenables, on obtient pour la première fois des estimations explicites de solution des équations polyhamoniques, via l'indice de Morse. Dans la seconde partie, on considère un système de Lane-Emden-Δu = ρ(x)vp; -Δv = ρ(x)u θ ; u; v > 0; dans RN; avec 1 < p< θ et un poids radial ρ strictement positif. Nous montrons la non-existence de solution stable en petites dimensions N. Nos résultats améliorent les travaux précédents de Cowan & Fazly [2012]; Fazly [2012]; Hu [2015], et fournissent notamment des résultats du type Liouville pour solution stable, en petites dimensions N, valables pour tout 1 < ρ min(4 3 ; θ) / The main concern of this thesis deals with the study of solutions of several elliptic partial differential equations via the Morse index, including the stable solutions, i.e. when the Morse index is zero. The thesis has two independent parts. In the first part, under suplinear and subcritical assumptions on f, we establish firstly some explicit estimation for the L1 norms of solutions to -Δu = f(u) avec u = 0 on the boundary, via its Morse index. We propose an approach more transparent and easier than the work of Yang [1998], which allow us to treat some nonlinearities very close to the critical growth. These results motivated us to consider the polyharmonic equations (-Δ)ku = f(x; u) with especially k = 2 and 3. With the hypothesis on f similar to Yang [1998] and appropriate boundary conditions, we obtain for the _rst time some explicit estimations of solution via its Morse index, for the polyharmonic equations.In the second part, we consider a Lane-Emden system -Δu = ρ(x)vp; -Δv = ρ(x)u_; u; v > 0; in RN; with 1 < p< θ and a radial positive weight ρ. We prove the non-existence of stable solution in small dimension case. Our results improve the previous works Cowan & Fazly [2012]; Fazly [2012]; Hu [2015], especially we prove some general Liouville type results for stable solutions in small dimension which hold true for any 1 < ρ min(4 3 ; θ)

Les conditions extrêmes de température et de densité d'énergie créées lors des collisions d'ions lourds ultra-relativistes au sein du Grand collisionneur de hadrons (LHC) fournissent une occasion unique d'étudier les propriétés de la matière. Une transition de phase de la matière hadronique vers un milieu déconfiné de quarks et de gluons (PQG) est prédite par la chromodynamique quantique et des efforts théoriques et expérimentaux considérables ont été investis pour étudier ses propriétés. Parmi les sondes éminentes du PQG, les quarks lourds jouent un rôle prépondérant car ils sont créés lors des processus durs initiaux, avant la formation du PQG, et leur nombre est conservé durant les phases partoniques et hadroniques de la collision. Les états liés de quarks lourds– quarkonium (charmonium pour cc et bottomonium pour bb), constituent des sondes remarquables du milieu. Les observations expérimentales de quarkonia aux énergies du LHC dans les collisions A-A sont reproduites au travers de deux mécanismes antagonistes : la suppression séquentielle, proposée très tôt comme signature du PQG, et la (re)génération de quarkonia par (re)combinaison de quarks déconfinés. Cependant des incertitudes importantes sont associées aux prédictions théoriques et de nombreuses inconnues demeurent. L'anisotropie azimuthale de l'espace des moments (désignée sous le terme de flot elliptique, v2) de la production de charmonium devrait permettre d'avoir une meilleure vue d'ensemble et de contraindre davantage les paramètres des modèles théoriques. Si les quarks charmés se (re)combinent en paires cc, les J/ψ produits devraient hériter de leur flot. Les études précédentes ont montré des premiers signes d'un v2positif du J/ψ dans les collisions Pb-Pb à √sNN=2.76TeV. Cette thèse porte sur la mesure de l'anisotropie azimuthale du J/ψ dans les collisions Pb-Pb à√sNN=5.02TeV où une (re)combinaison légèrement plus importante est prédite par rapport aux énergies inférieures et une augmentation statistique d'un facteur 3 a été enregistrée. L'étude du flot elliptique du J/ψ fournit des informations clés sur la magnitude et la dynamique des mécanismes de suppression et de (re)génération de charmonia. Par ailleurs elle offre un regard unique sur l'évolution et les interactions des quarks charmés au sein du milieu en expansion. / Extreme temperatures and energy densities produced in ultra-relativistic heavy-ion collisions at the Large Hadron Collider provide a unique opportunity to study the properties of matter. A phase transition of the hadronic matter to a deconfined medium of quarks and gluons, the Quark-Gluon Plasma (QGP), is predicted by Quantum Chromodynamics and considerable theoretical and experimental efforts have been invested to study its properties. Among the prominent probes of the QGP, heavy quarks play a crucial role since they are created in primary hard-scattering processes, before the QGP formation, and their number is conserved throughout the partonic and hadronic phases of the collision. Bound states of heavyquarks – quarkonium (charmonium for cc and bottomoniumbb) provide remarkable probes of the medium. At LHC energies, experimental observations of quarkonium in A-Acollisions are reproduced through two antagonist mechanisms: a sequential suppression of the quarkonium states, early suggested as a signature of the QGP, and quarkonium (re)generation by (re)combination of deconfined quarks. However, theoretical predictions carry large uncertainties and many unknows remain. The momentum space azimuthal anisotropy of charmonium production (referred as elliptic flow v2) should help to clarify the picture and to constrain the model parameters. If charm quarks (re)combine in the medium into cc pairs, the J/ψ originating from (re)combination should inherit their flow. Previous studies have shown first hints of a positive J/ψ v2 in Pb-Pb collisions at √sNN=2.76TeV. This thesis focuses on the measurement of J/ψ azimuthal anisotropy in Pb-Pb collisions at √sNN=5.02TeV where a slightly stronger (re)generation component is predicted with respect to lower collision energies, and a factor of 3more data were collected. The study of J/ψ v2 provides important information on the magnitude and dynamics of charmonium suppression and (re)generation mechanisms. In addition, it offers a unique insight on charm quark evolution and interactions in the expanding medium.