Cette thèse traite de la reconstruction d'hypersurfaces au sein de champs de normales en dimension quelconque et trouve des applications dans l’analyse des empreintes digitales (lignes dermiques), des images satellites météorologiques (lieux de turbulence) et astrophysiques (bras de galaxies) ainsi que dans l’analyse stratigraphique des images sismiques (horizons). Les méthodes développées s’appuient sur la minimisation d’une équation aux dérivées partielles non linéaire reliant une hypersurface au pendage déduit d’un champ de normales. Elles prennent en compte des contraintes diverses telles que des points de passages, des frontières, des bornes et des discontinuités. La contribution principale de la thèse réside dans l’introduction d’un changement d’espace du pendage qui permet de reconstruire aussi bien des hypersurfaces exprimées sous des formes implicites dans les repères de définition des champs de normales que des horizons sismiques de manière rapide et interactive. Deux schémas de reconstruction d’horizons sismiques unidimensionnels présentant une discontinuité d’amplitude et de lieu inconnus sont également proposés. / This thesis deals with the reconstruction of hypersurfaces from a finite-dimensional normal vector field. Application scopes can be found in the analysis of fingerprints (epidermal ridges), meteorological images (eddies and cyclones), astrophysical images (galaxy arms) and in the stratigraphic analysis of seismic images (horizons). The hypersurfaces are obtained by solving a non-linear partial derivative equation relied on the local dip deduced from a normal vector field. Several constraints such as boundaries, bounds, points belonging to the hypersurface or discontinuities can be considered.The major contribution of this thesis consists in a local dip transformation which allows to reconstruct implicit hypersurfaces as well as seismic horizons by a fast and interactive method. Two schemes dedicated to the reconstruction of discontinuous one-dimensional seismic horizons are also proposed when the discontinuity location and jump are unknown.

Cette thèse présente des travaux sur l’étude et la calibration de modèles d’équations aux dérivées partielles pour la croissance tumorale. La première partie porte sur l’analyse d’un modèle de croissance tumorale pour le cas de métastases au foie de tumeurs gastro-intestinales (GIST). Le modèle est un système d’équations aux dérivées partielles couplées et prend en compte plusieurs traitements dont un traitement anti-angiogénique. Le modèle permet de reproduire des données cliniques. La première partie de ce travail concerne la preuve d’existence/unicité de la solution du modèle. La seconde partie du travail porte sur l’étude du comportement asymptotique de la solution du modèle lorsqu’un paramètre du modèle, décrivant la capacité de la tumeur à évacuer la nécrose, converge vers 0. La seconde partie de la thèse concerne le développement d’un modèle de croissance pour des sphéroïdes tumoraux ainsi que sur la calibration de ce modèle à partir de données expérimentales in vitro. L’objectif est de développer un modèle permettant de reproduire quantitativement la distribution des cellules proliférantes à l’intérieur d’un sphéroïde en fonction de la concentration en nutriments. Le travail de modélisation et de calibration du modèle a été effectué à partir de données expérimentales permettant d’obtenir la répartition spatiale de cellules proliférantes dans un sphéroïde tumoral. / In this thesis, we present several works on the study and the calibration of partial differential equations models for tumor growth. The first part is devoted to the mathematical study of a model for tumor drug resistance in the case of gastro-intestinal tumor (GIST) metastases to the liver. The model we study consists in a coupled partial differential equations system and takes several treatments into account, such as a anti-angiogenic treatment. This model is able to reproduce clinical data. In a first part, we present the proof of the existence/uniqueness of the solution to this model. Then, in a second part, we study the asymptotic behavior of the solution when a parameter of this model, describing the capacity of the tumor to evacuate the necrosis, goes to 0. In the second part of this thesis, we present the development of model for tumor spheroids growth. We also present the model calibration thanks to in vitro experimental data. The main objective of this work is to reproduce quantitatively the proliferative cell distribution in a spheroid, as a function of the concentration of nutrients. The modeling and calibration of this model have been done thanks to experimental data consisting of proliferative cells distribution in a spheroid.

