We study analytical and numerical properties of the L 1 - L 2 minimization problem for sparse representation of a signal over a highly coherent dictionary. Though the L 1 - L 2 metric is non-convex, it is Lipschitz continuous. The difference of convex algorithm (DCA) is readily applicable for computing the sparse representation coefficients. The L 1 - L 2 minimization appears as an initialization step of DCA. We further integrate DCA with a non-standard simulated annealing methodology to approximate globally sparse solutions. Non-Gaussian random perturbations are more effective than standard Gaussian perturbations for improving sparsity of solutions. In numerical experiments, we conduct an extensive comparison among sparse penalties such as L 1 - L 2 for L 1 - L 2 based on data from three specific applications (over-sampled discreet cosine basis, differential absorption optical spectroscopy, and image denoising) where

How to properly specify boundary conditions for pressure is a longstanding problem for the incompressible Navier–Stokes equations with no-slip boundary conditions. An analytical resolution of this issue stems from a recently developed formula for the pressure in terms of the commutator of the Laplacian and Leray projection operators. Here we make use of this formula to (a) improve the accuracy of computing pressure in two kinds of existing time-discrete projection methods implicit in viscosity only, and (b) devise new higher-order accurate time-discrete projection methods that extend a slip-correction idea behind the well-known finite-difference scheme of Kim and Moin. We test these schemes for stability and accuracy using various combinations of C0 finite elements. For all three kinds of time discretization, one can obtain third-order accuracy for both pressure and velocity without a time-step stability restriction of diffusive type. Furthermore, two kinds of projection methods are found stable using piecewise-linear elements for both velocity and pressure.

No abstract uploaded!

In this paper, we design a new class of central compact schemes based on the cell-centered compact schemes of Lele {[}S.K. Lele, Compact finite difference schemes with spectral-like resolution, J. Comput. Phys. 103 (1992) 16-42{]}. These schemes equate a weighted sum of the nodal derivatives of a smooth function to a weighted sum of the function on both the grid points (cell boundaries) and the cell-centers. In our approach, instead of using a compact interpolation to compute the values on cell-centers, the physical values on these half grid points are stored as independent variables and updated using the same scheme as the physical values on the grid points. This approach increases the memory requirement but not the computational costs. Through systematic Fourier analysis and numerical tests, we observe that the schemes have excellent properties of high order, high resolution and low dissipation. It is an ideal class of schemes for the simulation of multi-scale problems such as aeroacoustics and turbulence.

In this paper we propose a finite element time-domain method for modeling the optical black holes (OBHs) coupled with the perfectly matched layer (PML) technique. Stability analysis is carried out for the proposed scheme. Simulations of cylindrical, elliptical and square black holes demonstrate that our method is quite effective in modeling OBHs in time domain. To our best knowledge, this is the first OBHs simulation realized by the finite element time-domain method.