Search

Lax-Wendroff flux reconstruction (LWFR) schemes have high order of accuracy in both space and time despite having a single internal time step. Here, we design a Jacobian-free LWFR type scheme to solve the special relativistic hydrodynamics equations. We then blend the scheme with a first-order finite volume scheme to control the oscillations near discontinuities. We also use a scaling limiter to preserve the physical admissibility of the solution after ensuring the scheme is admissible in means. A particular focus is given to designing a discontinuity indicator model to detect the local non-smoothness in the solution of the highly non-linear relativistic hydrodynamics equations. Finally, we present the numerical results of a wide range of test cases with fourth and fifth-order schemes to show their robustness and efficiency.

We extend the fourth order, two stage Multi-Derivative Runge Kutta (MDRK) scheme to the Flux Reconstruction (FR) framework by writing both stages in terms of a time averaged flux and then using the approximate Lax-Wendroff procedure to compute the time averaged flux. Numerical flux is carefully constructed to enhance Fourier CFL stability and accuracy. A subcell based blending limiter is developed for the MDRK scheme which ensures that the limited scheme is provably admissibility preserving. Along with being admissibility preserving, the blending scheme is constructed to minimize dissipation errors by using Gauss-Legendre solution points and performing MUSCL-Hancock reconstruction on subcells. The accuracy enhancement of the blending scheme is numerically verified on compressible Euler's equations, with test cases involving shocks and small-scale structures., Comment: arXiv admin note: text overlap with arXiv:2305.10781, arXiv:2207.02954; text overlap with arXiv:1609.04491 by other authors

ADER (Arbitrary high order by DERivatives) and Lax-Wendroff (LW) schemes are two high order single stage methods for solving time dependent partial differential equations. ADER is based on solving a locally implicit equation to obtain a space-time predictor solution while LW is based on an explicit Taylor's expansion in time. We cast the corrector step of ADER Discontinuous Galerkin (DG) scheme into an equivalent quadrature free Flux Reconstruction (FR) framework and then show that the obtained ADER-FR scheme is equivalent to the LWFR scheme with D2 dissipation numerical flux for linear problems. This also implies that the two schemes have the same Fourier stability limit for time step size. The equivalence is verified by numerical experiments.

Lax-Wendroff Flux Reconstruction (LWFR) is a single-stage, high order, quadrature free method for solving hyperbolic conservation laws. This work extends the LWFR scheme to solve conservation laws on curvilinear meshes with adaptive mesh refinement (AMR). The scheme uses a subcell based blending limiter to perform shock capturing and exploits the same subcell structure to obtain admissibility preservation on curvilinear meshes. It is proven that the proposed extension of LWFR scheme to curvilinear grids preserves constant solution (free stream preservation) under the standard metric identities. For curvilinear meshes, linear Fourier stability analysis cannot be used to obtain an optimal CFL number. Thus, an embedded-error based time step computation method is proposed for LWFR method which reduces fine-tuning process required to select a stable CFL number using the wave speed based time step computation. The developments are tested on compressible Euler's equations, validating the blending limiter, admissibility preservation, AMR algorithm, curvilinear meshes and error based time stepping.

This work introduces an extension of the high order, single stage Lax-Wendroff Flux Reconstruction (LWFR) of Babbar et al., JCP (2022) to solve second order time-dependent partial differential equations in conservative form on curvilinear meshes. The method uses BR1 scheme to reduce the system to first order so that the earlier LWFR scheme can be applied. The work makes use of the embedded error-based time stepping introduced in Babbar, Chandrashekar (2024) which becomes particularly relevant in the absence of CFL stability limit for parabolic equations. The scheme is verified to show optimal order convergence and validated with transonic flow over airfoil and unsteady flow over cylinder., Comment: 16 Pages, 8 figures

Lax-Wendroff Flux Reconstruction (LWFR) is a single-stage, high order, quadrature free method for solving hyperbolic conservation laws. We perform a cell average decomposition of the LWFR scheme that is similar to the one used in the admissibility preserving framework of Zhang and Shu (2010). By performing a flux limiting of the time averaged numerical flux, the decomposition is used to obtain an admissibility preserving LWFR scheme. The admissibility preservation framework is further extended to a newly proposed extension of LWFR scheme for conservation laws with source terms. This is the first extension of the high order LW scheme that can handle source terms. The admissibility and accuracy are verified by numerical experiments on the Ten Moment equations of Livermore et al., Comment: 16 pages, 11 figures, accepted for publication in proceedings of ICOSAHOM2023

Lax-Wendroff Flux Reconstruction (LWFR) is a single-stage, high order, quadrature free method for solving hyperbolic conservation laws. We develop a subcell based limiter by blending LWFR with a lower order scheme, either first order finite volume or MUSCL-Hancock scheme. While the blending with a lower order scheme helps to control oscillations, it may not guarantee admissibility of discrete solution, e.g., positivity property of quantities like density and pressure. By exploiting the subcell structure and admissibility of lower order schemes, we devise a strategy to ensure that the blended scheme is admissibility preserving for the mean values and then use a scaling limiter to obtain admissibility of the polynomial solution. For MUSCL-Hancock scheme on non-cell-centered subcells, we develop a slope limiter, time step restrictions and suitable blending of higher order fluxes, that ensures admissibility of lower order updates and hence that of the cell averages. By using the MUSCL-Hancock scheme on subcells and Gauss-Legendre points in flux reconstruction, we improve small-scale resolution compared to the subcell-based RKDG blending scheme with first order finite volume method and Gauss-Legendre-Lobatto points. We demonstrate the performance of our scheme on compressible Euler's equations, showcasing its ability to handle shocks and preserve small-scale structures.

The Lax-Wendroff method is a single step method for evolving time dependent solutions governed by partial differential equations, in contrast to Runge- Kutta methods that need multiple stages per time step. We develop a flux reconstruction version of the method in combination with a Jacobian-free Lax- Wendroff procedure that is applicable to general hyperbolic conservation laws. The method is of collocation type, is quadrature free and can be cast in terms of matrix and vector operations. Special attention is paid to the construction of numerical flux, including for non-linear problems, resulting in higher CFL numbers than existing methods, which is shown through Fourier analysis and yielding uniform performance at all orders. Numerical results up to fifth order of accuracy for linear and non-linear problems are given to demonstrate the performance and accuracy of the method.