Flood propagation over urban areas can cause an interaction between the free-surface flow and large underground pipe networks used for storm drainage and sewage, causing outflows and inflows at the bed. The associated waves may collide with each other and the surface waves. In this thesis the shallow water equations are used to model this type of wave interaction over dry or wet beds with bathymetry gradients and friction terms. The proposed shallow water scheme is solved based on finite volume high-resolution Godunov-type methods. The solver is well-balanced and can accurately balance the source terms and flux-gradients for the steady-state solutions. The solver also utilises a new type of Riemann wave speed to provide depth-positive results over nearly dry beds and dry states. Additionally a new type of source term is introduced in the continuity equation to model pipe inflow and outflow conditions at bed connections. For the standard one-dimensional shallow water equations the numerical results are validated with analytical solutions or other reference solutions provided in the literature. This includes the incipient Riemann problems for nearly dry and dry-states, steady flow over a hump in a rectangular channel and the wave propagation problem. Eventually, the generation of dry bed in the middle, over discontinuous topography is considered. Close agreement is achieved between the shallow water scheme and analytical or reference solutions for the above test cases. For the shallow water problems with influx/efflux source terms comparisons are made with STAR-CD, a commercial Navier-Stokes solver for general fluid flow prediction. The shallow water model is first used to simulate vertical flows through finite gaps in the bed. Next, the interaction of the vertical flows with a dam-break flow is considered for both dry and wet beds. An efflux number, En, is defined based on the vertical efflux velocity and the gap length. A parameter study is undertaken to investigate the effect of the one-dimensional approximation of the present model, for a range of non-dimensional efflux numbers. It is found that the shallow flow model gives sensible predictions at all times provided En<0.5, and for long durations for En>0.5. Dam break flow over an underground connecting pipe is also considered for the one-dimensional efflux problems. To solve two-dimensional problems the shallow water scheme uses the dimensional-splitting method which solves each one-dimensional Riemann problem in the x- and y-directions separately. The cross-derivative terms for second-order accuracy are incorporated by solving another Riemann problem in the orthogonal direction. For two-dimensional problems first the dam-break problems are considered over wet and dry beds. Then, flood propagation over complex terrain is demonstrated. Next, efflux discharge is modelled in isolation over a dry bed and then with dam-break interaction, comparing with STAR-CD results. Again very good agreement is shown between the two-dimensional shallow water model and STAR-CD for the efflux numbers of En<0.5. For modelling the inundation problem over an underground pipe network the solver is coupled with the general underground pipe network solver to calculate the efflux discharge as the flood waves pass through the pipe network. For analysing the pipe network with unknown effluxes an additional set of equations is incorporated into the solution of a general pipe network solver. The shallow water solver coupled to an underground pipe network is then used to simulate dam-break interaction with pipe networks with 9 and 25 nodes to demonstrate the versatility of the method.

This thesis describes the development, verification, and validation of a three dimensional time domain thermoacoustic solver. The purpose of the solver is to predict the frequencies, modeshapes, linear growth rates, and limit cycle amplitudes for combustion instability modes in gas turbine combustion chambers. The linearized Euler equations with nonlinear heat release source terms are solved using the finite volume method. The treatment of mean density gradients was found to be vital to the success of frequency and modeshape predictions due to the sharp density gradients that occur across deflagration waves. In order to treat mean density gradients with physical fidelity, a non-conservative finite volume method based on the wave propagation approach to the Riemann problem is applied. For modelling unsteady heat release, user input flexibility is maximized using a virtual class hierarchy within the OpenFOAM C++ library. Unsteady heat release based on time lag models are demonstrated. The solver gives accurate solutions compared with analytical methods for one-dimensional cases involving mean density gradients, cross-sectional area changes, uniform mean flow, arbitrary impedance boundary conditions, and unsteady heat release in a one-dimensional Rijke tube. The solver predicted resonant frequencies within 1% of the analytical solution for these verification cases, with the dominant component of the error coming from the finite time interval over which the simulation is performed. The linear iii growth rates predicted by the solver for the Rijke tube verification were within 5% of the theoretical values, provided that numerical dissipation effects were controlled. Finally, the solver is then used to predict the frequencies and limit cycle amplitudes for two lab scale experiments in which detailed acoustics data are available for comparison. For experiments at the University of Melbourne, an empirical flame describing function was provided. The present simulation code predicted a limit cycle of 0.21 times the mean pressure, which was in close agreement with the estimate of 0.25 from the experimental data. The experiments at Purdue University do not yet have an empirical flame model, so a general vortex-shedding model is proposed on physical grounds. It is shown that the coefficients of the model can be tuned to match the limit cycle amplitude of the 2L mode from the experiment with the same accuracy as the Melbourne case. The code did not predict the excitation of the 4L mode, therefore it is concluded that the vortex-shedding model is not sufficient and must be supplemented with additional heat release models to capture the entirety of the physics for this experiment.

