No abstract uploaded!

Kurdyka–Łojasiewicz (KL) exponent plays an important role in estimating the convergence rate of many contemporary first-order methods. In particular, a KL exponent of 1/2 for a suitable potential function is related to local linear convergence. Nevertheless, KL exponent is in general extremely hard to estimate. In this paper, we show under mild assumptions that KL exponent is preserved via inf-projection. Inf-projection is a fundamental operation that is ubiquitous when reformulating optimization problems via the lift-and-project approach. By studying its operation on KL exponent, we show that the KL exponent is 1/2 for several important convex optimization models, including some semidefinite-programming-representable functions and some functions that involve C2-cone reducible structures, under conditions such as strict complementarity. Our results are applicable to concrete optimization models such as group-fused Lasso and overlapping group Lasso. In addition, for nonconvex models, we show that the KL exponent of many difference-of-convex functions can be derived from that of their natural majorant functions, and the KL exponent of the Bregman envelope of a function is the same as that of the function itself. Finally, we estimate the KL exponent of the sum of the least squares function and the indicator function of the set of matrices of rank at most k.

Whether the three-dimensional incompressible Euler equations can develop a singularity in finite time from smooth initial data is one of the most challenging problems in mathematical fluid dynamics. This work attempts to provide an affirmative answer to this long-standing open question from a numerical point of view by presenting a class of potentially singular solutions to the Euler equations computed in axisymmetric geometries. The solutions satisfy a periodic boundary condition along the axial direction and a no-flow boundary condition on the solid wall. The equations are discretized in space using a hybrid 6th-order Galerkin and 6th-order finite difference method on specially designed adaptive (moving) meshes that are dynamically adjusted to the evolving solutions. With a maximum effective resolution of over $(3 × 10^{12})^2$ near the point of the singularity, we are able to advance the solution up to $\tau_2 = 0.003505$ and predict a singularity time of $t_s ≈ 0.0035056$, while achieving a pointwise relative error of $O(10^{−4})$ in the vorticity vector ω and observing a $(3 × 10^8)$-fold increase in the maximum vorticity $\|ω\|_{\infty}$. The numerical data are checked against all major blowup/non-blowup criteria, including Beale–Kato–Majda, Constantin–Fefferman–Majda, and Deng–Hou–Yu, to confirm the validity of the singularity. A local analysis near the point of the singularity also suggests the existence of a self-similar blowup in the meridian plane.

In this paper, we propose a new type of weighted essentially non-oscillatory (WENO) limiter, which belongs to the class of Hermite WENO (HWENO) limiters, for the Runge-Kutta discontinuous Galerkin (RKDG) methods solving hyperbolic conservation laws. This new HWENO limiter is a modification of the simple WENO limiter proposed recently by Zhong and Shu. Both limiters use information of the DG solutions only from the target cell and its immediate neighboring cells, thus maintaining the original compactness of the DG scheme. The goal of both limiters is to obtain high order accuracy and non-oscillatory properties simultaneously. The main novelty of the new HWENO limiter in this paper is to reconstruct the polynomial on the target cell in a least square fashion while the simple WENO limiter is to use the entire polynomial of the original DG solutions in the neighboring cells with an addition of a constant for conservation. The modification in this paper improves the robustness in the computation of problems with strong shocks or contact discontinuities, without changing the compact stencil of the DG scheme. Numerical results for both one and two dimensional equations including Euler equations of compressible gas dynamics are provided to illustrate the viability of this modified limiter.

Discontinuous Galerkin (DG) methods combine features in finite element methods (weak formulation, finite dimensional solution and test function spaces) and in finite volume methods (numerical fluxes, nonlinear limiters) and are particularly suitable for simulating convection dominated problems, such as linear and nonlinear waves including shock waves. In this article we will give a brief survey of DG methods, emphasizing their applications in computational fluid dynamics (CFD). We will discuss essential ingredients and properties of DG methods, and will also give a few examples of recent developments of DG methods which have facilitated their applications in CFD.