The class of hyperbolic conservation laws model the phenomena of non-linear wave propagation, including the presence and propagation of discontinuities and expansion waves. Such nonlinear systems can generate discontinuities in the so-lution even for smooth initial conditions. Presence of discontinuities results in break down of a solution in the classical sense and to show existence, weak for-mulation of a problem is required. Moreover, closed form solutions are di cult to obtain and in some cases such solutions are even unavailable. Thus, numerical algorithms play an important role in solving such systems. There are several dis-cretization techniques to solve hyperbolic systems numerically and Finite Volume Method (FVM) is one of such important frameworks. Numerical algorithms based on FVM are broadly classi ed into two categories, central discretization methods and upwind discretization methods. Various upwind and central discretization methods developed so far di er widely in terms of robustness, accuracy and ef-ciency and an ideal scheme with all these characteristics is yet to emerge. In this thesis, novel upwind and central schemes are formulated for various hyper-bolic systems, with the aim of maintaining right balance between accuracy and robustness. This thesis is divided into two parts. First part consists of the formulation of upwind methods to simulate genuine weakly hyperbolic (GWH) systems. Such systems do not possess full set of linearly independent (LI) eigenvectors and some of the examples include pressureless gas dynamics system, modi ed Burgers' sys-tem and further modi ed Burgers' system. The main challenge while formulating an upwind solver for GWH systems, using the concept of Flux Di erence Splitting (FDS), is to recover full set of LI eigenvectors, which is done through addition of generalized eigenvectors using the theory of Jordan Canonical Forms. Once the defective set of LI eigenvectors are completed, a novel (FDS-J) solver is for-mulated in such a manner that it is independent of generalized eigenvectors, as they are not unique. FDS-J solver is capable of capturing various shocks such as -shocks, 0-shocks and 00-shocks accurately. In this thesis, the FDS-J schemes are proposed for those GWH systems each of which have one particular repeated eigenvalue with arithmetic multiplicity (AM) greater than one. Moreover, each ux Jacobian matrix corresponding to such systems is similar to a unique Jordan matrix. After the successful treatment of genuine weakly hyperbolic systems, this strategy is further applied to those weakly hyperbolic subsystems which result on employ-ing various convection-pressure splittings to the Euler ux function. For example, Toro-Vazquez (TV) splitting and Zha-Bilgen (ZB) type splitting approaches to split the Euler ux function yield genuine weakly hyperbolic convective parts and strict hyperbolic pressure parts. Moreover, the ux Jacobian of each convective part is similar to a Jordan matrix with at least two lower order Jordan blocks. Based on the lines of FDS-J scheme, we develop two numerical schemes for Eu-ler equations using TV splitting and ZB type splitting. Both the new ZBS-FDS and TVS-FDS schemes are tested on various 1-D shock tube problems and out of two, contact capturing ZBS-FDS scheme is extended to 2-dimensional Euler system where it is tested successfully on various test cases including many shock instability problems. Second part of the thesis is associated with the development of simple, robust and accurate central solvers for systems of hyperbolic conservation laws. The idea of splitting schemes together with the notion of FDS is not easily extendable to systems such as shallow water equations. Thus, a novel central solver Convection Isolated Discontinuity Recognizing Algorithm (CIDRA) is formulated for shallow water equations. As the name suggests, the convective ux is isolated from the total ux in such a way that other ux, in present case other ux represents celerity part, must possess non-zero eigenvalue contribution. FVM framework is applied to each part separately and ux equivalence principle is used to x the coe cient of numerical di usion. CIDRA for SWE is computed on various 1-D and 2-D benchmark problems and extended to Euler systems e ortlessly. As a further improvement, a scalar di usion based algorithm CIDRA-1 is designed for v Euler systems. The scalar di usion coe cient depends on that particular part of the Rankine-Hugoniot (R-H) condition which involves total energy of the system as a direct contribution. This algorithm is applied to a variety of shock tube test cases including a class of low density ow problems and also to various 2-D test problems successfully. vi

