Maojun Li

Maojun Li, Jiancheng Wang, Zeyu Xia et al.

In this paper, we derive a new model to simulate the incompressible resistive magnetohydrodynamic (MHD) free surface flow. A thermodynamically consistent diffuse interface method is adopted to characterize the moving interface in the modeling process. The formal convergence of the proposed MHD free surface flow model to the sharp interface model is established via a matched asymptotic argument, and the model can be solved without the need for sophisticated free surface capturing schemes. We design a fully decoupled linear finite element scheme that preserves the divergence-free constraint of the magnetic field at a discrete level. The reliability and robustness of the proposed model and algorithm are validated through numerical investigations of the magnetic damping effect on bubble dynamics. In particular, we provide a quantitative numerical comparison of the present results with those obtained from an inductionless MHD model and a sharp interface arbitrary Lagrangian--Eulerian model.

Jun She, Haiyun Dong, Maojun Li et al.

In this paper, we present a numerical scheme designed for coupled systems of variable-topography shallow water flow and solute transport. By integrating a variable-density system with an expression for relative density of mixtures, a novel formulation of the coupled system is derived. To ensure the well-balanced property, auxiliary variables are introduced to reformulate the variable-density shallow water equations into a new form, which is then discretized using the discontinuous Galerkin (DG) method with the Lax-Friedrichs (LF) flux as the numerical flux. By selecting appropriate values for the auxiliary variables, we demonstrate that the proposed method accurately preserves steady-state solutions under still water conditions, thereby verifying its well-balanced nature. Furthermore, sufficient conditions for preserving the positivity of both water depth and concentration are proposed and rigorously proven. A positivity-preserving limiter is introduced to enforce these conditions. Finally, a series of numerical examples are conducted to validate the computational accuracy and effectiveness of the proposed method.