A novel discontinuous Galerkin (DG) method is developed to solve timedependent bi-harmonic type equations involving fourth derivatives in one and multiple space dimensions. We present the spatial DG discretization based on a mixed formulation and central interface numerical fluxes so that the resulting semi-discrete schemes are L^2 stable even without interior penalty. For time discretization, we use Crank–Nicolson so that the resulting scheme is unconditionally stable and second order in time. We present the optimal L^2 error estimate of O(h^{k+1}) for polynomials of degree k for semi-discrete DG schemes, and the L2 error of O(h^{k+1} + (\Delta t)^2) for fully discrete DG schemes. Extensions to more general fourth order partial differential equations and cases with non-homogeneous boundary conditions are provided. Numerical results are presented to verify the stability and accuracy of the schemes. Finally, an application to the one-dimensional Swift–Hohenberg equation endowed with a decay free energy is presented.

In this paper we discuss the derivation and use of local pressure boundary conditions for finite difference schemes for the unsteady incompressible Navier-Stokes equations in the velocity-pressure formulation. Their use is especially well suited for the computation of moderate to large Reynolds number flows. We explore the similarities between the implementation and use of local pressure boundary conditions and local vorticity boundary conditions in the design of numerical schemes for incompressible flow in 2D. In their respective formulations, when these local numerical boundary conditions are coupled with a fully explicit convectively stable time stepping procedure, the resulting methods are simple to implement and highly efficient. Unlike the vorticity formulation, the use of the local pressure boundary condition approach is readily applicable to 3D flows. The simplicity of the local pressure boundary condition approach and its easy adaptation to more general flow settings make the resulting scheme an attractive alternative to the more popular methods for solving the Navier-Stokes equations in the velocity-pressure formulation. We present numerical results of a second-order finite difference scheme on a nonstaggered grid using local pressure boundary conditions. Stability and accuracy of the scheme applied to Stokes flow is demonstrated using normal mode analysis. Also described is the extension of the method to variable density flows.

In this paper, we apply a simple finite element numerical scheme, proposed in an earlier work (Liu in Math Comput 70(234):579-593, 2000), to perform a high resolution numerical simulation of incompressible flow over an irregular domain and analyze its boundary layer separation. Compared with many classical finite element fluid solvers, this numerical method avoids a Stokes solver, and only two Poisson-like equations need to be solved at each time step/stage. In addition, its combination with the fully explicit fourth order Runge-Kutta (RK4) time discretization enables us to compute high Reynolds number flow in a very efficient way. As an application of this robust numerical solver, the dynamical mechanism of the boundary layer separation for a triangular cavity flow with Reynolds numbers \(Re=10^4\) and \(Re=10^5\) , including the precise values of bifurcation location and critical time, are reported in this paper. In addition, we provide a super-convergence analysis for the simple finite element numerical scheme, using linear elements over a uniform triangulation with right triangles.

A finite element method for computing viscous incompressible flows based on the gauge formulation introduced in [Weinan E, Liu J-G. Gauge method for viscous incompressible flows. Journal of Computational Physics (submitted)] is presented. This formulation replaces the pressure by a gauge variable. This new gauge variable is a numerical tool and differs from the standard gauge variable that arises from decomposing a compressible velocity field. It has the advantage that an additional boundary condition can be assigned to the gauge variable, thus eliminating the issue of a pressure boundary condition associated with the original primitive variable formulation. The computational task is then reduced to solving standard heat and Poisson equations, which are approximated by straightforward, piecewise linear (or higher-order) finite elements. This method can achieve high-order accuracy at a cost comparable with that of solving standard heat and Poisson equations. It is naturally adapted to complex geometry and it is much simpler than traditional finite element methods for incompressible flows. Several numerical examples on both structured and unstructured grids are presented. Copyright

We observe that Sturm’s error bounds readily imply that for semidefinite feasibility problems, the method of alternating projections converges at a rate of O(k^{-1/(2^{d+1}-2)}), where d is the singularity degree of the problem—the minimal number of facial reduction iterations needed to induce Slater’s condition. Consequently, for almost all such problems (in the sense of Lebesgue measure), alternating projections converge at a worst-case rate of O(1/k^{0.5}).