Consider the image restoration from incomplete noisy frequency data with total variation and sparsity regularizing penalty terms. Firstly, we establish an unconstrained optimization model with dierent smooth approximations on the regularizing terms. Then, to weaken the amount of computations for cost functional with total variation term, the alternating iterative scheme is developed to obtain the exact solution through shrinkage thresholding in inner loop, while the nonlinear Euler equation is appropriately linearized at each iteration in exterior loop, yielding a linear system with diagonal coefficient matrix in frequency domain. Finally the linearized iteration is proven to be convergent in generalized sense for suitable regularizing parameters, and the error between the linearized iterative solution and the one gotten from the exact nonlinear Euler equation is rigorously estimated, revealing the essence of the proposed alternative iteration scheme. Numerical tests for different configurations show the validity of the proposed scheme, compared with some existing algorithms.

In this article, we give a brief overview on high order accurate shock capturing schemes with the aim of applications in compressible turbulence simulations. The emphasis is on the basic methodology and recent algorithm developments for two classes of high order methods: the weighted essentially non-oscillatory (WENO) and discontinuous Galerkin (DG) methods.

In this paper, we investigate the superconvergence property of the local discontinuous Galerkin (LDG) methods for solving one-dimensional linear time dependent fourth order problems. We prove that the error between the LDG solution and a particular projection of the exact solution, $\bar{e}_u$, achieves $(k+\frac32)$-th order superconvergence when polynomials of degree $k$ ($k\ge 1$) are used. Numerical experiments of $P^k$ polynomials, with $1\le k \le3$, are displayed to demonstrate the theoretical results, which show that the error $\bar{e}_u$ actually achieves $(k+2)$-th order superconvergence, indicating that the error bound for $\bar{e}_u$ obtained in this paper is sub-optimal. Initial-boundary value problems, nonlinear equations and solutions having singularities are numerically investigated to verify that the conclusions hold true for very general cases.

A spectral collocation scheme for the three-dimensional incompressible \(({\varvec{u}},p)\) formulation of the Navier鈥揝tokes equations, in domains \(\varOmega \) with a non-periodic boundary condition, is described. The key feature is the high order approximation, by means of a local Hermite interpolant, of a Neumann boundary condition for use in the numerical solution of the pressure Poisson system. The time updates of the velocity \({\varvec{u}}\) and pressure \(p\) are decoupled as a result of treating the pressure gradient in the momentum equation explicitly in time. The pressure update is computed from a pressure Poisson equation. Extension of the overall methodology to the Boussinesq system is also described. The uncoupling of the pressure and velocity time updates results in a highly efficient scheme that is simple to implement and well suited for simulating moderate to high Reynolds and Rayleigh number flows. Accuracy checks are presented, along with simulations of the lid-driven cavity flow and a differentially heated cavity flow, to demonstrate the scheme produces accurate three-dimensional results at a reasonable computational cost.

We consider an incompressible viscous flow without surface tension in a finite-depth domain of two dimensions, with free top boundary and fixed bottom boundary. This system is governed by the Navier-Stokes equations in this moving domain and the transport equation on the moving boundary. In this paper, we construct a stable numerical scheme to simulate the evolution of this system by discontinuous Galerkin method, and discuss the error analysis of the fluid under certain assumptions. Our formulation is mainly based on the geometric structure introduced earlier in the literature, and the natural energy estimate, which is rarely used in the numerical study of this system before.