Siddhartha Mishra

Optimal entropy solutions for conservation laws with discontinuous flux-functions

We deal with a single conservation law in one space dimension whose flux function is discontinuous in the space variable and we introduce a proper framework of entropy solutions. We consider a large class of fluxes, namely, fluxes of the convex-convex type and of the concave-convex (mixed) type. The alternative entropy framework that is proposed here is based on a two step approach. In the first step, infinitely many classes of entropy solutions are defined, each associated with an interface connection. We show that each of these class of entropy solutions form a contractive semigroup in

Arbitrarily High-order Accurate Entropy Stable Essentially Nonoscillatory Schemes for Systems of Conservation Laws

L^1

Well-balanced and energy stable schemes for the shallow water equations with discontinuous topography

and is hence unique. Godunov type schemes based on solutions of the Riemann problem are designed and shown to converge to each class of these entropy solutions. The second step is to choose one of these classes of solutions. This choice depends on the Physics of the problem being considered and we concentrate on the model of two-phase flows in a heterogeneous porous medium. We define an optimizationproblem on the set of admissible interface connections and show the existence of an unique optimal connection and its corresponding optimal entropy solution. The optimal entropy solution is consistent with the expected solutions for two-phase flows in heterogeneous porous media.

Well-balanced schemes for the Euler equations with gravitation

We design arbitrarily high-order accurate entropy stable schemes for systems of conservation laws. The schemes, termed TeCNO schemes, are based on two main ingredients: (i) high-order accurate entropy conservative fluxes and (ii) suitable numerical diffusion operators involving ENO reconstructed cell-interface values of scaled entropy variables. Numerical experiments in one and two space dimensions are presented to illustrate the robust numerical performance of the TeCNO schemes.

Entropy stable shock capturing space---time discontinuous Galerkin schemes for systems of conservation laws

We consider the shallow water equations with non-flat bottom topography. The smooth solutions of these equations are energy conservative, whereas weak solutions are energy stable. The equations possess interesting steady states of lake at rest as well as moving equilibrium states. We design energy conservative finite volume schemes which preserve (i) the lake at rest steady state in both one and two space dimensions, and (ii) one-dimensional moving equilibrium states. Suitable energy stable numerical diffusion operators, based on energy and equilibrium variables, are designed to preserve these two types of steady states. Several numerical experiments illustrating the robustness of the energy preserving and energy stable well-balanced schemes are presented.

Scaling morphogen gradients during tissue growth by a cell division rule

Well-balanced high-order finite volume schemes are designed to approximate the Euler equations with gravitation. The schemes preserve discrete equilibria, corresponding to a large class of physically stable hydrostatic steady states. Based on a novel local hydrostatic reconstruction, the derived schemes are well-balanced for any consistent numerical flux for the Euler equations. The form of the hydrostatic reconstruction is both very simple and computationally efficient as it requires no analytical or numerical integration. Moreover, as required by many interesting astrophysical scenarios, the schemes are applicable beyond the ideal gas law. Both first- and second-order accurate versions of the schemes and their extension to multi-dimensional equilibria are presented. Several numerical experiments demonstrating the superior performance of the well-balanced schemes, as compared to standard finite volume schemes, are also presented.

Construction of Approximate Entropy Measure-Valued Solutions for Hyperbolic Systems of Conservation Laws

We present a streamline diffusion shock capturing spacetime discontinuous Galerkin (DG) method to approximate nonlinear systems of conservation laws in several space dimensions. The degrees of freedom are in terms of the entropy variables and the numerical flux functions are the entropy stable finite volume fluxes. We show entropy stability of the (formally) arbitrarily high order accurate method for a general system of conservation laws. Furthermore, we prove that the approximate solutions converge to the entropy measure valued solutions for nonlinear systems of conservation laws. Convergence to entropy solutions for scalar conservation laws and for linear symmetrizable systems is also shown. Numerical experiments are presented to illustrate the robustness of the proposed schemes.

ENO Reconstruction and ENO Interpolation Are Stable

Morphogen gradients guide the patterning of tissues and organs during the development of multicellular organisms. In many cases, morphogen signaling is also required for tissue growth. The consequences of this interplay between growth and patterning are not well understood. In the Drosophila wing imaginal disc, the morphogen Dpp guides patterning and is also required for tissue growth. In particular, it was recently reported that cell division in the disc correlates with the temporal increase in Dpp signaling. Here we mathematically model morphogen gradient formation in a growing tissue, accounting also for morphogen advection and dilution. Our analysis defines a new scaling mechanism, which we term the morphogen-dependent division rule (MDDR): when cell division depends on the temporal increase in morphogen signaling, the morphogen gradient scales with the growing tissue size, tissue growth becomes spatially uniform and the tissue naturally attains a finite size. This model is consistent with many properties of the wing disc. However, we find that the MDDR is not consistent with the phenotype of scaling-defective mutants, supporting the view that temporal increase in Dpp signaling is not the driver of cell division during late phases of disc development. More generally, our results show that local coupling of cell division with morphogen signaling can lead to gradient scaling and uniform growth even in the absence of global feedbacks. The MDDR scaling mechanism might be particularly beneficial during rapid proliferation, when global feedbacks are hard to implement.

Multilevel Monte Carlo Finite Volume Methods for Shallow Water Equations with Uncertain Topography in Multi-dimensions

Entropy solutions have been widely accepted as the suitable solution framework for systems of conservation laws in several space dimensions. However, recent results in De Lellis and Székelyhidi Jr (Ann Math 170(3):1417–1436, 2009) and Chiodaroli et al. (2013) have demonstrated that entropy solutions may not be unique. In this paper, we present numerical evidence that state-of-the-art numerical schemes need not converge to an entropy solution of systems of conservation laws as the mesh is refined. Combining these two facts, we argue that entropy solutions may not be suitable as a solution framework for systems of conservation laws, particularly in several space dimensions. We advocate entropy measure-valued solutions, first proposed by DiPerna, as the appropriate solution paradigm for systems of conservation laws. To this end, we present a detailed numerical procedure which constructs stable approximations to entropy measure-valued solutions, and provide sufficient conditions that guarantee that these approximations converge to an entropy measure-valued solution as the mesh is refined, thus providing a viable numerical framework for systems of conservation laws in several space dimensions. A large number of numerical experiments that illustrate the proposed paradigm are presented and are utilized to examine several interesting properties of the computed entropy measure-valued solutions.

On the upstream mobility scheme for two-phase flow in porous media

We prove that the ENO reconstruction and ENO interpolation procedures are stable in the sense that the jump of the reconstructed ENO point values at each cell interface has the same sign as the jump of the underlying cell averages across that interface. Moreover, we prove that the size of these jumps after reconstruction relative to the jump of the underlying cell averages is bounded. Similar sign properties and the boundedness of the jumps hold for the ENO interpolation procedure. These estimates, which are shown to hold for ENO reconstruction and interpolation of arbitrary order of accuracy and on nonuniform meshes, indicate a remarkable rigidity of the piecewise polynomial ENO procedure.