This paper presents entropy symmetrization and high-order accurate entropy stable schemes for the relativistic magnetohydrodynamic (RMHD) equations. It is shown that the conservative RMHD equations are not symmetrizable and do not admit a thermodynamic entropy pair. To address this issue, a symmetrizable RMHD system, equipped with a convex thermodynamic entropy pair, is first proposed by adding a source term into the equations, providing an analogue to the nonrelativistic Godunov--Powell system. Arbitrarily high-order accurate entropy stable finite difference schemes are developed on Cartesian meshes based on the symmetrizable RMHD system. The crucial ingredients of these schemes include (i) affordable explicit entropy conservative fluxes which are technically derived through carefully selected parameter variables, (ii) a special high-order discretization of the source term in the symmetrizable RMHD system, and (iii) suitable high-order dissipative operators based on essentially nonoscillatory reconstruction to ensure the entropy stability. Several numerical tests demonstrate the accuracy and robustness of the proposed entropy stable schemes.

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.

Suppose that $\eta$ is a Schramm-Loewner evolution ($\SLE_\kappa$) in a smoothly bounded simply connected domain $D \subset \C$ and that $\phi$ is a conformal map from $\D$ to a connected component of $D \setminus \eta([0,t])$ for some $t>0$. The multifractal spectrum of $\eta$ is the function $(-1,1) \to [0,\infty)$ which, for each $s \in (-1,1)$, gives the Hausdorff dimension of the set of points $x \in \partial \D$ such that $|\phi'( (1-\epsilon) x)| = \epsilon^{-s+o(1)}$ as $\epsilon \to 0$. We rigorously compute the a.s.\ multifractal spectrum of $\SLE$, confirming a prediction due to Duplantier. As corollaries, we confirm a conjecture made by Beliaev and Smirnov for the a.s.\ bulk integral means spectrum of $\SLE$, we obtain the optimal H\"older exponent for a conformal map which uniformizes the complement of an $\SLE$ curve, and we obtain a new derivation of the a.s.\ Hausdorff dimension of the $\SLE$ curve for $\kappa \leq 4$. Our results also hold for the $\SLE_\kappa(\ul\rho)$ processes with general vectors of weight $\ul\rho$.

The height fluctuations of the models in the KPZ class are expected to converge to a universal process. The spatial process at equal time is known to converge to the Airy process or its variations. However, the temporal process, or more generally the two-dimensional space-time fluctuation field, is less well understood. We consider this question for the periodic TASEP (totally asymmetric simple exclusion process). For a particular initial condition, we evaluate the multitime and multilocation distribution explicitly in terms of a multiple integral involving a Fredholm determinant. We then evaluate the large-time limit in the so-called relaxation time scale.

No abstract uploaded!