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 structured bisection method with adaptive randomized sampling for finding selected or all of the eigenvalues of certain real symmetric matrices A. For A with a low-rank property, we construct a hierarchically semiseparable (HSS) approximation and show how to quickly evaluate and update its inertia in the bisection method. Unlike some existing randomized HSS constructions, the methods here do not require the knowledge of the off-diagonal (numerical) ranks in advance. Moreover, for A with a weak rank property or slowly decaying off-diagonal singular values, we show an idea for aggressive low-rank inertia evaluation, which means that a compact HSS approximation can preserve the inertia for certain shifts. This is analytically justified for a special case, and numerically shown for more general ones. A generalized LDL factorization of the HSS approximation is then designed for the fast evaluation of the inertia. A significant advantage over standard LDL factorizations is that the HSS LDL factorization (and thus the inertia) of A − sI can be quickly updated with multiple shifts s in bisection. The factorization with each new shift can reuse about 60% of the work. As an important application, the structured eigensolver can be applied to symmetric Toeplitz matrices, and the cost to find one eigenvalue is nearly linear in the order of the matrix. The numerical examples demonstrate the efficiency and the accuracy of our methods, especially the benefit of low-rank inertia evaluations. The ideas and methods can be potentially adapted to other HSS computations where shifts are involved and to more problems without a significant low-rank property.

Partial differential equations (PDE) on manifolds arise in many areas, including mathematics and many applied fields. Due to the complicated geometrical structure of the manifold, it is difficult to get efficient numerical method to solve PDE on manifold. In the paper, we propose a method called point integral method (PIM) to solve the Poisson-type equations from point clouds. Among different kinds of PDEs, the Poisson-type equations including the standard Poisson equation and the related eigenproblem of the Laplace-Beltrami operator are one of the most important. In PIM, the key idea is to derive the integral equations which approximates the Poisson-type equations and contains no derivatives but only the values of the unknown function. This feature makes the integral equation easy to be discretized from point cloud. In the paper, we explain the derivation of the integral equations, describe the point integral method and its implementation, and present the numerical experiments to demonstrate the convergence of PIM.

In this paper, we propose and analyze a new stabilized finite element method using continuous piecewise linear (or bilinear) elements for solving 2D reaction-convection-diffusion equations. The equation under consideration involves a small diffusivity $\varepsilon$ and a large reaction coefficient $\sigma$, leading to high P\'eclet number and high Damk\"{o}hler number. In addition to giving error estimates of the approximations in $L^2$ and $H^1$ norms, we explicitly establish the dependence of error bounds on the diffusivity, the $L^{\infty}$ norm of convection field, the reaction coefficient and the mesh size. Our analysis shows that the proposed method is particularly suitable for problems with a small diffusivity and a large reaction coefficient, or more precisely, with a large mesh P\'eclet number and a large mesh Damk\"{o}hler number. Several numerical examples exhibiting boundary or interior layers are given to illustrate the high accuracy and stability of the proposed method. The results obtained are also compared with those of existing stabilization methods.

In this paper, we propose and analyze a novel stabilized finite element method (FEM) for the system of generalized Stokes equations arising from the time-discretization of transient Stokes problem. The system involves a small viscosity, which is proportional to the inverse of large Reynolds number, and a large reaction coefficient, which is the inverse of small time step. The proposed stabilized FEM employs the $C^0$ piecewise linear elements for both velocity field and pressure on the same mesh and uses the residuals of the momentum equation and the divergence-free equation to define the stabilization terms. The stabilization parameters are fixed and element-independent, without a comparison of the viscosity, the reaction coefficient and the mesh size. Using the finite element solution of an auxiliary boundary value problem as the interpolating function for velocity and the $H^1$-seminorm projection for pressure, instead of the usual nodal interpolants, we derive error estimates for the stabilized finite element approximations to velocity and pressure in the $L^2$ and $H^1$ norms and most importantly, we explicitly establish the dependence of error bounds on the viscosity, the reaction coefficient and the mesh size. Our analysis reveals that this stabilized FEM is particularly suitable for the generalized Stokes system with a small viscosity and a large reaction coefficient, which has never been achieved before in the error analysis of other stabilization methods in the literature. We then numerically confirm the effectiveness of the proposed stabilized FEM. Comparisons made with other existing stabilization methods show that the newly proposed method can attain better accuracy and stability.