Dans cette thèse, on étudie des systèmes de particules en interaction dont le comportement est lié à certaines équations aux dérivées partielles lorsque le nombre de particules tend vers l’infini. L’équation de Kuramoto-Sivashinsky modélise par exemple la propagation de certains fronts de flamme, la topographie de la surface d’une couche mince en cours de croissance, et fait apparaître des structures macroscopiques. Un modèle de particules en interaction par un couplage harmonique des vitesses, attractif aux premières vitesses voisines, répulsive aux secondes voisines, associée à des collisions élastiques, produit des profils de vitesses analogues aux fronts de flamme. On observe également la création et l’annihilation d’agrégats de particules. Un autre modèle, où les particules fusionnent lors des collisions en préservant masse et quantité de mouvement, et avec uniquement attraction au plus proche voisin, permet de retrouver un modèle de type gaz sans pression avec viscosité. Ces modèles sont étudiés théoriquement, en particulier les facteurs de mise à l’échelle des forces d’interaction sont précisés pour obtenir les équations correctes dans la limite du grand nombre de particules. Des simulations numériques confirment la validité et la pertinence des modèles. / This work is concerned by systems of interacting particles, which are linked to partial derivative equations when the particle number becomes large enough. The Kuramoto-Sivashinsky equation is actually modeling as well the front flame propagation as the morphology of growing interfaces, in deposition, for example. Moreover, surface periodical macroscopic structuring is occurring. An interacting particle model through an harmonic velocity coupling, attractive with the first velocity-neighbor and repulsive for the second neighbors, associated with elestic collisions. This model thus provides us with velocity profiles close to those of front flame propagation. Creation and annihilation of particle clusters is also observed. Another model, where particle are merging during collisions, while retaining mass and momentum conservation and with only nearest neighbor attraction, allows to recover a viscous pressureless gas model. These models are studied using mathematical tools. Especially interaction scaling factors are determined for obtaining the suitable equations in the large particle number limit. The numerical simulations confirm the relevance of the models.

<p>In the first part, the pressureless gas dynamic system with source (body force) is examined and solved by using Shadow Waves. The source represents gravity and Shadow Wave solution (containing the delta function) shows acceleration (contrary to shocks, for example). In the second part, one will nd numerical calculations that conrms the above results.</p> / <p>Rad je posvecen analizi modela gasa bez pritiska uz dodatak izvora. Model je resen koriscenjem senka talasa. U ovom slucaju, izvor predstavlja uticaj gravitacije na cestice u modelu. Za razliku od udarnih talasa, talasi senke koje sadrze delta funkciju, krecu se ubrzano pod gravitacionim uticajem. U drugom delu rada su naprevljeni numericki eksperimenti koji potvrdjuju teoijske rezultate.</p>

The Eulerian description of dispersed two-phase flows results in a system of partial differential equations describing characteristics of the flow, namely volume fraction, density and velocity of the two phases, around any point in space over time. When pressure forces are neglected or a same pressure is considered for both phases, the resulting system is weakly hyperbolic and solutions may exhibit vacuum states (regions void of the dispersed phase) or localized unbounded singularities (delta shocks) that are not physically desirable. Therefore, it is crucial to find a physical way for preventing the formation of such undesirable solutions in weakly hyperbolic Eulerian two-phase flow models. This thesis focuses on the mathematical analysis of an Eulerian model for air- droplet flows, here called the Eulerian droplet model. This model can be seen as the sticky particle system with a source term and is successfully used for the prediction of droplet impingement and more recently for the prediction of particle flows in air- ways. However, this model includes only one-way momentum exchange coupling, and develops delta shocks and vacuum states. The main goal of this thesis is to improve this model, especially for the prevention of delta shocks and vacuum states, and the adjunction of two-way momentum exchange coupling. Using a characteristic analysis, the condition for loss of regularity of smooth solutions of the inviscid Burgers equation with a source term is established. The same condition applies to the droplet model. The Riemann problems associated, respectively, to the Burgers equation with a source term and the droplet model are solved. The characteristics are curves that tend asymptotically to straight lines. The existence of an entropic solution to the generalized Rankine-Hugoniot conditions is proven. Next, a way for preventing the formation of delta shocks and vacuum states in the model is identified and a new Eulerian droplet model is proposed. A new hierarchy of two-way coupling Eulerian models is derived. Each model is analyzed and numerical comparisons of the models are carried out. Finally, 2D computations of air-particle flows comparing the new Eulerian droplet model with the standard Eulerian droplet model are presented.