High-order residual distribution scheme for the time-dependent Euler equations of fluid dynamics

In the present work, a high order finite element type residual distribution scheme is designed in the framework of multidimensional compressible Euler equations of gas dynamics. The strengths of the proposed approximation rely on the generic spatial discretization of the model equations using a continuous finite element type approximation technique, while avoiding the solution of a large linear system with a sparse mass matrix which would come along with any standard ODE solver in a classical finite element approach to advance the solution in time. In this work, we propose a new Residual Distribution (RD) scheme, which provides an arbitrary explicit high order approximation of the smooth solutions of the Euler equations both in space and time. The design of the scheme via the coupling of the RD formulation \cite{mario,abg} with a Deferred Correction (DeC) type method \cite{shu-dec,Minion2}, allows to have the matrix associated to the update in time, which needs to be inverted, to be diagonal. The use of Bernstein polynomials as shape functions, guarantees that this diagonal matrix is invertible and ensures strict positivity of the resulting diagonal matrix coefficients. This work is the extension of \cite{enumath,Abgrall2017} to multidimensional systems. We have assessed our method on several challenging benchmark problems for one- and two-dimensional Euler equations and the scheme has proven to be robust and to achieve the theoretically predicted high order of accuracy on smooth solutions

Space-time POD-Galerkin approach for parametric flow control

In this contribution we propose reduced order methods to fast and reliably solve parametrized optimal control problems governed by time dependent nonlinear partial differential equations. Our goal is to provide a tool to deal with the time evolution of several nonlinear optimality systems in many-query context, where a system must be analysed for various physical and geometrical features. Optimal control can be used in order to fill the gap between collected data and mathematical model and it is usually related to very time consuming activities: inverse problems, statistics, etc. Standard discretization techniques may lead to unbearable simulations for real applications. We aim at showing how reduced order modelling can solve this issue. We rely on a space-time POD-Galerkin reduction in order to solve the optimal control problem in a low dimensional reduced space in a fast way for several parametric instances. The proposed algorithm is validated with a numerical test based on environmental sciences: a reduced optimal control problem governed by viscous Shallow Waters Equations parametrized not only in the physics features, but also in the geometrical ones. We will show how the reduced model can be useful in order to recover desired velocity and height profiles more rapidly with respect to the standard simulation, not losing accuracy

High-order well-balanced schemes and applications to non-equilibrium flow

The appearance of the source terms in modeling non-equilibrium flow problems containing finite-rate chemistry or combustion poses additional numerical difficulties beyond that for solving non-reacting flows. A well-balanced scheme, which can preserve certain non-trivial steady state solutions exactly, may help minimize some of these difficulties. In this paper, a simple one-dimensional non-equilibrium model with one temperature is considered. We first describe a general strategy to design high-order well-balanced finite-difference schemes and then study the well-balanced properties of the high-order finite-difference weighted essentially non-oscillatory (WENO) scheme, modified balanced WENO schemes and various total variation diminishing (TVD) schemes. The advantages of using a well-balanced scheme in preserving steady states and in resolving small perturbations of such states will be shown. Numerical examples containing both smooth and discontinuous solutions are included to verify the improved accuracy, in addition to the well-balanced behavior

The Simple Finite Volume Lax-Wendroff Weighted Essentially Nonoscillatory Schemes for Shallow Water Equations with Bottom Topography

A Lax-Wendroff-type procedure with the high order finite volume simple weighted essentially nonoscillatory (SWENO) scheme is proposed to simulate the one-dimensional (1D) and two-dimensional (2D) shallow water equations with topography influence in source terms. The system of shallow water equations is discretized using the simple WENO scheme in space and Lax-Wendroff scheme in time. The idea of Lax-Wendroff time discretization can avoid part of characteristic decomposition and calculation of nonlinear weights. The type of simple WENO was first developed by Zhu and Qiu in 2016, which is more simple than classical WENO fashion. In order to maintain good, high resolution and nonoscillation for both continuous and discontinuous flow and suit problems with discontinuous bottom topography, we use the same idea of SWENO reconstruction for flux to treat the source term in prebalanced shallow water equations. A range of numerical examples are performed; as a result, comparing with classical WENO reconstruction and Runge-Kutta time discretization, the simple Lax-Wendroff WENO schemes can obtain the same accuracy order and escape nonphysical oscillation adjacent strong shock, while bringing less absolute truncation error and costing less CPU time for most problems. These conclusions agree with that of finite difference Lax-Wendroff WENO scheme for shallow water equations, while finite volume method has more flexible mesh structure compared to finite difference method