Riemann problems are of prime importance in computational fluid dynamics simulations using finite elements or finite volumes discretizations. In some applications, billions of Riemann problems might need to be solved in a single simulation, therefore it is important to have reliable and computationally efficient algorithms to do so. Given the nonlinearity of the flux function in most systems considered in practice, to obtain an exact solution for the Riemann problem explicitly is often not possible, and iterative solvers are required. However, because of issues found with existing iterative solvers like lack of convergence and high computational cost, their use is avoided and approximate solvers are preferred. In this thesis work, motivated by the advances in computer hardware and algorithms in the last years, we revisit the possibility of using iterative solvers to compute the exact solution for Riemann problems. In particular, we focus on the development, implementation, and performance comparison of iterative Riemann solvers for the shallow water and Euler equations. In a one-dimensional homogeneous framework for these systems, we consider several initial guesses and iterative methods for the computation of the Riemann solution. We find that efficient and reliable iterative solvers can be obtained by using recent estimates on the Riemann solution to modify and combine well-known methods. Finally, we consider the application of these solvers in finite volume simulations using the wave propagation algorithms implemented in Clawpack.

An axisymmetric, Navier-Stokes finite volume flow solver, which uses Harten, Lax and van Leer (HLL) and Harten, Lax and van Leer&ndash / Contact (HLLC) upwind flux differencing scheme for spatial and uses Runge-Kutta explicit multi-stage time stepping scheme for temporal discretization on unstructured meshe is developed. Developed solver can solve the compressible axisymmetric flow. The spatial accuracy of the solver can be first or second order accurate. Second order accuracy is achieved by piecewise linear reconstruction. Gradients of flow variables required for piecewise linear reconstruction are calculated by Green-Gauss theorem. Baldwin-Lomax turbulent model is used to compute the turbulent viscosity. Approximate Riemann solver of HLL and HLLC implemented in solver are validated by solving a cylindrical explosion case. Also the solver&rsquo / s capability of solving unstructured, multi-zone domain is investigated by this problem. First and second order results of solver are compared by solving the flow over a circular bump. Axisymmetric flow in solid propellant rocket motor is solved in order to validate the axisymmetric feature of solver. Laminar flow over flat plate is solved for viscous terms validation. Turbulent model is studied in the flow over flat plate and flow with mass injection test cases.

L’objectif principal de cette thèse concerne l’étude, la conception et la mise en œuvre numérique de schémas volumes finis associés aux solveurs de type Godunov. On s’intéresse à des systèmes hyperboliques de lois de conservation non linéaires, avec une attention particulière sur les équations d’Euler et les équations MHD idéale. Tout d’abord, nous dérivons un solveur de Riemann simple et véritablement multidimensionnelle, pouvant s’appliquer à tout système de lois de conservation. Ce solveur peut être considéré comme une généralisation 2D de l’approche HLL. Les ingrédients de base de la dérivation sont : la consistance avec la formulation intégrale et une utilisation adéquate des relations de Rankine-Hugoniot. Au final nous obtenons des expressions assez simples et applicables dans les contextes des maillages structurés et non structurés. Dans un second temps, nous nous intéressons à la préservation, au niveau discret, de la contrainte de divergence nulle du champ magnétique pour les équations de la MHD idéale. Deux stratégies sont évaluées et nous montrons comment le solveur de Riemann multidimensionnelle peut être utilisé pour obtenir des simulations robustes à divergence numérique nulle. Deux autres points sont abordés dans cette thèse : la méthode de relaxation pour un système Euler-Poisson pour des écoulements gravitationnels en astrophysique, la formulation volumes finis en coordonnées curvilignes. Tout au long de la thèse, les choix numériques sont validés à travers de nombreux résultats numériques. / The main objective of this thesis concerns the study, design and numerical implementation of finite volume schemes based on the so-Called Godunov-Type solvers for hyperbolic systems of nonlinear conservation laws, with special attention given to the Euler equations and ideal MHD equations. First, we derive a simple and genuinely two-Dimensional Riemann solver for general conservation laws that can be regarded as an actual 2D generalization of the HLL approach, relying heavily on the consistency with the integral formulation and on the proper use of Rankine-Hugoniot relations to yield expressions that are simple enough to be applied in the structured and unstructured contexts. Then, a comparison between two methods aiming to numerically maintain the divergence constraint of the magnetic field for the ideal MHD equations is performed and we show how the 2D Riemann solver can be employed to obtain robust divergence-Free simulations. Next, we derive a relaxation scheme that incorporates gravity source terms derived from a potential into the hydrodynamic equations, an important problem in astrophysics, and finally, we review the design of finite volume approximations in curvilinear coordinates, providing a fresher view on an alternative discretization approach. Throughout this thesis, numerous numerical results are shown.

