Search

In this work we develop a new framework to deal numerically with discontinuous solutions in nonconservative hyperbolic systems. First an extension of the MOOD methodology to nonconservative systems based on Taylor expansions is presented. This extension combined with an in-cell discontinuous reconstruction operator are the key points to develop a new family of high-order methods that are able to capture exactly isolated shocks. Several test cases are proposed to validate these methods for the Modified Shallow Water equations and the Two-Layer Shallow Water system., Comment: Preprint. arXiv admin note: text overlap with arXiv:2105.00424

The aim of this work is to design implicit and semi-implicit high-order well-balanced finite-volume numerical methods for 1D systems of balance laws. The strategy introduced by two of the authors in a previous paper for explicit schemes based on the application of a well-balanced reconstruction operator has been applied. The well-balanced property is preserved when quadrature formulas are used to approximate the averages and the integral of the source term in the cells. Concerning the time evolution, this technique is combined with a time discretization method of type RK-IMEX or RK-implicit. The methodology will be applied to several systems of balance laws.

We are interested in the numerical approximation of discontinuous solutions in non conservative hyperbolic systems. An extension to second-order of a new strategy based on in-cell discontinuous reconstructions to deal with this challenging topic is presented. This extension is based on the combination of the first-order in-cell reconstruction and the MUSCL-Hancock reconstruction. The first-order strategy allowed in particular to capture exactly the isolated shocks and this new second-order extension keep this property. We also set the basis of an extension to high-order methods following this in-cell methodology and a way for capturing exactly more than one shock. Several numerical tests are proposed to validate the methods for the Coupled-Burgers system, Gas dynamics equations in Lagrangian coordinates and the modified shallow water system.

For the evolution of a compressible fluid in spherical symmetry on a Schwarzschild curved background, we design a class of well-balanced numerical algorithms with first-order or second-order of accuracy. We treat both the relativistic Burgers-Schwarzschild model and the relativistic Euler-Schwarzschild model and take advantage of the explicit or implicit forms available for the stationary solutions of these models. Our schemes follow the finite volume methodology and preserve the stationary solutions. Importantly, they allow us to investigate the global asymptotic behavior of such flows and determine the asymptotic behavior of the mass density and velocity field of the fluid., Comment: 40 pages

The goal of this work is to introduce new families of shock-capturing high-order numerical methods for systems of conservation laws that combine Fast WENO (FWENO) and Optimal WENO (OWENO) reconstructions with Approximate Taylor methods for the time discretization. FWENO reconstructions are based on smoothness indicators that require a lower number of calculations than the standard ones. OWENO reconstructions are based on a definition of the nonlinear weights that allows one to unconditionally attain the optimal order of accuracy regardless of the order of critical points. Approximate Taylor methods update the numerical solutions by using a Taylor expansion in time in which, instead of using the Cauchy-Kovalevskaya procedure, the time derivatives are computed by combining spatial and temporal numerical differentiation with Taylor expansions in a recursive way. These new methods are compared between them and against methods based on standard WENO implementations and/or SSP-RK time discretization. A number of test cases are considered ranging from scalar linear 1d problems to nonlinear systems of conservation laws in 2d.

In this paper, high order well-balanced finite difference weighted essentially non-oscillatory methods to solve general systems of balance laws are presented. Two different families are introduced: while the methods in the first one preserve every stationary solution, those in the second family only preserve a given set of stationary solutions that depend on some parameters. The accuracy, well-balancedness, and conservation properties of the methods are discussed, as well as their application to systems with singular source terms. The strategy is applied to derive third and fifth order well-balanced methods for a linear scalar balance law, Burgers' equation with a nonlinear source term, and for the shallow water model. In particular, numerical methods that preserve every stationary solution or only water at rest equilibria are derived for the latter.

A new family of high order methods for systems of conservation laws are introduced: the Compact Approximate Taylor (CAT) methods. These methods are based on centered (2p + 1)-point stencils where p is an arbitrary integer. We prove that the order of accuracy is 2p and that CAT methods are an extension of high-order Lax-Wendroff methods for linear problems. Due to this, they are linearly L2-stable under a CFL<1 condition. In order to prevent the spurious oscillations that appear close to discontinuities two shock-capturing techniques have been considered: a fux-limiter technique (FL-CAT methods) and WENO reconstruction for the frst time derivative (WENO-CAT methods). We follow the WENO-Lax Wendroff Approximate Taylor method of Zorio, Baeza and Mullet (2017) in the second approach. A number of test cases are considered to compare these methods with other WENO-based schemes: the linear transport equation, Burgers equation, and the 1D compressible Euler system are considered. Although CAT methods present an extra computational cost due to the local character, this extra cost is compensated by the fact that they still give good solutions with CFL values close to 1., Comment: 10 figures, 33 pages

We are interested in nonlinear hyperbolic systems in nonconservative form arising in fluid dynamics, and, for solutions containing shock waves, we investigate the convergence of finite difference schemes applied to such systems. According to Dal Maso, LeFloch, and Murat's theory, a shock wave theory for a given nonconservative system requires prescribing a priori a family of paths in the phase space. In the present paper, we consider schemes that are formally consistent with a given family of paths, and we investigate their limiting behavior as the mesh is refined. We generalize to systems a property established earlier by Hou and LeFloch for scalar conservation laws, and we prove that nonconservative schemes generate, at the level of the limiting hyperbolic system, a "convergence error" source-term which, provided the total variation of the approximations remains uniformly bounded, is a locally bounded measure. We discuss the role of the equivalent equation associated with a difference scheme; here, the distinction between scalar equations and systems appears most clearly since, for systems, the equivalent equation of a scheme that is formally path-consistent depends upon the prescribed family of paths. The core of this paper is devoted to investigate numerically the approximation of several models arising in fluid dynamics. For systems having nonconservative products associated with linearly degenerate characteristic fields, the convergence error vanishes. For some other models, this measure is evaluated very accurately, especially by plotting the shock curves associated with each scheme under consideration., Comment: 34 pages, 16 figures