A Well-Balanced HLL Scheme for Hyperbolic Conservation Systems With Source Terms

Mathematical balancing of flux gradient and source terms prior to using Roe’s approximate Riemann solver

The exact Riemann Solver to the Shallow Water equations for natural channels with bottom steps

On many occasions the Navier-Stokes equations are integrated over one or more space directions to remove small-scale features from the solution, such as the size of individual rotor blades in large wind farms, or the micro-pores in flow through porous media. The result of these manipulations is the appearance of additional terms on the right-hand side of the conservation laws. These additional terms may exhibit discontinuous or singular behavior, even if the underlying behavior they represent is regular. Naïve treatment of these additional terms results in non-physical features in the solution or inaccurate predictions of mass, momentum, or energy transfer. This paper proposes a rigorous and straightforward approach for accurately handling source terms with singularities or discontinuities for systems governed by conservation laws, discretized using finite volume schemes with Godunov-type fluxes. We distinguish between volume and surface source terms. Both types can be handled by carefully constructed upwind procedures. This discretization exhibits second-order accuracy, for formally first-order numerical models. The method is validated using the Euler equations with large source terms, with both preconditioned and non-preconditioned flux functions. We also present a novel, elegant method for introducing jump conditions in the solution via singular source terms. This ensures global compatibility relations are met without the need for internal boundaries, as commonly done in actuator disk models. The robustness of the method is demonstrated on a Riemann problem with abrupt area change, matching analytical solutions for 10 different wave patterns. Finally, the new methods are applied to a multi-dimensional, open rotor problem and compared to other results in the literature.

High order well-balanced finite volume WENO schemes and discontinuous Galerkin methods for a class of hyperbolic systems with source terms

Second-order accurate finite volume schemes with the discrete maximum principle for solving Richards’ equation on unstructured meshes

Self-gravitating isothermal supersonic turbulence is analyzed in the asymptotic limit of large Reynolds numbers. Based on the inviscid invariance of total energy, an exact relation is derived for homogeneous (not necessarily isotropic) turbulence. A modified definition for the two-point energy correlation functions is used to comply with the requirement of detailed energy equipartition in the acoustic limit. In contrast to the previous relations (S. Galtier and S. Banerjee, Phys. Rev. Lett. 107, 134501 (2011)PRLTAO0031-900710.1103/PhysRevLett.107.134501; S. Banerjee and S.Galtier, Phys. Rev. E 87, 013019 (2013)PLEEE81539-375510.1103/PhysRevE.87.013019), the current exact relation shows that the pressure dilatation terms play practically no role in the energy cascade. Both the flux and source terms are written in terms of two-point differences. Sources enter the relation in a form of mixed second-order structure functions. Unlike the kinetic and thermodynamic potential energies, the gravitational contribution is absent from the flux term. An estimate shows that, for the isotropic case, the correlation between density and gravitational acceleration may play an important role in modifying the energy transfer in self-gravitating turbulence. The exact relation is also written in an alternative form in terms of two-point correlation functions, which is then used to describe scale-by-scale energy budget in spectral space.

Adaptive quadtree simulation of shallow flows with wet–dry fronts over complex topography

In this work we present some moment preserving finite volume schemes (FVS) for solving population balance equations. We are considering unified numerical methods to simultaneous aggregation, breakage, growth and source terms, e.g. for nucleation. The criteria for the preservation of different moments are given. The property of conservation is a special case of preservation. First we present a FVS which shows the preservation with respect to one-moment depending upon the processes under consideration. In case of the aggregation and breakage it satisfies first-moment preservation whereas for the growth and nucleation we observe zeroth-moment preservation. This is due to the well-known property of conservativity of FVS. However, coupling of all the processes shows no preservation for any moments. To overcome this issue, we reformulate the cell average technique into a conservative formulation which is coupled together with a modified upwind scheme to give moment preservation with respect to the first two-moments for all four processes under consideration. This allows for easy coupling of these processes. The preservation is proven mathematically and verified numerically. The numerical results for the first two-moments are verified for various coupled processes where analytical solutions are available.

2D well-balanced augmented ADER schemes for the Shallow Water Equations with bed elevation and extension to the rotating frame

Scalar hyperbolic balance laws in several space dimensions play a central role in this thesis. First, we deal with a new class of mixed parabolic–hyperbolic systems on all Rn: we obtain the basic well–posedness theorems, devise an ad hoc numerical algorithm, prove its convergence and investigate the qualitative properties of the solutions. The extension of these results to bounded domains requires a deep understanding of the initial boundary value problem (IBVP) for hyperbolic balance laws. The last part of the thesis provides rigorous estimates on the solution to this IBVP, under precise regularity assumptions. In Chapter 1 we introduce a predator–prey model. A non local and non linear balance law is coupled with a parabolic equation: the former describes the evolution of the predator density, the latter that of prey. The two equations are coupled both through the convective part of the balance law and the source terms. The drift term is a non local function of the prey density. This allows the movement of predators to be directed towards the regions where the concentration of prey is higher. We prove the well–posedness of the system, hence the existence and uniqueness of solution, the continuous dependence on the initial data and various stability estimates. In Chapter 2 we devise an algorithm to compute approximate solutions to the mixed system introduced above. The balance law is solved numerically by a Lax–Friedrichs type method via dimensional splitting, while the parabolic equation is approximated through explicit finite–differences. Both source terms are integrated by means of a second order Runge–Kutta scheme. The key result in Chapter 2 is the convergence of this algorithm. The proof relies on a careful tuning between the parabolic and the hyperbolic methods and exploits the non local nature of the convective part in the balance law. This algorithm has been implemented in a series of Python scripts. Using them, we obtain information about the possible order of convergence and we investigate the qualitative properties of the solutions. Moreover, we observe the formation of a striking pattern: while prey diffuse, predators accumulate on the vertices of a regular lattice. The analytic study of the system above is on all Rn. However, both possible biological applications and numerical integrations suggest that the boundary plays a relevant role. With the aim of studying the mixed hyperbolic–parabolic system in a bounded domain, we noticed that for balance laws known results lack some of the estimates necessary to deal with the coupling. In Chapter 3 we then focus on the IBVP for a general balance law in a bounded domain. We prove the well–posedness of this problem, first with homogeneous boundary condition, exploiting the vanishing viscosity technique and the doubling of variables method, then for the non homogeneous case, mainly thanks to elliptic techniques. We pay particular attention to the regularity assumptions and provide rigorous estimates on the solution.

