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.

A finite element method for computing viscous incompressible flows based on the gauge formulation introduced in [Weinan E, Liu J-G. Gauge method for viscous incompressible flows. Journal of Computational Physics (submitted)] is presented. This formulation replaces the pressure by a gauge variable. This new gauge variable is a numerical tool and differs from the standard gauge variable that arises from decomposing a compressible velocity field. It has the advantage that an additional boundary condition can be assigned to the gauge variable, thus eliminating the issue of a pressure boundary condition associated with the original primitive variable formulation. The computational task is then reduced to solving standard heat and Poisson equations, which are approximated by straightforward, piecewise linear (or higher-order) finite elements. This method can achieve high-order accuracy at a cost comparable with that of solving standard heat and Poisson equations. It is naturally adapted to complex geometry and it is much simpler than traditional finite element methods for incompressible flows. Several numerical examples on both structured and unstructured grids are presented. Copyright

We propose a generalized multiscale finite element method (GMsFEM) based on clustering algorithm to study the elliptic PDEs with random coefficients in the multi-query setting. Our method consists of offline and online stages. In the offline stage, we construct a small number of reduced basis functions within each coarse grid block, which can then be used to approximate the multiscale finite element basis functions. In addition, we coarsen the corresponding random space through a clustering algorithm. In the online stage, we can obtain the multiscale finite element basis very efficiently on a coarse grid by using the pre-computed multiscale basis. The new GMsFEM can be applied to multiscale SPDE starting with a relatively coarse grid, without requiring the coarsest grid to resolve the smallest-scale of the solution. The new method offers considerable savings in solving multiscale SPDEs. Numerical results are

No abstract uploaded!

In this paper, we apply a simple finite element numerical scheme, proposed in an earlier work (Liu in Math Comput 70(234):579-593, 2000), to perform a high resolution numerical simulation of incompressible flow over an irregular domain and analyze its boundary layer separation. Compared with many classical finite element fluid solvers, this numerical method avoids a Stokes solver, and only two Poisson-like equations need to be solved at each time step/stage. In addition, its combination with the fully explicit fourth order Runge-Kutta (RK4) time discretization enables us to compute high Reynolds number flow in a very efficient way. As an application of this robust numerical solver, the dynamical mechanism of the boundary layer separation for a triangular cavity flow with Reynolds numbers \(Re=10^4\) and \(Re=10^5\) , including the precise values of bifurcation location and critical time, are reported in this paper. In addition, we provide a super-convergence analysis for the simple finite element numerical scheme, using linear elements over a uniform triangulation with right triangles.