We investigate the impact of cash reserves upon the optimal behaviour of a modelled firm that has uncertain future revenues. To achieve this, we build up a corporate financing model of a firm from a Real Options foundation, with the option to close as a core business decision maintained throughout. We model the firm by employing an optimal stochastic control mathematical approach, which is based upon a partial differential equations perspective. In so doing, we are able to assess the incremental impacts upon the optimal operation of the cash constrained firm, by sequentially including: an optimal dividend distribution; optimal equity financing; and optimal debt financing (conducted in a novel equilibrium setting between firm and creditor). We present efficient numerical schemes to solve these models, which are generally built from the Projected Successive Over Relaxation (PSOR) method, and the Semi-Lagrangian approach. Using these numerical tools, and our gained economic insights, we then allow the firm the option to also expand the operation, so they may also take advantage of favourable economic conditions.

Cette thèse présente l'étude et le développement de méthodes numériques pour la résolution d'équations de transport, en particulier d'une méthode de remapping bidimensionnel dont un avantage important par rapport aux algorithmes existants est la propriété de conservation de la masse. De nombreux cas-tests permettront de comparer ces approches entre elles ainsi qu'à des méthodes de référence. On s'intéressera en particulier aux équations dites de Vlasov-Poisson et du Centre-Guide, qui apparaissent très classiquement dans le cadre de la physique des plasmas. / This PhD thesis presents the study and development of numerical methods for the resolution of transport equations, in particular a bidimensional remapping method whose main advantage over existing algorithms is the property of mass conservation. Numerous test cases are presented in order to compare these approaches with regard to the others and with reference methods. Focus is made on the so-called Vlasov-Poisson and Center-Guide equations, that appear very classically in the domain of plasma physics.

Cette thèse propose et analyse des méthodes numériques pour la résolution de l'équation de Vlasov. Cette équation modélise l'évolution d'une espèce de particules chargées sous l'effet d'un champ électromagnétique. La première partie est consacrée à une analyse mathématique de schémas semi-Lagrangiens résolvant l'équation de transport linéaire qui constituent la brique de base des méthodes de splitting directionnel.Des méthodes de résolution de l'équation de Vlasov couplée à l'équation de Poisson, dans le cas où uniquement le champ électrique est considéré, sont optimisées dans la seconde partie. Il s'agit d'optimisation en temps de calcul par l'utilisation de cartes graphiques (GPU) et l'utilisation d'un maillage non homogène.Dans la troisième et dernière partie, nous étudierons une méthode numérique de calcul de l'opérateur de gyromoyenne intervenant dans la théorie gyrocinétique que nous appliquerons à l'équation de quasi-neutralité. / This thesis proposes and analyzes numerical methods for solving the Vlasov equation. This equation models the evolution of a species of charged particles under the effet of an electromagnetic field. The first part is devoted to a mathematical analysis of semi-Lagrangian schemes solving the linear transport equation which is the basic building block of directional splitting methods.Solving methods for the Vlasov equation coupled to the Poisson equation, in the case where only the electric field is considered, are optimized in the second part. This optimization relates to the time of calculation by the use of Graphics Processing Unit (GPU) and the use of an inhomogeneous mesh.In the third and final part, we study a numerical method for calculating the gyroaverage operator involved in gyrokinetic theory. This method will be applied to solve the quasi-neutrality equation.

This work makes basic analysis of historical and current algorithms for smoke animation. Modern approaches to rendering volumetric data are briefly described. We choose algorithms for implementation on basis of this analysis. These algorithms are described in detail and we make emphasis on their important properties according to dedication of this work. Detailed description of implementation follows along with performance measurement. Conclusion evaluates results of work and proposes possible extensions.

Beaucoup de méthodes numériques ont été développées pour résoudre l'équation de Vlasov, car obtenir des simulations numériques précises en un temps raisonnable pour cette équation est un véritable défi. Cette équation décrit en effet l'évolution de la fonction de distribution de particules (électrons/ions) qui dépend de 3 variables d'espace, 3 variables de vitesse et du temps. L'idée principale de cette thèse est de réécrire l'équation de Vlasov sous forme d'un système hyperbolique par semi-discrétisation en vitesse. Cette semi-discrétisation est effectuée par méthode d'éléments finis. Le modèle ainsi obtenu est appelé équation de Vlasov réduite. Nous proposons différentes méthodes numériques pour résoudre efficacement ce modèle: méthodes des volumes finis, méthodes semi-Lagrangiennes et méthodes Galerkin discontinus. / Many numerical methods have been developed in order to selve the Vlasov equation, because computing precise simulations in a reasonable time is a real challenge. This equation describes the time evolution of the distribution function of charged particles (electrons/ions), which depends on 3 variables in space, 3 in velocity and time. The main idea of this thesis is to rewrite the Vlasov equation in the form of a hyperbolic system using a semi-discretization of the velocity. This semi-discretization is achieved using the finite element method. The resulting model is called the reduced Vlasov equation. We propose different numerical methods to salve this new model efficiently: finite volume methods, semi-Lagrangian methods and discontinuous Galerkin methods.

Cette thèse porte sur le comportement en temps long de solutions d’équations de type Vlasov, principalement le modèle Vlasov-HMF. On s’intéresse en particulier au phénomène d’amortissement Landau, prouvé mathématiquement dans divers cadres, pour plusieurs équations de type Vlasov, comme l’équation de Vlasov-Poisson ou le modèle Vlasov-HMF, et présentant certaines analogies avec le phénomène d’amortissement non visqueux pour l’équation d’Euler 2D. Les résultats qui y sont décrits sont les suivants. Le premier est un théorème d’amortissement Landau pour des solutions numériques du modèle Vlasov-HMF, obtenues par discrétisation en temps de ce dernier via des méthodes de splitting. Nous prouvons en outre la convergence des schémas numériques. Le second est un théorème d’amortissment Landau pour des solutions du modéle Vlasov-HMF linéarisé autour d’états stationnaires inhomogènes. Ce théorème est accompagné de nombreuses simulations numériques destinées à étudier numériquement le cas non-linéaire, et semblant mettre en lumière de nouveaux phénomènes. Enfin, le dernier résultat porte sur la discrétisation en temps de l’équation d’Euler 2D par un intégrateur de Crouch-Grossman symplectique. Nous prouvons la convergence du schéma. / This thesis concerns the long time behavior of certain Vlasov equations, mainly the Vlasov- HMF model. We are in particular interested in the celebrated phenomenon of Landau damp- ing, proved mathematically in various frameworks, foar several Vlasov equations, such as the Vlasov-Poisson equation or the Vlasov-HMF model, and exhibiting certain analogies with the inviscid damping phenomenon for the 2D Euler equation. The results described in the document are the following.The first one is a Landau damping theorem for numerical solutions of the Vlasov-HMF model, constructed by means of time-discretizations by splitting methods. We prove more- over the convergence of the schemes. The second result is a Landau damping theorem for solutions of the Vlasov-HMF model linearized around inhomogeneous stationary states. We provide moreover a quite large amount of numerical simulations, which are designed to study numerically the nonlinear case, and which seem to show new phenomenons. The last result is the convergence of a scheme that discretizes in time the 2D Euler equation by means of a symplectic Crouch-Grossmann integrator.