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.

No abstract uploaded!

In this paper, we give an integral approximation for the elliptic operators with anisotropic coecients on smooth manifold. Using the integral approximation, the elliptic equation is tranformed to an integral equation. The integral approximation preserves the symmetry and coercivity of the original elliptic operator. Based on these good properties, we prove the convergence between the solutions of the integral equation and the original elliptic equation.

Abstract BDDC and FETI-DP algorithms are developed for three-dimensional elliptic problems with adaptively enriched coarse components. It is known that these enriched components are necessary in the development of robust preconditioners. To form the adaptive coarse components, carefully designed generalized eigenvalue problems are introduced for each faces and edges, and the coarse components are formed by using eigenvectors with their corresponding eigenvalues larger than a given tolerance T O L. Upper bounds for condition numbers of the preconditioned systems are shown to be C T O L, with the constant C depending only on the maximum number of edges and faces per subdomain, and the maximum number of subdomains sharing an edge. Numerical results are presented to test the robustness of the proposed approach.

In this paper, we propose a simple bound-preserving sweeping procedure for conservative numerical approximations. Conservative schemes are of importance in many applications, yet for high order methods, the numerical solutions do not necessarily satisfy maximum principle. This paper constructs a simple sweeping algorithm to enforce the bound of the solutions. It has a very general framework acting as a postprocessing step accommodating many point-based or cell average-based discretizations. The method is proven to preserve the bounds of the numerical solution while conserving a prescribed quantity designated as a weighted average of values from all points. The technique is demonstrated to work well with a spectral method, high order finite difference and finite volume methods for scalar conservation laws and incompressible flows. Extensive numerical tests in 1D and 2D are provided to verify the accuracy of the sweeping procedure.