In the thesis work, a code to solve the two dimensional compressible Euler equations for external flows around arbitrary geometries have been developed. A Cartesianmesh generator is incorporated to the solver. Hence the pre-processing can be performed together with the solution within a single code. The code is written in the C++ programming language and its object oriented capabilities have been exploited to save memory in the data structure developed. The Cartesian mesh is formed by dividing squares successively into its four quadrants. The main advantage of using this type of a mesh is the ability to generate meshes around geometries of arbitrary complexity quickly and to adapt the mesh easily based on the solution. The main disadvantage of this method is that the treatment of the cells that are cut by the geometry. For the solution procedure Roe&rsquo / s method as well as flux vector splitting methods are used for the flux evaluation. The flux vector splitting schemes used are van Leer, AUSM, AUSMD and AUSMV methods. Time discretization is performed using a multi-stage method. To increase the accuracy least squares reconstruction is employed. The code is validated by performing calculations around a NACA0012 airfoil profile. The effect of reconstruction is demonstrated by plotting the pressure coefficient on the airfoil. The distribution obtained using reconstruction is very close to the experimental one while there is a considerable deviation for the case without reconstruction. Also the shock capturing capabilities of different methods have been investigated. In addition the performance of each method is analyzed for flow around an NLR 7301 airfoil with a flap.

The purpose of this thesis is to implement Finite Volume Weighted Essentially Non-Oscillatory (FV-WENO) scheme to solution of one and two-dimensional discretised Euler equations with different flux algorithms. The effects of the different fluxes on the solution have been tested and discussed. Beside, the effect of the grid on these fluxes has been investigated. Weighted Essentially Non-Oscillatory (WENO) schemes are high order accurate schemes designed for problems with piecewise smooth solutions that involve discontinuities. WENO schemes have been successfully used in applications, especially for problems containing both shocks and complicated smooth solution structures. Fluxes are used as building blocks in FV-WENO scheme. The efficiency of the scheme is dependent on the fluxes used in scheme The applications tested in this thesis are the 1-D Shock Tube Problem, Double Mach Reflection, Supersonic Channel Flow, and supersonic Staggered Wedge Cascade. The numerical solutions for 1-D Shock Tube Problem and the supersonic channel flow are compared with the analytical solutions. The results for the Double Mach Reflection and the supersonic staggered cascade are compared with results from literature.

No

Le présent travail est consacré à l’approximation numérique des solutions du modèle aux moments M1 pour le transfert radiatif. Il s’agit, ici, de développer des solveurs numériques performants et précis capables de prédire avec précision et robustesse des écoulements où le transfert radiatif joue un rôle essentiel. Dans ce sens, plusieurs méthodes numériques ont été envisagées pour la dérivation des schémas numériques de type solveur de Godunov. Une attention particulière a été portée sur les solveurs préservant les ondes de contact stationnaires. En particulier, un schéma de relaxation et un solveur HLLC sont présentés dans ce travail. Pour chacun de ces solveurs, la robustesse de la méthode a été établie (positivité de l’énergie radiative et limitation du flux radiatif). La validation et l’intérêt des méthodes abordées sont exhibés à travers de nombreuses expériences numériques mono et multidimensionelles. / The present work is dedicated to the numerical approximation of the M1 moments model solutions for radiative transfer. The objective is to develop efficient and accurate numerical solvers, able to provide with precise and robust computations of flows where radiative transfer effects are important. With this aim, several numerical methods have been considered in order to derive numerical schemes based on Godunov type solvers. A particular attention has been paid to solvers preserving the stationary contact waves. Namely, a relaxation scheme and a HLLC solver are presented in this thesis. The robustness of each of these solvers has been established (radiative energy positivity and radiative flux limitation). Several numerical experiments in one and two space dimensions validate the developed methods and outline their interest.

Ce sujet concerne un travail pluridisciplinaire mathématique et astrophysique. Le but de cette thèse est l'étude des modèles d'hydrodynamique radiative dont l'application est bien évidemment très vaste en physique et astrophysique. Les modèles M1-multigroupes sont explorés pour décrire le transfert radiatif sans faire à priori d'hypothèse sur la profondeur optique du milieu. L'intérêt qui découle directement de ce travail est le développement du code d'hydrodynamique radiative HADES 2D permettant le calcul massivement parallèle. Il autorise des simulations dans des configurations astrophysiques réalistes en termes de nombre de Mach et de contraste de densité et de température entre les différents milieux. Nous nous sommes concentrés sur deux applications intéressantes : les jets d'étoiles jeunes et les chocs radiatifs dont les premières simulations seront présentées. / This topic is a multidisciplinary work between mathematics and astrophysics. The aim of this thesis is the study of radiation hydrodynamic models of which application is obviously very broad in physics and astrophysics. M1-multigroup models are explored to describe the radiative transfer without a priori assumption on the optical depth of the medium. The interest ensuing directly from this work is the development of a radiation hydrodynamic code, namely HADES 2D, for massively parallel computing. It allows simulations in realistic astrophysical configurations in terms of Mach number, density and temperature contrasts between different environments. We focused on two interesting applications: the jets from young stars and the radiative shocks of which first simulations will be presented.