We propose a numerically efficient algorithm for simulating the multi-color optical self-focusing phenomena in nematic liquid crystals. The propagation of the nematicon is modeled by a parabolic wave equation coupled with a nonlinear elliptic partial differential equation governing the angle between the crystal and the direction of propagation. Numerically, the paraxial parabolic wave equation is solved by a fast Huygens sweeping method, while the nonlinear elliptic PDE is handled by the alternating direction explicit (ADE) method. The overall algorithm is shown to be numerically efficient for computing high frequency beam propagations.

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.

In this paper, we adapt a simple weighted essentially non-oscillatory (WENO) limiter, originally designed for discontinuous Galerkin (DG) schemes, to the correction procedure via reconstruction (CPR) framework for solving conservation laws. The objective of this simple WENO limiter is to simultaneously maintain uniform high order accuracy of the CPR framework in smooth regions and control spurious numerical oscillations near discontinuities. The WENO limiter we adopt in this paper is particularly simple to implement and will not harm the conservativeness of the CPR framework. Also, it uses information only from the target cell and its immediate neighbors, thus can maintain the compactness of the CPR framework. Since the CPR framework with the WENO limiter does not in general preserve positivity of the solution, we also extend the positivity preserving limiters to the CPR framework. Numerical results in one and two dimensions are provided to illustrate the good behavior of this procedure.

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.

The authors establish error estimates for recently developed finite-element methods for incompressible viscous flow in domains with no-slip boundary conditions. The methods arise by discretization of a well-posed extended Navier-Stokes dynamics for which pressure is determined from current velocity and force fields. The methods use C1 elements for velocity and C0 elements for pressure. A stability estimate is proved for a related finite-element projection method close to classical time-splitting methods of Orszag, Israeli, DeVille and Karniadakis.