Cette thèse conclut 3 années d'études doctorales sur le "dépôt de couches atomiques (ALD) de nitrure de bore (BN)". Le but de ce travail a été d'adapter la voie des céramiques dérivées de polymères (PDC) à la technique ALD pour la croissance de films minces de h-BN et l'élaboration de nanostructures fonctionnelles. Tout d'abord, un nouveau procédé d'ALD sans ammoniac en deux étapes, comprenant une croissance par ALD à basse température (80 °C) de polyborazine (PBN) à partir de 2,4,6-trichloroborazine et d'hexaméthyldisilazane suivi un traitement thermique à haute température sous atmosphère contrôlée a été développé. Ainsi, des films minces uniformes et homogènes de BN ont pu être déposés sur divers substrats. Le caractère autolimité des réactions mises en jeu ainsi que l'homogénéité des films sur des supports très structurés ont été vérifiés. De ce fait des nanostructures fonctionnelles BN ont été réalisées à partir de substrats ou de templates de dimensionnalité variée. Leurs applications en tant que revêtements protecteurs et comme filtres et éponges absorbantes pour purifier les eaux polluées par des hydrocarbures ont en particulier été étudiées. Enfin, un deuxième procédé ALD basse température (85-150°C) utilisant le tri(isopropylamino)borane et la méthylamine comme précurseurs a été préalablement étudié afin de confirmer l'adaptabilité de la voie PDC et la technique ALD. Des films minces de BN ont été obtenus sur des substrats plans et il a été prouvé que les vapeurs de tri(isopropylamino)borane peuvent infiltrer des fibres de polyacrylonitrile électrofilées.Ce travail a été entièrement réalisé à l'Université de Lyon et a reçu le soutien financier du China Scholarship Council (CSC) pour la bourse de doctorat ainsi que de l'Agence Nationale de la Recherche (projet n° ANR-16-CE08-0021-01) / This thesis achieves 3 years of PhD studies on “Atomic layer deposition (ALD) of boron nitride (BN)”. The aim of this PhD work is to adapt the polymer derived ceramics (PDCs) route to the ALD technique for h-BN thin film growth and elaboration of functional nanostructures. A novel two-step ammonia-free ALD process, which includes ALD deposition of polyborazine at low temperature (80 °C) from 2,4,6-trichloroborazine and hexamethyldisilazane followed by post heat treatment under controlled atmosphere, has been established. Conformal and homogeneous BN thin films have been deposited onto various substrates. The self-limitation of the reactions on flat substrates and the conformality of the films on structured substrates have been verified. Functional BN nanostructures have thus been fabricated using substrates or templates with different dimensionalities. In particular, their applications as protective coatings as well as filter and absorber to purify polluted water from organic/oil hav e been investigated. Finally, a second low temperature (85-150 °C) ALD process using tri(isopropylamine)borane and methylamine as precursors has preliminary been studied in order to confirm the adaptability of PDCs route to ALD technique. BN thin films have been grown onto flat substrate and it has been proven that tri(isopropylamino)borane vapor can infiltrate into electrospun polyacrylonitrile fibers.This work was carried out at University of Lyon and financially supported by the National Research Agency (project n° ANR-16-CE08-0021-01)

