Discrete updates of numerical partial differential equations (PDEs) rely on two branches temporal integration. The first branch is the widely-adopted, traditionally popular approach method-of-lines (MOL) formulation, in which multi-stage Runge-Kutta (RK) methods have shown great success solving ordinary (ODEs) at high-order accuracy. clear separation between and spatial discretizations governin...

Compact Approximate Taylor (CAT) methods for systems of conservation laws were introduced by Carrillo and Pares in 2019. These methods, based on a strategy that allows one to extend high-order Lax-Wendroff nonlinear without using the Cauchy-Kovalevskaya procedure, have arbitrary even order accuracy 2p use (2p + 1)-point stencils, where p is an positive integer. More recently 2021 Carrillo, Macc...

A mathematical model is developed of a spring-loaded pressure relief valve connected to a reservoir of compressible fluid via a single, straight pipe. The valve is modelled using Newtonian mechanics, under assumptions that the reservoir pressure is sufficient to ensure choked flow conditions. The usual assumptions of ideal gas theory lead to a system of first-order partial differential equation...

Many of the recently developed high-resolution schemes for hyperbolic conservation laws are based on upwind di erencing. The building block of these schemes is the averaging of an approximate Godunov solver; its time consuming part involves the eld-byeld decomposition which is required in order to identify the \direction of the wind." Instead, we propose to use as a building block the more robu...

We develop a high order numerical boundary condition for compressible inviscid flows involving complex moving geometries. It is based on finite difference methods on fixed Cartesian meshes which pose a challenge that the moving boundaries intersect the grid lines in an arbitrary fashion. Our method is an extension of the so-called inverse Lax-Wendroff procedure proposed in [16] for conservation...

This paper is devoted to the construction of numerical fluxes for hyperbolic systems. We first present a GFORCE numerical flux, which is a weighted average of the Lax-Friedrichs and Lax-Wendroff fluxes. For the linear advection equation with constant coefficient, the new flux reduces identically to that of the Godunov first order upwind method. Then we incorporate GFORCE in the framework of the...

We have developed a new computer program in Fortran 90, in order to obtain numerical solutions of a system of Relativistic Magnetohydrodynamics partial differential equations with predetermined gravitation (GRMHD), capable of simulating the formation of relativistic jets from the accretion disk of matter up to his ejection. Initially we carried out a study on numerical methods of unidimensional...

Classical and recent numerical schemes for solving hyperbolic conservation laws were analyzed for computational efficiency and application to nonideal gas flows. The Roe-Pike approximate Riemann solver with entropy correction, the Harten second-order scheme and the extension of the Roe-Pike method to secondorder by the MUSCL strategy were compared for one-dimensional flows of an ideal gas. Thes...

In this paper, we propose a high order residual distribution conservative finite difference scheme for solving convection– diffusion equations on non-smooth Cartesian meshes. WENO (weighted essentially non-oscillatory) integration and linear interpolation for the derivatives are used to compute the numerical fluxes based on the point values of the solution. The objective is to obtain a high ord...

Gas-kinetic schemes based on the BGK model are proposed as an alternative evolution model which can cure some of the limitations of current Riemann solvers. To analyse the schemes, simple advection equations are reconstructed and solved using the gas-kinetic BGK model. Results for gas-dynamic application are also presented. The ®nal ̄ux function derived in this model is a combination of a gas-ki...