he incompressibility constraint makes Navier-Stokes equations difficult. A reformulation to a better posed problem is needed before solving it numerically. The sequential regularization method (SRM) is a reformulation which combines the penalty method with a stabilization method in the context of constrained dynamical systems and has the benefit of both methods. In the paper, we study the existence and uniqueness for the solution of the SRM and provide a simple proof of the convergence of the solution of the SRM to the solution of the Navier-Stokes equations. We also give error estimates for the time discretized SRM formulation.

Poisson equation on point cloud with Dirichlet boundary condition plays important role in many problems. In this paper, we use the volume constraint proposed by Du et.al to handle the Dirichlet boundary condition in the point integral method for Poisson equation on point cloud. We prove that the solution given by volume constraint converges to the true solution as the point cloud converges to the underlying smooth manifold.

This paper first studies the admissible state set $\mathcal G$ of relativistic magnetohydrodynamics (RMHD). It paves a way for developing physical-constraints-preserving (PCP) schemes for the RMHD equations with the solutions in $\mathcal G$. To overcome the difficulties arising from the extremely strong nonlinearities and no explicit formulas of the primitive variables and the flux vectors with respect to the conservative vector, two equivalent forms of $\mathcal G$ with explicit constraints on the conservative vector are skillfully discovered. The first is derived by analyzing roots of several polynomials and transferring successively them, and further used to prove the convexity of $\mathcal G$ with the aid of semi-positive definiteness of the second fundamental form of a hypersurface. While the second is derived based on the convexity, and then used to show the orthogonal invariance of $\mathcal G$. The Lax–Friedrichs (LxF) splitting property does not hold generally for the nonzero magnetic field, but by a constructive inequality and pivotal techniques, we discover the generalized LxF splitting properties, combining the convex combination of some LxF splitting terms with a discrete divergence-free condition of the magnetic field. Based on the above analyses, several 1D and 2D PCP schemes are then studied. In the 1D case, a first-order accurate LxF-type scheme is first proved to be PCP under the Courant–Friedrichs–Lewy (CFL) condition, and then the high-order accurate PCP schemes are proposed via a PCP limiter. In the 2D case, the discrete divergence-free condition and PCP property are analyzed for a first-order accurate LxF-type scheme, and two sufficient conditions are derived for high-order accurate PCP schemes. Our analysis reveals in theory for the first time that the discrete divergence-free condition is closely connected with the PCP property. Several numerical examples demonstrate the theoretical findings and the performance of numerical schemes.

A fourth order finite difference method is presented for the 2D unsteady viscous incompressible Boussinesq equations in vorticity-stream function formulation. The method is especially suitable for moderate to large Reynolds number flows. The momentum equation is discretized by a compact fourth order scheme with the no-slip boundary condition enforced using a local vorticity boundary condition. Fourth order long-stencil discretizations are used for the temperature transport equation with one-sided extrapolation applied near the boundary. The time stepping scheme for both equations is classical fourth order Runge–Kutta. The method is highly efficient. The main computation consists of the solution of two Poisson-like equations at each Runge–Kutta time stage for which standard FFT based fast Poisson solvers are used. An example of Lorenz flow is presented, in which the full fourth order accuracy is checked. The numerical simulation of a strong shear flow induced by a temperature jump, is resolved by two perfectly matching resolutions. Additionally, we present benchmark quality simulations of a differentially-heated cavity problem. This flow was the focus of a special session at the first MIT conference on Computational Fluid and Solid Mechanics in June 2001.

No abstract uploaded!