The numerical modelling of the VKI-Longshot facility remains a challeng-ing task as it requires multi-physical numerical methods in order to simulate all the components. In the current dissertation, numerical tools were developed in order to study each component of the facility separately and a deep investigations of each stage of the shot were performed. This helped to better understand the different processes involved in the flow development inside this hypersonic wind tunnel. However the numerical computation of different regions of the facility treated as independent from each others remains an approximation at best.The accuracy of the rebuilding code for determining the free stream conditions and the total enthalpy in the VKI-Longshot facility was investigated by using a series of unsteady numerical computations of axisymmetric hypersonic flow over a heat flux probe. Good agreement was obtained between the numerical results and the measured data for both the stagnation pressure and the heat flux dur- ing the useful test time.The driver-driven part of the Longshot facility was modelled using the quasi one-dimensional Lagrangian solver L1d2. The three main conditions used for the experiments —low, medium and high Reynolds number —were considered.The chambrage effect due to the junction between the driver and the driven tubes in the VKI-Longshot facility was investigated. The computation showed great ben- efit of the chambrage in increasing the speed of the piston and thus the final compression ratio of the test gas.Two dimensional simulations of the flow in the driver and the driven tube were performed using Arbitrary Lagrangian Eulerian (ALE) solver in COOLFLuiD. A parallel multi-domain strategy was developed in order to integrate the moving piston within the computational domain.The computed pressure in the reservoir is compared to the one provided by the experiment and good agreement was obtained for both con- editions.Finally, an attempt was made to compute the starting process of the flow in the contoured nozzle. The transient computation of the flow showed how the primary shock initiates the flow in the nozzle before reaching the exit plan at time of 1.5 [ms] after the diaphragm rupture. The complex interactions of the reflected shocks in the throat raise the temperature above 9500 [K] which was not expected. Chemical dissociation of Nitrogen was not taken into account during this transient investigation which may play a key role considering the range of temperature reached near the throat. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Dans cette thèse, on construit des systèmes d’Euler à partir des unités (conjecturales) de Stark et celles de Rubin-Stark d’un corps de nombres K, pour décrire l’idéal caractéristique du X-quotient du module d’Iwasawa standard X∞ pour certains caractères p-adiques irréductibles X. Ici X∞ est le groupe de Galois de la pro-p-extension abélienne non ramifiée maximale de K∞, où K∞ est une Zp-extension adéquate de K. Plus précisément, on démontre des résultats de divisibilité formulée par la conjecture principale de la théorie d’Iwasawa. Nos démonstrations reposent essentiellement sur la théorie des systèmes d’Euler. / In this thesis, we construct Euler systems coming from the (conjectural) Stark units and those of Rubin-Stark of a number field K, to describe the characteristic ideal of the X-quotient of the standard Iwasawa module X∞, for some p-adic irreducible characters X. Here X∞ is the Galois group of the maximal unramified abelian pro-p-extension of K∞, where K∞ is an adequate Zp-extension of K. Precisely, we demonstrate a divisibility results formulated by the main conjecture of Iwasawa theory. Our demonstrations essentially are based on the theory of Euler systems.

Etude des propriétés de symétrie des solutions de quelques problèmes aux limites de type elliptique. / Doctorat en sciences, Spécialisation mathématiques / info:eu-repo/semantics/nonPublished

