Search

We present a novel method of analysis and prove finite time asymptotically self-similar blowup of the De Gregorio model \cite{DG90,DG96} for some smooth initial data on the real line with compact support. We also prove self-similar blowup results for the generalized De Gregorio model \cite{OSW08} for the entire range of parameter on $\mathbb{R}$ or $S^1$ for H\"older continuous initial data with compact support. Our strategy is to reformulate the problem of proving finite time asymptotically self-similar singularity into the problem of establishing the nonlinear stability of an approximate self-similar profile with a small residual error using the dynamic rescaling equation. We use the energy method with appropriate singular weight functions to extract the damping effect from the linearized operator around the approximate self-similar profile and take into account cancellation among various nonlocal terms to establish stability analysis. We remark that our analysis does not rule out the possibility that the original De Gregorio model is well posed for smooth initial data on a circle. The method of analysis presented in this paper provides a promising new framework to analyze finite time singularity of nonlinear nonlocal systems of partial differential equations., Comment: Added discussion in Section 2.3 and made some minor edits. Main paper 57 pages, Supplementary material 29 pages. In previous arXiv versions, the hyperlinks of the equation number in the main paper are linked to the supplementary material, which is fixed in this version

In connection with the recent proposal for possible singularity formation at the boundary for solutions of three-dimensional axisymmetric incompressible Euler's equations (Luo and Hou, Proc. Natl. Acad. Sci. USA (2014)), we study models for the dynamics at the boundary and show that they exhibit a finite-time blowup from smooth data.

We investigate the stabilizing effect of convection in three-dimensional incompressible Euler and Navier-Stokes equations. The convection term is the main source of nonlinearity for these equations. It is often considered destabilizing although it conserves energy due to the incompressibility condition. In this paper, we show that the convection term together with the incompressibility condition actually has a surprising stabilizing effect. We demonstrate this by constructing a new three-dimensional model that is derived for axisymmetric flows with swirl using a set of new variables. This model preserves almost all the properties of the full three-dimensional Euler or Navier-Stokes equations except for the convection term, which is neglected in our model. If we added the convection term back to our model, we would recover the full Navier-Stokes equations. We will present numerical evidence that seems to support that the three-dimensional model may develop a potential finite time singularity. We will also analyze the mechanism that leads to these singular events in the new three-dimensional model and how the convection term in the full Euler and Navier-Stokes equations destroys such a mechanism, thus preventing the singularity from forming in a finite time.

In this paper, we study the dynamic stability of the three-dimensional axisymmetric Navier-Stokes Equations with swirl. To this purpose, we propose a new one-dimensional model that approximates the Navier-Stokes equations along the symmetry axis. An important property of this one-dimensional model is that one can construct from its solutions a family of exact solutions of the three-dimensionaFinal Navier-Stokes equations. The nonlinear structure of the one-dimensional model has some very interesting properties. On one hand, it can lead to tremendous dynamic growth of the solution within a short time. On the other hand, it has a surprising dynamic depletion mechanism that prevents the solution from blowing up in finite time. By exploiting this special nonlinear structure, we prove the global regularity of the three-dimensional Navier-Stokes equations for a family of initial data, whose solutions can lead to large dynamic growth, but yet have global smooth solutions.

We consider the motion of a two-dimensional interface separating an inviscid, incompressible, irrotational fluid, influenced by gravity, from a region of zero density. We show that under certain conditions the equations of motion, linearized about a presumed time-dependent solution, are wellposed; that is, linear disturbances have a bounded rate of growth. If surface tension is neglected, the linear equations are well-posed provided the underlying exact motion satisfies a condition on the acceleration of the interface relative to gravity, similar to the criterion formulated by G. I. Taylor. If surface tension is included, the linear equations are well-posed without qualifications, whether the fluid is above or below the interface. An interesting qualitative structure is found for the linear equations. A Lagrangian approach is used, like that of numerical work such as [3], except that the interface is assumed horizontal at infinity. Certain integral equations which occur, involving double layer potentials, are shown to be solvable in the present case. © 1993 John Wiley & Sons, Inc.

Multi-valued solutions are constructed for 2 × 2 first-order systems using a generalization of the hodograph transformation. The solution is found as a complex analytic function on a complex Riemann surface for which the branch points move as part of the solution. The branch point singularities are envelopes for the characteristics and thus move at the characteristic speeds. We perform an analysis of stability of these singularities with respect to perturbations of the initial data. The generic singularity types are folds, cusps, and nondegenerate umbilic points with non-zero 3-jet. An isolated singularity is generically a square root branch point corresponding to a fold. Two types of collisions between singularities are generic: At a “tangential” collision between two singularities moving at the same characteristic speed, a cube root branch point is formed, corresponding to a cusp. A “non-tangential” collision, between two square root branch points moving at different characteristic speeds, remains a square root branch point at the collision and corresponds to a nondegenerate umbilic point. These results are also valid for a diagonalizable n-th order system for which there are exactly two speeds. © 1993 John Wiley & Sons, Inc.

We have compared the results of a von Neumann-Richtmyer calculation with the weak limit of calculations performed by von Neumann's original method; although there are discernible differences, they are surprisingly small. We analyze von Newmann's algorithm for solving the initial value problem for the Lagrangian equations of compressible flow

We prove consistency, stability, and convergence of a point vortex approximation to the 3-D incompressible Euler equations with smooth solutions. The 3-D algorithm we consider here is similar to the corresponding 3-D vortex blob algorithm introduced by Beale and Majda; see [3]. We first show that the discretization error is second-order accurate. Then we show that the method is stable in lp norm for the particle trajectories and in w−1.p norm for discrete vorticity. Consequently, the method converges up to any time for which the Euler equations have a smooth solution. One immediate application of our convergence result is that the vortex filament method without smoothing also converges.

The 2-D incompressible Euler equations with oscillatory vorticity fields are studied. A homogenization result for 2-D Euler equations in velocity-vorticity formulation is obtained and weak continuity of the equations is proved. Convergence of the vortex method is analyzed in the case when the continuous vorticity is not well resolved by the computational particles. Numerical results are given. Comparisons are made with the corresponding finite difference approximation.

A grid free Lagrangian finite element method is introduced for the 2-D incompressible Euler equations. The method is derived based on the observation that the product of the Biot-Savart kernel and a polynomial can be integrated analytically over any triangle. This enables us to obtain a numerically stable Lagrangian method without using numerical smoothing. Moreover, we show that the method converges uniformly with second-order accuracy. Actually, we establish a l∞ stability result which applies to kernels that are more singular than the Biot-Savart kernel, as long as the kernel is L integrable. Another useful result is that we prove convergence of our method when using local regridding, which allows the method to run for longer time even with a fixed mesh. The second-order convergence is also illustrated by our numerical experiments.

We prove consistency, stability and convergence of the point vortex approximation to the 2-D incompressible Euler equations with smooth solutions. We first show that the discretization error is second-order accurate. Then we show that the method is stable in l p norm. Consequently the method converge in l p norm for all time. The convergence is also illustrated by a numerical experiment

The behavior of multi-dimensional discrete Boltzmann systems with highly oscillatory data is studied. Homogenized equations for the mean solutions are obtained. Uniform convergence of the oscillatory solutions of the discrete Boltzmann equations to the solutions of the corresponding homogenized equations is established. Moreover, we find that the weak limits of the oscillatory solutions for a model of Broadwell type are not continuous functions of the discrete velocities. Generalization of the above results to problems with multiple-scale initial data is also established. Reprints.