  Inspired by the numerical evidence of a potential 3D Euler singularity by
Luo-Hou [30,31] and the recent breakthrough by Elgindi [11] on the singularity
formation of the 3D Euler equation without swirl with $C^{1,\alpha}$ initial
velocity, we prove the finite time singularity for the 2D Boussinesq and the 3D
axisymmetric Euler equations with boundary and $C^{1,\alpha}$ initial data for
the velocity (and density in the case of Boussinesq equations). Our finite time
blowup solution for the 3D Euler equations and the singular solution considered
in [30,31] share many essential features, including the symmetry properties of
the solution, the flow structure, and the sign of the solution in each
quadrant, except that we use $C^{1,\alpha}$ initial velocity. We use a dynamic
rescaling formulation and follow the general framework of analysis developed in
[11]. We also use some strategy proposed in our recent joint work with Huang in
[7] and adopt several methods of analysis in [11] to establish the linear and
nonlinear stability of an approximate self-similar profile. The nonlinear
stability enables us to prove that the solution of the 3D Euler or the 2D
Boussinesq equations with $C^{1,\alpha}$ initial data will develop a finite
time singularity. Moreover, the velocity field has finite energy before the
singularity time.
  In the previous version of this paper, we proved the blowup results for the
3D axisymmetric Euler equations with initial data $(u_0^{\theta})^2, u_0^r,
u_0^z \in C^{1,\alpha}$. Though the velocity $u^r, u^z$ in the axisymmetric
setting is $C^{1,\alpha}$, our interpretation that the velocity is
$C^{1,\alpha}$ is not correct since the velocity in 3D also depends on
$u^{\theta}$, which is not $C^{1,\alpha}$. This oversight can be fixed easily
with a minor change in the construction of the approximate steady state and a
minor modification to localize the approximate steady state.
