Browsing by Author "Toro, E. F."
Now showing 1 - 11 of 11
Results Per Page
Sort Options
Item Open Access An adaptive gridding technique for conservation laws on complex domains(Cranfield University, 1997-05) Boden, E. P.; Toro, E. F.Obtaining accurate solutions to flows that involve discontinuous features still re- mains one of the most difficult tasks in computational fluid dynamics today. Some discontinuous features, such as shear waves and material interfaces, are quite deli- cate, yet they have a profound effect on the rest of the flow field. The accuracy of the numerical scheme and the quality of the grid discretisation of the flow domain, are both critical when computing multi-dimensional discontinuous solutions. Here, the second order WAF scheme is used in conjuction with an adaptive grid algorithm, which is able to automatically modify the grid in regions of discontinuous features and solid boundaries. The grid algorithm is a combination of two successful ap- proaches, namely Chimera and Cartesian grid Adaptive Mesh Refinement (AMR). The Chimera approach is able to accurately represent non-Cartesian boundaries, whilst the AMR approach yields significant savings in memory storage and cPu time. The combined algorithm has been thoroughly validated for convection test problems in gas dynamics. The computed solutions compare well with other numerical and experimental results. These tests have also been used to assess the efficiency of the grid adaption algorithms. Finally, the approach is applied to axi-symmetric, two- dimensional, two-phase, reactive flows in the context of internal ballistics problems. Again, the computed results are compared with other numerical and experimental results.Item Open Access Fully discrete arbitrary-order schemes for a model hyperbolic conservation law(1993) Shi, Jian; Toro, E. F.We investigate the fully discrete methodology and establish a formula from which two-level explicit fully discrete arbitrary-order (both in space and time) conservative numerical schemes for a model hyperbolic conservation law can be derived. To illustrate this approach fully discrete second, third and fourth order numerical schemes are presented.Item Open Access Fully discrete high resolution schemes for systems of conservation laws(Cranfield University, 1994-09) Shi, Jian; Toro, E. F.Effective and robust high resolution schemes are of vital importance for simulation of viscous and inviscid flows. Since second-order high resolution schemes in practice are inadquate for many applications, large efforts have been put towards developing higher- order accurate schemes in the past. Although some progress has been made, the efforts were frustrated by the lack of effective and robust new schemes. Therefore this thesis is aimed at challenging this difficult but very important issue. Some new theories and methodologies were established during this research, which covers the linear stability analysis for high-order numerical schemes; the fully discrete techniques for model equations; the formulation of conservative high-order schemes and the high-order Total Variation Diminishing (TVD) schemes. According to these theories arbitrary-order high resolution schemes can be developed. To illustrate the methodologies second-, third-, fourth-, and 20th-order schemes are presented. These high resolution schemes were tested and validated by solving some popular test problems for one and two dimensional Euler and incompressible Navier-Stokes equations. The efficiency and robustness are the features of these high-order schemes.Item Open Access Fully discrete high-order TVD schemes for a scalar hyperbolic conservation law(1993) Shi, Jian; Toro, E. F.In this paper we investigate fully discrete high-order TVD schemes for a scalar hyper- bolic conservation law using flux limiters . Formulae which define Courant number dependent TVD regions for second and third-order TVD schemes are established. A semi-empirical TVD procedure for an m-th order scheme (m ≥ 4) are proposed and tested.Item Open Access A linearised Reimann solver for Gudonov-type methods(1991) Toro, E. F.Very simple linearisations for the solution to the Riemann problem for the time-dependent and for the steady supersonic Euler equations are presented. When used locally in conjunction with Godunov-type methods, computing savings by a factor of about four, relative to the use of exact Riemann solvers, can be achieved. For severe flow regimes however, the linearisation looses accuracy and robustness. We then propose the use of a Riemann-solver adaptation procedure. This retains the accuracy and robustness of the exact Riemann solver and the computational efficiency of the cheap linearised Riemann solver. Also, reliable and simple switching criteria are presented. Numerical results for one, two and three-dimensional test problems suggest that the resulting numerical methods are competitive for practical applications, in terms of robustness, accuracy and computational efficiency.Item Open Access A linearised Reimann solver for the time-dependent Euler equations of gas dynamics(1991) Toro, E. F.The time-dependent Euler equations of Gas Dynamics are a set of non-linear hyperbolic conservation laws that admit discontinuous solutions (e.g. shocks). In this paper we are concerned with Riemann-problem based numerical methods for solving the general initial-value problem for these equations. We present an approximate, linearised Riemann solver for the time-dependent Euler equations. The solution is direct and involves few and simple arithmetic operations. The Riemann solver is then used, locally, in conjunction with the WAF numerical method to solve the time-dependent Euler equations in one and two space dimensions with general initial data. For flows with shocks waves of moderate strength the computed results are very accurate. For severe flow regimes we advocate the use of the present linearised Riemann solver in combination with the exact Riemann solver in an adaptive fashion. Numerical experiments demonstrate that such an approach can be very successful. One and two-dimensional test problems show that the linearised Riemann solver is used in over 99% of the flow field producing net computing savings by a factor of about 2. A reliable and simple switching criterion is also presented. Results show that the adaptive approach effectively provides the resolution and robustness of the exact Riemann solver at the computing cost of the simple linearised Riemann solver. The relevance of the present methods concerns the numerical solution of multi-dimensional problems accurately and economically.Item Open Access A numerical investigation of the artificial compressibility method for the solution of the Navier-Stokes equations(1992) Elsworth, D. T.; Toro, E. F.The Artificial Compressibility approach is an important numerical method for solving the incompressible Navier-Stokes Equations. The application of high resolution numerical methods to the equations of the artificial compressibility approach is a relatively new phenomenon and deserves further investigation. In this paper we examine the performance of five Riemann solvers: an exact Riemann solver, and four approximate solvers. The application of reflective boundary conditions is investigated, as well as the way in which the artificial compressibility coefficient is chosen.Item Open Access Riemann problems and the WAF method for two-dimensional shallow water equations(Cranfield Institute of Technology; College of Aeronautics, 1990) Toro, E. F.Item Open Access Some aspects of shock capturing methods for gas dynamics(1991) Toro, E. F.Three topics on modern shock capturing methods for the time-dependent Euler equations of Gas Dynamics are addressed. First we present the Weighted Average Flux Method (or WAF), one of several Riemann-problem based shock capturing methods. Then we deal with the Riemann problem. We present an efficient exact Riemann solver, a robust non-iterative Riemann solver based on the behaviour of the exact solver, and an improved version of the Harten-Lax-van Leer Riemann solver. Also, a very simple linearised Riemann solver is presented together with a Riemann solver adaptation procedure. We also present a Riemann-solver adaptation procedure that has proved successful. Applications of the WAF method with the various Riemann solvers are presented.Item Open Access TVD regions for the weighted average flux (WAF) method as applied to a model hyperbolic conservation law(Cranfield Institute of Technology; College of Aeronautics, 1989) Toro, E. F.Item Open Access Viscous flux limiters(1991) Toro, E. F.We present Numerical Viscosity Functions, or NVFs, for use with Riemann-problem based shock-capturing methods as applied to viscous flows. In particular, viscous flux limiters are derived. The analysis pertain to a linear convection-diffusion model equation. Our NVFs combine the physical viscosity, the role of which is maximised, with numerical viscosity, whose role is minimised, to capture TVD solutions to viscous flows