We present a new formulaiton of the incompressible Navier-Stokes equation in terms of an auxiliary field that differs from the velocity by a gauge transformation. The gauge freedom allows us to assign simple and specific boundary conditions for both the auxiliary field and the gauge field, thus eliminating the issue of pressure boundary condition in the usual primitive variable formulation. The resulting dynamic and kinematic equations can then be solved by standard methods for heat and Poisson equations. A normal mode analysis suggests that in contrast to the classical projection method, the gauge method does not suffer from the problem of numerical boundary layers. Thus the subtleties on the spatial discretization for the projection methods are removed. Consequently, the projection step in the gauge method can be accomplished by standard Poissin solves. We demonstrate the efficiency and accuracy of the gauge method by several numerical examples, including the flow past cylinder.

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.

We use a positive parameter to develop a dierentiable perturbed reconstruction model to solve the L1–L2 magnetic resonance image (MRI) reconstruction problem. We use Bregman iterative formulation to solve the dierentiable perturbed L1–L2 model, and lagged diusivity xed-point iteration to solve the minimization problem in the Bregman iteration. Two Fourier transforms and an inverse Fourier transform are used to accelerate L1–L2 MRI reconstruction. Real MR images are used to test the method in numerical experiments. The results demonstrate that the proposed method is very ecient for L1–L2MRI reconstruction.

In this paper, we propose new spectral viscosity methods based on the generalized Hermite functions for the solution of nonlinear scalar conservation laws in the whole line. It is shown rigorously that these schemes converge to the unique entropy solution by using compensated compactness arguments, under some conditions. The numerical experiments of the inviscid Burger's equation support our result, and it verifies the reasonableness of the conditions.

Eigenvectors and eigenvalues of discrete Laplacians are often used for manifold learning and nonlinear dimensionality reduction. Graph Laplacian is one widely used discrete laplacian on point cloud. It was previously proved by Belkin and Niyogi [4] that the eigenvectors and eigenvalues of the graph Laplacian converge to the eigenfunctions and eigenvalues of the Laplace-Beltrami operator of the manifold in the limit of innitely many data points sampled independently from the uniform distribution over the manifold. Recently, we introduced Point Integral method (PIM) [15] to solve elliptic equations and corresponding eigenvalue problem on point clouds. In this paper, we prove that the eigenvectors and eigenvalues obtained by PIM converge in the limit of innitely many random samples. Moreover, one estimate of the rate of the convergence is given.