Cette thèse est dédiée à l'analyse mathématique de quelques problèmes variationnels motivés par le modèle de Ginzburg-Landau en théorie de la supraconductivité. Dans la première partie on étudie l'existence de solutions pour les équations de Ginzburg-Landau sans champ magnétique et avec données au bord de type semi-rigides. Ces données consistent à prescrire le module de la fonction sur le bord du domaine ainsi que son degré topologique. C'est un cas particulier de problèmes à bord libre, ou la donnée complète de la fonction sur le bord est une inconnue du problème. L'existence de solutions à ce problème n'est pas assurée. En effet la méthode directe du calcul des variations ne peut pas s'appliquer car le degré sur le bord n'est pas continu pour la convergence faible dans l'espace de Sobolev adapté. On dit que c'est un problème sans compacité. En étudiant le phénomène de "bubbling" qui apparaît dans l'étude de tels problèmes on donne des résultats d'existence et de non existence de solutions. Dans le Chapitre 1 on étudie des conditions qui permettent d'affirmer que la différence entre deux niveaux d'énergie est strictement optimale. Pour cela on adapte une technique due à Brezis-Coron. Ceci nous permet de redémontrer un résultat (précédemment obtenu par Berlaynd Rybalko et Dos Santos) d'existence de solutions stables pour les équations de Ginzburg-Landau dans des domaines multiplement connexes. Dans le Chapitre 2 on considère les applications harmoniques a valeurs dans $R^2$ avec des conditions au bord de type degrés prescrits sur un anneau. On fait un lien entre ce problème et la théorie des surfaces minimales dans $R^3$ grâce à la différentielle quadratique de Hopf. Ceci nous conduit à l'étude des surfaces minimales bordées par deux cercles dans des plans parallèles. On prouve l'existence de telles surfaces qui ne sont pas des catenoides grâce a un résultat de bifurcation. On utilise alors les résultats obtenus pour déduire des théorèmes d'existence et de non existence de minimiseurs de l'énergie de Ginzburg-Landau à degrés prescrits dans un anneau. Dans ce troisième Chapitre on obtient des résultats pour une valeur du paramètre " grand. Le Chapitre 4 a pour objet l'étude des problèmes a degrés prescrits en dimension n3. On y montre la non existence des minimiseurs de la n-énergie de Ginzburg-Landau a degrés prescrits dans un domaine simplement connexe. On étudie ensuite des points critiques de type min-max pour une énergie perturbée. La deuxième partie est consacrée a l'analyse asymptotique des solutions des équations deGinzburg-Landau lorsque " tend vers zero. Sandier et Serfaty ont étudié le comportement asymptotique des mesures de vorticité associées aux équations. Ils ont notamment trouvé des conditions critiques sur les mesures limites dans le cas des équations avec et sans champ magnétique. Nous nous intéressons alors à ces conditions critiques dans le cas sans champ magnétique. Le problème de la régularité locale des mesures limites se ramène ainsi a l'étude de la régularité des fonctions stationnaires harmoniques dont le Laplacien est une mesure. Nous montrons que localement de telles mesures sont supportées par une union de lignes appartenant à l'ensemble des zéros d'une fonction harmonique / This thesis is devoted to the mathematical analysis of some variational problems. These problem sare motivated by the Ginzburg-Landau model related to the super conductivity. In the first part we study existence of solutions of the Ginzburg-Landau equations without magnetic eld but with semi-sti boundary conditions. These conditions are obtained by prescribing the modulus of the function on the boundary of the domain along with its topological degree. This is a particular case of free boundary problems, where the function on the boundary is an unknown of the problem. Existence of solutions of that problem does not necessary hold. Indeed we can not apply the direct method of the calculus of variations since the degree on the boundaryis not continuous with respect to the weak convergence in an appropriated Sobolev space. This is problem with loss of compactness. By studying the bublling" phenomenon which come upin such problems we obtain some existence and non existence results .In Chapter 1 we study conditions under which the dierence between two energy levels is strictly optimal. In order to do that we adapt a technique due to Brezis-Coron. This allow us to recover known existence results (previously obtained by Berlyand and Rybalko and DosSantos) for stable solutions of the Ginzburg-Landau equations in multiply connected domains. In Chapter 2 we are interested in harmonic maps with values in $R^2$ with prescribed degree boundary condition in an annulus. We make a link between this problem and the minimal surface theory in $R^3$ thanks to the so-called Hopf quadratic differential. This leads us to study immersed minimal surfaces bounded by two circles in parallel planes. We prove the existence of such surfaces die rent from catenoids by using a bifurcation argument. We then apply the results obtained to deduce existence and non existence results for minimizers of the Ginzburg-Landau energy with prescribed degrees. This is done in Chapter 3 where the results are obtained for large ".Chapter 4 is devoted to prescribed degree problems in dimension n3 . We prove the non existence of minimizers of the Ginzburg-Landau energy in simply connected domains. We then study min-max critical points of a perturbed energy. The second part is devoted to the asymptotic analysis of solutions of the Ginzburg-Landau equations when "goes to zero. Sandier and Serfaty studied the asymptotic behavior of the vorticity measures associated to these equations. They derived critical conditions on the limiting measures both with and without magnetic Field. We are interested by these conditions when there is no magnetic Field. The problem of the local regularity of the limiting measures is then equivalent to the study of regularity of stationary harmonic functions whose Laplacianis a measure. We show that locally such measures are concentrated on a union of lines which belong to the zero set of an harmonic function

Une modélisation numérique des ondes poroélastiques, décrites par le modèle de Biot, est proposée dans le domaine temporel. La dissipation visqueuse à l'intérieur des pores est décrite par le modèle de perméabilité dynamique de Johnson-Koplik-Dashen (JKD). Certains coefficients du modèle de Biot-JKD sont proportionnels à la racine carrée de la fréquence, introduisant dans le domaine temporel des dérivées fractionnaires décalées d'ordre 1/2, revenant à un produit de convolution. Basé sur une représentation diffusive, le produit de convolution est remplacé par un nombre fini de variables de mémoire satisfaisant une équation différentielle ordinaire locale en temps, menant au modèle de Biot-DA (diffusive approximation). Les propriétés des deux modèles sont analysées : hyperbolicité, décroissance de l'énergie, dispersion. On montre que la meilleure méthode de détermination des coefficients de l'approximation diffusive - quadratures de Gauss, optimisation linéaire ou non-linéaire sur la plage de fréquence d'intérêt - est l'optimisation non-linéaire. Une méthode de splitting est utilisée numériquement : la partie propagative est discrétisée par un schéma aux différences finies ADER d'ordre 4, et la partie diffusive est intégrée exactement. Les conditions de saut aux interfaces sont discrétisées avec une méthode d'interface immergée. Des simulations numériques sont présentées pour des milieux isotropes et isotropes transverses. Des comparaisons avec des solutions analytiques montrent l'efficacité et la précision de cette approche. Des simulations numériques en milieux complexes sont réalisées : influence de la porosité d'os spongieux, diffusion multiple en milieu aléatoire. / A time-domain numerical modeling of Biot poroelastic waves is proposed. The viscous dissipation in the pores is described using the dynamic permeability model of Johnson-Koplik-Dashen (JKD). Some of the coefficients in the Biot-JKD model are proportional to the square root of the frequency: in the time-domain, these coefficients introduce shifted fractional derivatives of order 1/2, involving a convolution product. Based on a diffusive representation, the convolution product is replaced by a finite number of memory variables that satisfy local-in-time ordinary differential equations, resulting in the Biot-DA (diffusive approximation). The properties of the two models are analyzed: hyperbolicity, decrease of energy, dispersion. To determine the coefficients of the diffusive approximation, different methods of quadrature are analyzed: Gaussian quadratures, linear or nonlinear optimization procedures in the frequency range of interest. The nonlinear optimization is shown to be the best way of determination. A splitting strategy is applied numerically: the propagative part is discretized using a fourth-order ADER scheme on a Cartesian grid, and the diffusive part is solved exactly. An immersed interface method is implemented to discretize the jump conditions at interfaces. Numerical experiments are presented for isotropic and transversely isotropic media. Comparisons with analytical solutions show the efficiency and the accuracy of this approach. Some numerical experiments are performed in complex media: influence of the porosity of a cancellous bone, multiple scattering across a set of random scatterers.

