A high-order nonconservative approach for hyperbolic equations in fluid dynamics

It is well known, thanks to Lax-Wendroff theorem, that the local conservation of a numerical scheme for a conservative hyperbolic system is a simple and systematic way to guarantee that, if stable, a scheme will provide a sequence of solutions that will converge to a weak solution of the continuous problem. In [1], it is shown that a nonconservative scheme will not provide a good solution. The question of using, nevertheless, a nonconservative formulation of the system and getting the correct solution has been a long-standing debate. In this paper, we show how get a relevant weak solution from a pressure-based formulation of the Euler equations of fluid mechanics. This is useful when dealing with nonlinear equations of state because it is easier to compute the internal energy from the pressure than the opposite. This makes it possible to get oscillation free solutions, contrarily to classical conservative methods. An extension to multiphase flows is also discussed, as well as a multidimensional extension

Arbitrary-Lagrangian-Eulerian One-Step WENO Finite Volume Schemes on Unstructured Triangular Meshes

In this article we present a new class of high order accurate Arbitrary-Eulerian-Lagrangian (ALE) one-step WENO finite volume schemes for solving nonlinear hyperbolic systems of conservation laws on moving two dimensional unstructured triangular meshes. A WENO reconstruction algorithm is used to achieve high order accuracy in space and a high order one-step time discretization is achieved by using the local space-time Galerkin predictor. For that purpose, a new element--local weak formulation of the governing PDE is adopted on moving space--time elements. The space-time basis and test functions are obtained considering Lagrange interpolation polynomials passing through a predefined set of nodes. Moreover, a polynomial mapping defined by the same local space-time basis functions as the weak solution of the PDE is used to map the moving physical space-time element onto a space-time reference element. To maintain algorithmic simplicity, the final ALE one-step finite volume scheme uses moving triangular meshes with straight edges. This is possible in the ALE framework, which allows a local mesh velocity that is different from the local fluid velocity. We present numerical convergence rates for the schemes presented in this paper up to sixth order of accuracy in space and time and show some classical numerical test problems for the two-dimensional Euler equations of compressible gas dynamics.Comment: Accepted by "Communications in Computational Physics

Observations on the radiative corrections to pion beta-decay

We find that, in the local V-A theory, the radiative corrections to pion beta-decay involving the weak vector current, when evaluated in the current algebra (CA) formulation in which quarks are the fundamental underlying fields, show a small difference with the more elementary calculations based directly on the pion fields. We show that this difference arises from a specific short-distance effect that depends on the algebra satisfied by the weak and electromagnetic currents. On the other hand, we present a simple theoretical argument that concludes that this difference does not occur when the CA formulation is compared with the chiral perturbation theory (chiPT) approach. Comparisons with previous studies, and with a more recent calculation based on chiPT, are included. We also briefly review the important differences between the results in the local V-A theory and the Standard Model.Comment: 5 pages, 1 figure. V2: two paragraphs have been added in Section III. Final version on PR

Hamiltonian structure of peakons as weak solutions for the modified Camassa-Holm equation

The modified Camassa-Holm (mCH) equation is a bi-Hamiltonian system possessing $N$-peakon weak solutions, for all $N\geq 1$, in the setting of an integral formulation which is used in analysis for studying local well-posedness, global existence, and wave breaking for non-peakon solutions. Unlike the original Camassa-Holm equation, the two Hamiltonians of the mCH equation do not reduce to conserved integrals (constants of motion) for $2$-peakon weak solutions. This perplexing situation is addressed here by finding an explicit conserved integral for $N$-peakon weak solutions for all $N\geq 2$. When $N$ is even, the conserved integral is shown to provide a Hamiltonian structure with the use of a natural Poisson bracket that arises from reduction of one of the Hamiltonian structures of the mCH equation. But when $N$ is odd, the Hamiltonian equations of motion arising from the conserved integral using this Poisson bracket are found to differ from the dynamical equations for the mCH $N$-peakon weak solutions. Moreover, the lack of conservation of the two Hamiltonians of the mCH equation when they are reduced to $2$-peakon weak solutions is shown to extend to $N$-peakon weak solutions for all $N\geq 2$. The connection between this loss of integrability structure and related work by Chang and Szmigielski on the Lax pair for the mCH equation is discussed.Comment: Minor errata in Eqns. (32) to (34) and Lemma 1 have been fixe

Beyond local effective material properties for metamaterials

To discuss the properties of metamaterials on physical grounds and to consider them in applications, effective material parameters are usually introduced and assigned to a given metamaterial. In most cases, only weak spatial dispersion is considered. It allows to assign local material properties, i.e. a permittivity and a permeability. However, this turned out to be insufficient. To solve this problem, we study here the effective properties of metamaterials with constitutive relations beyond a local response and take strong spatial dispersion into account. The isofrequency surfaces of the dispersion relation are investigated and compared to those of an actual metamaterial. The significant improvement provides evidence for the necessity to use nonlocal material laws in the effective description of metamaterials. The general formulation we choose here renders our approach applicable to a wide class of metamaterials