Numerical schemes provably preserving the positivity of density and pressure are highly desirable for ideal magnetohydrodynamics (MHD), but the rigorous positivity-preserving (PP) analysis remains challenging. The difficulties mainly arise from the intrinsic complexity of the MHD equations as well as the indeterminate relation between the PP property and the divergence-free condition on the magnetic field. This paper presents the first rigorous PP analysis of conservative schemes with the Lax--Friedrichs (LF) flux for 1D and multidimensional ideal MHD. The significant innovation is the discovery of the theoretical connection between the PP property and a discrete divergence-free (DDF) condition. This connection is established through the generalized LF splitting properties, which are alternatives to the usually expected LF splitting property that does not hold for ideal MHD. The generalized LF splitting properties involve a number of admissible states strongly coupled by the DDF condition, making their derivation very difficult. We derive these properties via a novel equivalent form of the admissible state set and an important inequality, which is skillfully constructed by technical estimates. Rigorous PP analysis is then presented for finite volume and discontinuous Galerkin schemes with the LF flux on uniform Cartesian meshes. In the 1D case, the PP property is proved for the first-order scheme with proper numerical viscosity, and also for arbitrarily high-order schemes under conditions accessible by a PP limiter. In the 2D case, we show that the DDF condition is necessary and crucial for achieving the PP property. It is observed that even slightly violating the proposed DDF condition may cause failure to preserve the positivity of pressure. We prove that the 2D LF type scheme with proper numerical viscosity preserves both the positivity and the DDF condition. Sufficient conditions are derived for 2D PP high-order schemes, and extension to 3D is discussed. Numerical examples provided in the supplementary material further confirm the theoretical findings.

We establish theoretical recovery guarantees of a family of Riemannian optimization algorithms for low rank matrix recovery, which is about recovering an mimes n rank mimes n matrix from mimes n number of linear measurements. The algorithms are first interpreted as iterative hard thresholding algorithms with subspace projections. Based on this connection, we show that provided the restricted isometry constant mimes n of the sensing operator is less than mimes n, the Riemannian gradient descent algorithm and a restarted variant of the Riemannian conjugate gradient algorithm are guaranteed to converge linearly to the underlying rank mimes n matrix if they are initialized by one step hard thresholding. Empirical evaluation shows that the algorithms are able to recover a low rank matrix from nearly the minimum number of measurements necessary.

We propose a class of simple and efficient numerical scheme for incompressible fluid equations with coordinate symmetry. By introducing a generalized vorticity-stream formulation, the divergence free constraints are automatically satisfied. In addition, with explicit treatment of the nonlinear terms and local vorticity boundary condition, the Navier–Stokes (MHD, respectively) equation essentially decouples into 2 (4, respectively) scalar equation and thus the scheme is very efficient. Moreover, with proper discretization of the nonlinear terms, the scheme preserves both energy and helicity identities numerically. This is achieved by recasting the nonlinear terms (convection, vorticity stretching, geometric source, Lorentz force and electro-motive force) in terms of Jacobians. This conservative property is valid even in the presence of the pole singularity for axisymmetric flows. The exact conservation of energy and helicity has effectively eliminated excessive numerical viscosity. Numerical examples have demonstrated both accuracy and efficiency of the scheme. Finally, local mesh refinement near the boundary can also be easily incorporated into the scheme without extra cost.

No abstract uploaded!

No abstract uploaded!