The sequential regularization method is a reformulation of the unsteady Navier–Stokes equations from the viewpoint of constrained dynamical systems or the approximate Helmholtz–Hodge projection. In this paper we study the long time behavior of the sequential regularization formulation. We give a uniform-in-time estimate between the solution of the reformulated system and that of the Navier–Stokes equations. We also conduct an error analysis for the temporal discrete system and show that the error bound is independent of time. A couple of long time flow examples are computed to demonstrate this method.

No abstract uploaded!

This paper analyzes the random fluctuations obtained by a heterogeneous multi-scale first-order finite element method applied to solve elliptic equations with a random potential. Several multiscale numerical algorithms have been shown to correctly capture the homogenized limit of solutions of elliptic equations with coefficients modeled as stationary and ergodic random fields. Because theoretical results are available in the continuum setting for such equations, we consider here the case of a secondorder elliptic equations with random potential in two dimensions of space. We show that the random fluctuations of such solutions are correctly estimated by the heterogeneous multi-scale algorithm when appropriate fine-scale problems are solved on subsets that cover the whole computational domain. However, when the fine-scale problems are solved over patches that do not cover the entire domain, the random fluctuations may or may not be estimated accurately. In the case of random potentials with short-range interactions, the variance of the random fluctuations is amplified as the inverse of the fraction of the medium covered by the patches. In the case of random potentials with long-range interactions, however, such an amplification does not occur and random fluctuations are correctly captured independent of the (macroscopic) size of the patches. These results are consistent with those obtained in [9] for more general equations in the one-dimensional setting and provide indications on the loss in accuracy that results from using coarser, and hence computationally less intensive, algorithms.

In this paper, high order central finite difference schemes in a finite interval are analyzed for the diffusion equation. Boundary conditions of the initial-boundary value problem (IBVP) are treated by the simplified inverse Lax-Wendroff (SILW) procedure. For the fully discrete case, a third order explicit Runge-Kutta method is used as an example for the analysis. Stability is analyzed by both the GKS (Gustafsson, Kreiss and Sundstr\"om) theory and the eigenvalue visualization method on both semi-discrete and fully discrete schemes. The two different analysis techniques yield consistent results. Numerical tests are performed to demonstrate and validate the analysis results.

In this paper, we develop a three-dimensional parallel solver using the fifth order high-resolution weighted essentially non-oscillatory (WENO) finite difference scheme to perform extensive simulation for three-dimensional gaseousdetonations. A careful study is conducted for the propagation modes of three-dimensional gaseous detonation wave-front structures in a long square duct. The numerical results indicate that, the instability of detonation, overdrive factor, transverse dimension and initial perturbation are the main factors that influence the appearance of spinning detonation.