Cette thèse est consacrée à l'étude théorique de quelques modèles d'évolution non-linéaires issus de la mécanique des fluides. Nous distinguons trois parties indépendantes. La première partie de la thèse traite essentiellement de l'existence des poches de tourbillon en rotation uniforme (appelées aussi V-states) pour un modèle quasi-géostrophique non visqueux. Notre étude est répartie sur deux chapitres où les poches présentent des structures topologiques différentes. Dans le premier chapitre nous étudions le cas simplement connexe et nous validons l'existence de ces structures dans un voisinage du tourbillon de Rankine en utilisant des techniques de bifurcation. Dans le deuxième chapitre nous abordons le cas doublement connexe où la poche admet un seul trou. Plus précisément, proche d'un anneau donné, nous décrivons cette famille par des branches dénombrables bifurquant de cet anneau à certaines valeurs explicites des vitesses angulaires liées aux fonctions de Bessel. Notre étude théorique a été complétée par des simulations numériques portant sur les V-states limites et un bon nombre de constatations ont été formulées ouvrant la porte à de nouvelles perspectives de recherche. La seconde partie concerne l'étude du problème de Cauchy pour le système de Boussinesq non visqueux 2D avec des données initiales de type Yudovich. Le problème est dans un certain sens critique à cause de quelques termes comportant la transformée de Riesz dans la formulation tourbillon-densité. Nous donnons une réponse positive pour une sous-classe comprenant les poches de tourbillon régulières et singulières. Dans la dernière partie nous analysons le problème de la limite incompressible pour les équations d'Euler isentropiques 2D associées à des données initiales très mal préparées et pour lesquelles les tourbillons ne sont pas forcément bornés mais appartiennent plutôt à des espaces de type ''BMO'' à poids. On utilise principalement deux ingrédients: d'un côté les estimations de Strichartz pour contrôler la partie acoustique. D'un autre côté, on se sert de la structure de transport compressible du tourbillon et on démontre une estimation de propagation linéaire dans l'esprit d'un travail récent de Bernicot et Keraani mené dans le cas incompressible. / In this dissertation, we are concerned with the study of some non-linear evolution models arising in fluid mechanics. We distinguish three independent parts. The first part of the thesis deals with the existence of the rotating vortex patches (called also V-states) for an inviscid quasi-geostrophic model. Our study is divided into two chapters dealing with different topological structures of the V-states. In the first chapter we study the simply connected case and we prove the existence of such structures in a neighborhood of the Rankine vortices by using the bifurcation theory. In the second chapter we discuss the doubly connected case where the patches admit only one hole. More precisely, close to a given annulus we describe this family by countable branches bifurcating from this annulus at some explicit angular velocities related to Bessel functions of the first kind. Our theoretical study was completed by numerical simulations on the limiting V-states and a number of interesting numerical observation were formulated opening new research perspectives. The second part of the thesis concerns the local well-posedness theory for the inviscid Boussinesq system with rough initial data. The problem is in some sense critical due to some terms involving Riesz transforms in the vorticity-density formulation. We give a positive answer for a special sub-class of Yudovich data including smooth and singular vortex patches. In the last part we address the problem of the incompressible limit for the 2D isentropic fluids associated to ill-prepared initial data and for which the vortices are not necessarily bounded and belong to some weighted BMO spaces. We mainly use two ingredients: On one hand, the Strichartz estimates to control the acoustic part and prove that it does not contribute for low Mach number. On the other hand, we use the transport compressible structure of the vorticity and we establish a linear propagation estimate in the spirit of a recent work of Bernicot and Keraani conducted in the incompressible case. The first part of the thesis deals with the existence of the rotating vortex patches (called also V-states) for an inviscid quasi-geostrophic model. Our study is divided into two chapters dealing with different topological structures of the V-states. In the first chapter we study the simply connected case and we prove the existence of such structures in a neighborhood of the Rankine vortices by using the bifurcation theory. In the second chapter we discuss the doubly connected case where the patches admit only one hole. More precisely, close to a given annulus we describe this family by countable branches bifurcating from this annulus at some explicit angular velocities related to Bessel functions of the first kind. Our theoretical study was completed by numerical simulations on the limiting V-states and a number of interesting numerical observation were formulated opening new research perspectives. The second part of the thesis concerns the local well-posedness theory for the inviscid Boussinesq system with rough initial data. The problem is in some sense critical due to some terms involving Riesz transforms in the vorticity-density formulation. We give a positive answer for a special sub-class of Yudovich data including smooth and singular vortex patches. In the last part we address the problem of the incompressible limit for the 2D isentropic fluids associated to ill-prepared initial data and for which the vortices are not necessarily bounded and belong to some weighted BMO spaces. We mainly use two ingredients: On one hand, the Strichartz estimates to control the acoustic part and prove that it does not contribute for low Mach number. On the other hand, we use the transport compressible structure of the vorticity and we establish a linear propagation estimate in the spirit of a recent work of Bernicot and Keraani conducted in the incompressible case.

