The problem of inverting the total divergence operator is central to finding components of a given conservation law. This might not be taxing for a low-order conservation law of a scalar partial differential equation, but integrable systems have conservation laws of arbitrarily high order that must be found with the aid of computer algebra. Even low-order conservation laws of complex systems can be hard to find and invert. This paper describes a new, efficient approach to the inversion problem. Two main tools are developed: partial Euler operators and partial scalings. These lead to a line integral formula for the inversion of a total derivative and a procedure for inverting a given total divergence concisely.

The present paper concerns the system ut + [ϕ(u)]x = 0, vt + [ψ(u)v]x = 0 having distributions as initial conditions. Under certain conditions, and supposing ϕ, ψ: ℝ → ℝ functions, we explicitly solve this Cauchy problem within a convenient space of distributions u,v. For this purpose, a consistent extension of the classical solution concept defined in the setting of a distributional product (not constructed by approximation processes) is used. Shock waves, δ-shock waves, and also waves defined by distributions that are not measures are presented explicitly as examples. This study is carried out without assuming classical results about conservation laws. For reader's convenience, a brief survey of the distributional product is also included.

We prove the stability with respect to the flux of solutions to initial – boundary value problems for scalar non autonomous conservation laws in one space dimension. Key estimates are obtained through a careful construction of the solutions.

Given a vector field on a manifold $M$, we define a globally conserved quantity to be a differential form whose Lie derivative is exact. Integrals of conserved quantities over suitable submanifolds are constant under time evolution, the Kelvin circulation theorem being a well-known special case. More generally, conserved quantities are well behaved under transgression to spaces of maps into $M$. We focus on the case of multisymplectic manifolds and Hamiltonian vector fields. Our main result is that in the presence of a Lie group of symmetries admitting a homotopy co-momentum map, one obtains a whole family of globally conserved quantities. This extends a classical result in symplectic geometry. We carry this out in a general setting, considering several variants of the notion of globally conserved quantity.

We propose a fully conservative and less oscillatory multi-moment scheme for the approximation of hyperbolic conservation laws. The proposed scheme (CIP-CSL3ENO) is based on two CIP-CSL3 schemes and the essentially non-oscillatory (ENO) scheme. In this paper, we also propose an ENO indicator for the multimoment framework, which intentionally selects non-smooth stencil but can efficiently minimize numerical oscillations. The proposed scheme is validated through various benchmark problems and a comparison with an experiment of two droplets collision/separation. The CIP-CSL3ENO scheme shows approximately fourth-order accuracy for smooth solution, and captures discontinuities and smooth solutions simultaneously without numerical oscillations for solutions which include discontinuities. The numerical results of two droplets collision/separation (3D free surface flow simulation) show that the CIP-CSL3ENO scheme can be applied to various types of fluid problems not only compressible flow problems but also incompressible and 3D free surface flow problems.

The nonlinear Dirac equation is an important model in quantum physics with a set of conservation laws and a multi-symplectic formulation. In this paper, we propose energy-preserving and multi-symplectic wavelet algorithms for this model. Meanwhile, we evidently improve the efficiency of these algorithms in computations via splitting technique and explicit strategy. Numerical experiments are conducted during long-term simulations to show the excellent performances of the proposed algorithms and verify our theoretical analysis.

The nonlinear and weakly dispersive Serre equations contain higher-order dispersive terms. These include mixed spatial and temporal derivative flux terms which are difficult to handle numerically. These terms can be replaced by an alternative combination of equivalent temporal and spatial terms, so that the Serre equations can be written in conservation law form. The water depth and new conserved quantities are evolved using a second-order finite-volume scheme. The remaining primitive variable, the depth-averaged horizontal velocity, is obtained by solving a second-order elliptic equation using simple finite differences. Using an analytical solution and simulating the dam-break problem, the proposed scheme is shown to be accurate, simple to implement and stable for a range of problems, including flows with steep gradients. It is only slightly more computationally expensive than solving the shallow water wave equations.

In [SIAM J. Sci. Comput., 35(2)(2013), A1049-A1072], a class of multi-domain hybrid DG and WENO methods for conservation laws was introduced. Recent applications of this method showed that numerical instability may encounter if the DG flux with Lagrangian interpolation is applied as the interface flux during the moment of conservative coupling. In this continuation paper, we present a more robust approach in the construction of DG flux at the coupling interface by using WENO procedures of reconstruction. Based on this approach, such numerical instability is overcome very well. In addition, the procedure of coupling a DG method with a WENO-FD scheme on hybrid meshes is disclosed in detail. Typical testing cases are employed to demonstrate the accuracy of this approach and the stability under the flexibility of using either WENO-FD flux or DG flux at the moment of requiring conservative coupling.

For scalar conservation laws in one space dimension with a flux function discontinuous inspace, there exist infinitely many classes of solutions which are L1 contractive.Each class is characterized by a connection (A,B) which determines the interface entropy. Forsolutions corresponding to a connection (A,B), there exists convergent numerical schemesbased on Godunov or Engquist−Osher schemes. The natural question is how to obtain schemes,corresponding to computationally less expensive monotone schemes like Lax−Friedrichs etc., usedwidely in applications. In this paper we completely answer this question for more general(A,B)stable monotone schemes using a novel construction of interface flux function. Then fromthe singular mapping technique of Temple and chain estimate of Adimurthi and Gowda, weprove the convergence of the schemes.

The local one-dimensional multisymplectic scheme (LOD-MS) is developed for the three-dimensional (3D) Gross-Pitaevskii (GP) equation in Bose-Einstein condensates. The idea is originated from the advantages of multisymplectic integrators and from the cheap computational cost of the local one-dimensional (LOD) method. The 3D GP equation is split into three linear LOD Schrödinger equations and an exactly solvable nonlinear Hamiltonian ODE. The three linear LOD Schrödinger equations are multisymplectic which can be approximated by multisymplectic integrator (MI). The conservative properties of the proposed scheme are investigated. It is mass-preserving. Surprisingly, the scheme preserves the discrete local energy conservation laws and global energy conservation law if the wave function is variable separable. This is impossible for conventional MIs in nonlinear Hamiltonian context. The numerical results show that the LOD-MS can simulate the original problems very well. They are consistent with the numerical analysis.