In this paper, the problem on formation and construction of a shock wave for three dimensional compressible Euler equations with the small perturbed spherical initial data is studied. If the given smooth initial data satisfies certain nondegenerate condition, then from the results in [20], we know that there exists a unique blowup point at the blowup time such that the first order derivates of smooth solution blow up meanwhile the solution itself is still continuous at the blowup point. From the blowup point, we construct a weak entropy solution which is not uniformly Lipschitz continuous on two sides of shock curve, moreover the strength of the constructed shock is zero at the blowup point and then gradually increases. Additionally, some detailed and precise estimates on the solution are obtained in the neighbourhood of the blowup point.

In the last half a century, Computational Fluid Dynamics (CFD) has been established as an important complementary part and some times a significant alternative to Experimental and Theoretical Fluid Dynamics. Development of efficient computational algorithms for digital simulation of fluid flows has been an ongoing research effort in CFD. An accurate numerical simulation of compressible Euler equations, which are the gov-erning equations of high speed flows, is important in many engineering applications like designing of aerospace vehicles and their components. Due to nonlinear nature of governing equations, such flows admit solutions involving discontinuities like shock waves and contact discontinuities. Hence, it is nontrivial to capture all these essential features of the flows numerically. There are various numerical methods available in the literature, the popular ones among them being the Finite Volume Method (FVM), Finite Difference Method (FDM), Finite Element Method (FEM) and Spectral method. Kinetic theory based algorithms for solving Euler equations are quite popular in finite volume framework due to their ability to connect Boltzmann equation with Euler equations. In kinetic framework, instead of dealing directly with nonlinear partial differential equations one needs to deal with a simple linear partial differential equation. Recently, FEM has emerged as a significant alternative to FVM because it can handle complex geometries with ease and unlike in FVM, achieving higher order accuracy is easier. High speed flows governed by compressible Euler equations are hyperbolic partial differential equations which are characterized by preferred directions for information propagation. Such flows can not be solved using traditional FEM methods and hence, stabilized methods are typically introduced. Various stabilized finite element methods are available in the literature like Streamlined-Upwind Petrov Galerkin (SUPG) method, Galerkin-Least Squares (GLS) method, Taylor-Galerkin method, Characteristic Galerkin method and Discontinuous Galerkin Method. In this thesis a novel stabilized finite element method called as Kinetic Streamlined-Upwind Petrov Galerkin (KSUPG) method is formulated. Both explicit and implicit versions of KSUPG scheme are presented. Spectral stability analysis is done for explicit KSUPG scheme to obtain the stable time step. The advantage of proposed scheme is, unlike in SUPG scheme, diffusion vectors are obtained directly from weak KSUPG formulation. The expression for intrinsic time scale is directly obtained in KSUPG framework. The accuracy and robustness of the proposed scheme is demonstrated by solving various test cases for hyperbolic partial differential equations like Euler equations and inviscid Burgers equation. In the KSUPG scheme, diffusion terms involve computationally expensive error and exponential functions. To decrease the computational cost, two variants of KSUPG scheme, namely, Peculiar Velocity based KSUPG (PV-KSUPG) scheme and Circular distribution based KSUPG (C-KSUPG) scheme are formulated. The PV-KSUPG scheme is based on peculiar velocity based splitting which, upon taking moments, recovers a convection-pressure splitting type algorithm at the macroscopic level. Both explicit and implicit versions of PV-KSUPG scheme are presented. Unlike KSUPG and PV-KUPG schemes where Maxwellian distribution function is used, the C-KUSPG scheme uses a simpler circular distribution function instead of a Maxwellian distribution function. Apart from being computationally less expensive it is less diffusive than KSUPG scheme.

This work develops a discontinuous Galerkin finite element discretization of non- linear hyperbolic conservation equations with efficient and robust high order stabilization built on an entropy-based artificial viscosity approximation. The solutions of equations are represented by elementwise polynomials of an arbitrary degree p > 0 which are continuous within each element but discontinuous on the boundaries. The discretization of equations in time is done by means of high order explicit Runge-Kutta methods identified with respective Butcher tableaux. To stabilize a numerical solution in the vicinity of shock waves and simultaneously preserve the smooth parts from smearing, we add some reasonable amount of artificial viscosity in accordance with the physical principle of entropy production in the interior of shock waves. The viscosity coefficient is proportional to the local size of the residual of an entropy equation and is bounded from above by the first-order artificial viscosity defined by a local wave speed. Since the residual of an entropy equation is supposed to be vanishingly small in smooth regions (of the order of the Local Truncation Error) and arbitrarily large in shocks, the entropy viscosity is almost zero everywhere except the shocks, where it reaches the first-order upper bound. One- and two-dimensional benchmark test cases are presented for nonlinear hyperbolic scalar conservation laws and the system of compressible Euler equations. These tests demonstrate the satisfactory stability properties of the method and optimal convergence rates as well. All numerical solutions to the test problems agree well with the reference solutions found in the literature. We conclude that the new method developed in the present work is a valuable alternative to currently existing techniques of viscous stabilization.