High-Order Residual-Distribution Hyperbolic Advection-Diffusion Schemes: 3rd-, 4th-, and 6th-Order

In this paper, spatially high-order Residual-Distribution (RD) schemes using the first-order hyperbolic system method are proposed for general time-dependent advection-diffusion problems. The corresponding second-order time-dependent hyperbolic advection- diffusion scheme was first introduced in [NASA/TM-2014-218175, 2014], where rapid convergences over each physical time step, with typically less than five Newton iterations, were shown. In that method, the time-dependent hyperbolic advection-diffusion system (linear and nonlinear) was discretized by the second-order upwind RD scheme in a unified manner, and the system of implicit-residual-equations was solved efficiently by Newton's method over every physical time step. In this paper, two techniques for the source term discretization are proposed; 1) reformulation of the source terms with their divergence forms, and 2) correction to the trapezoidal rule for the source term discretization. Third-, fourth, and sixth-order RD schemes are then proposed with the above techniques that, relative to the second-order RD scheme, only cost the evaluation of either the first derivative or both the first and the second derivatives of the source terms. A special fourth-order RD scheme is also proposed that is even less computationally expensive than the third-order RD schemes. The second-order Jacobian formulation was used for all the proposed high-order schemes. The numerical results are then presented for both steady and time-dependent linear and nonlinear advection-diffusion problems. It is shown that these newly developed high-order RD schemes are remarkably efficient and capable of producing the solutions and the gradients to the same order of accuracy of the proposed RD schemes with rapid convergence over each physical time step, typically less than ten Newton iterations

A simple and efficient unstructured finite volume scheme for solving the shallow water equations in overland flow applications

An edited version of this paper was published by AGU. Copyright (2015) American Geophysical Union.This paper presents the Decoupled Hydrological Discretisation (DHD) scheme for solving the shallow water equations in hydrological applications involving surface runoff in rural and urban basins. The name of the scheme is motivated by the fact that the three equations which form the two-dimensional shallow water system are discretised independently from each other and thus, the numerical scheme is decoupled in a mathematical sense. Its main advantages compared to other classic finite volume schemes for the shallow water equations are its simplicity to code and the lower computational cost per time step. The validation of the scheme is presented in five test cases involving overland flow and rainfall-runoff transformation over topographies of different complexity. The scheme is compared to the finite volume scheme ofRoe [1986], to the simple inertia formulation [Bates et al., 2010], and to the diffusive wave model. The test cases show that the DHD scheme is able to compute subcritical and supercritical flows in rural and urban environments, and that in overland flow applications it gives similar results to the second order scheme of Roe with a lower computational cost. The results obtained with the simple inertia and diffusive wave models are very similar to those obtained with the DHD scheme in rural basins in which the bed friction and topography dominate the flow hydrodynamics but they deteriorate in typical urban configurations in which the presence of supercritical flow conditions and small scale patterns boost the relevance of the inertial terms in the momentum equations.Peer ReviewedPostprint (published version

Application of conservative residual distribution schemes to the solution of the shallow water equations on unstructured meshes

International audienc

Multi-dimensional higher resolution methods for flow in porous media.

Currently standard first order single-point upstream weighting methods are employed in reservoir simulation for integrating the essentially hyperbolic system components. These methods introduce both coordinate-line numerical diffusion (even in 1-D) and cross-wind diffusion into the solution that is grid and geometry dependent. These effects are particularly important when steep fronts and shocks are present and for cases where flow is across grid coordinate lines. In this thesis, families of novel edge-based and cell-based truly multidimensional upwind formulations that upwind in the direction of the wave paths in order to minimise crosswind diffusion are presented for hyperbolic conservation laws on structured and unstructured triangular and quadrilateral grids in two dimensions. Higher resolution as well as higher order multidimensional formulations are also developed for general structured and unstructured grids. The schemes are coupled with existing consistent and efficient continuous CVD (MPFA) Darcy flux approximations. They are formulated using an IMPES (Implicit in Pressure Explicit in Saturation) strategy for solving the coupled elliptic (pressure) and hyperbolic (saturation) system of equations governing the multi-phase multi-component flow in porous media. The new methods are compared with single point upstream weighting for two-phase and three-component two-phase flow problems. The tests arc conducted on both structured and unstructured grids and involve full-tensor coefficient velocity fields in homogeneous and heterogeneous domains. The comparisons demonstrate the benefits of multidimensional and higher order multidimensional schemes in terms of improved front resolution together with significant reduction in cross-wind diffusion