Les travaux exposés dans cette thèse sont consacrés à l’étude de méthodesprécises pour approcher des équations différentielles stochastiques ou deséquations aux dérivées partielles (EDP) déterministes. La première parties’inscrit dans le cadre du développement de méthodes visant à corriger le biaisdans les processus de diffusion paramétrique. Trois modèles sont étudiés enparticulier : Ornstein-Uhlenbeck, Auto-régressif et Moyenne mobile. A l’issuede ce travail, plusieurs approximations de biais ont été proposées suivant deuxapproches : la première consiste en un développement de Taylor del’estimateur obtenu alors que la seconde s'appuie sur une expansionstochastique de celui-ci.La deuxième partie de cette thèse porte sur l’approximation de l’équation de lachaleur obtenue après changement de variables à partir du modèle de Black etScholes. En général, on préfère utiliser des méthodes implicites pour résoudredes EDP paraboliques mais depuis quelques années, les méthodes dites deRunge-Kutta explicites stabilisées, sont de plus en plus utilisées. Nousmontrons que l’utilisation de ce type de méthodes explicites et notamment lesschémas ROCK donnent de très bons résultats même si les conditions initialessont peu régulières, ce qui est le cas dans les modèles financiers / The work presented in this thesis is devoted to the study of precise methods forapproximating stochastic differential equations (SDE) or deterministic partialdifferential equations (PDE). The first part is devoted to the development ofbias correction methods in parametric diffusion processes. Three models arestudied in particular : Ornstein-Uhlenbeck, auto-regressive and Movingaverage. At the end of this work, several approximations of bias have beenproposed following two approaches : the first consists in a Taylor developmentof the obtained estimator while the second one relies on a stochastic expansionof the latter.The second part of this thesis deals with the approximation of the heatequation obtained after changing variables from the Black-Scholes model. Likethe vast majority of PDE, this equation does not have an exact solution, sosolutions must be approached using explicit or implicit time schemes. Itis often customary to prefer the use of implicit methods to solve parabolic PDEsuch as the heat equation, but in the past few years, the stabilized explicitRunge-Kutta methods which have the largest possible domains of stabilityalong the negative real axis, are increasingly used. We show that the useof this type of explicit methods and in particular the ROCK (Runge-Orthogonal-Chebyshev-Kutta) schemes give very good results even if the initial conditionsare not very regular, which is the case in the financial models