A simple yet precise relation between the flux gradient and the bed slope source term is presented, which produces a net force within the cell with an inclined water surface, but ensures still water condition when there is no flow across the boundaries. The proposed method consists of calculating the pressure term based on the water depths at the cell vertices, which may be computed by a higher order scheme and the bed slope source term by a centered discretization technique. The methodology is demonstrated with a Godunov-type upwind finite volume formulation. The inviscid fluxes are calculated using Roe’s approximate Riemann solver and a second-order spatial accuracy is obtained by implementing multidimensional gradient reconstruction and slope limiting techniques. The accuracy and applicability of the numerical model is verified with a couple of test problems and a real flow example of tidal water movement in a stretch of River Hooghly in India.

An approximate Riemann solver for real gas parabolized Navier–Stokes equations

Treating network junctions in finite volume solution of transient gas flow models

We further investigate the high order positivity-preserving discontinuous Galerkin (DG) methods for linear hyperbolic and radiative transfer equations developed in Yuan et al. (SIAM J Sci Comput 38:A2987–A3019, 2016). The DG methods in Yuan et al. (2016) can maintain positivity and high order accuracy, but they rely both on the scaling limiter in Zhang and Shu (J Comput Phys 229:8918–8934, 2010) and a rotational limiter, the latter may alter cell averages of the unmodulated DG scheme, thereby affecting conservation. Even though a Lax–Wendroff type theorem is proved in Yuan et al. (2016), guaranteeing convergence to weak solutions with correct shock speed when such rotational limiter is applied, it would still be desirable if a conservative DG method without changing the cell averages can be obtained which has both high order accuracy and positivity-preserving capability. In this paper, we develop and analyze such a DG method for both linear hyperbolic equations and radiative transfer equations. In the one-dimensional case, the method uses traditional DG space $$P^k$$ of piecewise polynomials of degree at most k. A key result is proved that the unmodulated DG solver in this case can maintain positivity of the cell average if the inflow boundary value and the source term are both positive, therefore the positivity-preserving framework in Zhang and Shu (2010) can be used to obtain a high order conservative positivity-preserving DG scheme. Unfortunately, in two-dimensions this is no longer the case. We show that the unmodulated DG solver based either on $$P^k$$ or $$Q^k$$ spaces (piecewise kth degree polynomials or piecewise tensor-product kth degree polynomials) could generate negative cell averages. We augment the DG space with additional functions so that the positivity of cell averages from the unmodulated DG solver can be restored, thereby leading to high order conservative positivity-preserving DG scheme based on these augmented DG spaces following the framework in Zhang and Shu (2010). Computational results are provided to demonstrate the good performance of our DG schemes.

In this work, we develop an improved shock-capturing and high-resolution Lagrangian–Eulerian method for hyperbolic systems and balance laws. This is a new method to deal with discontinuous flux and complicated source terms having concentrations for a wide range of applications science and engineering, namely, 1D shallow-water equations, sedimentation processes, Geophysical flows in 2D, N-Wave models, and Riccati-type problems with forcing terms. We also include numerical simulations of a 1D two-phase flow model in porous media, with gravity and a nontrivial singular $$\delta $$ -source term representing an injection point. Moreover, we present approximate solutions for 2D nonlinear systems (Compressible Euler Flows and Shallow-Water Equations) for distinct benchmark configurations available in the literature aiming to present convincing and robust numerical results. In addition, for the linear advection model in 1D and for a smooth solution of the nonlinear Burgers’ problem, second order approximations were obtained. We also present a high-resolution approximation of the nonlinear non-convex Buckley–Leverett problem. Based on the work of A. Harten, we derive a convergent Lagrangian–Eulerian scheme that is total variation diminishing and second-order accurate, away from local extrema and discontinuous data. Additionally, using a suitable Kružkov’s entropy definition, introduced by K. H. Karlsen and J. D. Towers, we can verify that our improved Lagrangian–Eulerian scheme converges to the unique entropy solution for conservation laws with a discontinuous space-time dependent flux. A key hallmark of our method is the dynamic tracking forward of the no-flow curves, which are locally conservative and preserve the natural setting of weak entropic solutions related to hyperbolic problems that are not reversible systems in general. In the end, we have a general procedure to construct a class of Lagrangian–Eulerian schemes to deal with hyperbolic problems with or without forcing terms. The proposed scheme is free of Riemann problem solutions and no adaptive space-time discretizations are needed. The numerical experiments verify the efficiency and accuracy of our new Lagrangian–Eulerian method.