Source author record

Thomas Y. Hou

Thomas Y. Hou appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

33works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2026arXiv

$L^2$-based stability of blowup with log correction for semilinear heat equation

We propose an alternative proof of the classical result of Type-I blowup with log correction for the semilinear heat equation. Compared with previous proofs, we use a novel idea of enforcing stable normalizations for perturbations around the approximate profile and we establish a weighted $H^k$ stability, thereby avoiding the use of a topological argument and the analysis of a linearized spectrum. Consequently, this approach can be adopted even if we only have a numerical profile and do not have explicit information on the spectrum of its linearized operator. This result generalizes the $L^2$-based stability framework beyond exactly self-similar blowup and can be adapted to higher dimensions. Numerical results corroborate the effectiveness of our normalization, even in the large perturbation regime beyond our theoretical setting.

preprint2022arXiv

Asymptotic Escape of Spurious Critical Points on the Low-rank Matrix Manifold

We show that on the manifold of fixed-rank and symmetric positive semi-definite matrices, the Riemannian gradient descent algorithm almost surely escapes some spurious critical points on the boundary of the manifold. Our result is the first to partially overcome the incompleteness of the low-rank matrix manifold without changing the vanilla Riemannian gradient descent algorithm. The spurious critical points are some rank-deficient matrices that capture only part of the eigen components of the ground truth. Unlike classical strict saddle points, they exhibit very singular behavior. We show that using the dynamical low-rank approximation and a rescaled gradient flow, some of the spurious critical points can be converted to classical strict saddle points in the parameterized domain, which leads to the desired result. Numerical experiments are provided to support our theoretical findings.

preprint2022arXiv

Finite time blowup of 2D Boussinesq and 3D Euler equations with $C^{1,α}$ velocity and boundary

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,α}$ initial velocity, we prove the finite time singularity for the 2D Boussinesq and the 3D axisymmetric Euler equations with boundary and $C^{1,α}$ 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,α}$ 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,α}$ 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^θ)^2, u_0^r, u_0^z \in C^{1,α}$. Though the velocity $u^r, u^z$ in the axisymmetric setting is $C^{1,α}$, our interpretation that the velocity is $C^{1,α}$ is not correct since the velocity in 3D also depends on $u^θ$, which is not $C^{1,α}$. 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.

preprint2022arXiv

On stability and instability of $C^{1,α}$ singular solutions to the 3D Euler and 2D Boussinesq equations

Singularity formation of the 3D incompressible Euler equations is known to be extremely challenging. In [18], Elgindi proved that the 3D axisymmetric Euler equations with no swirl and $C^{1,α}$ initial velocity develops a finite time singularity. Inspired by Elgindi's work, we proved that the 3D axisymmetric Euler and 2D Boussinesq equations with $C^{1,α}$ initial velocity and boundary develop a stable asymptotically (or approximately) self-similar finite time singularity [8]. On the other hand, the authors of [35,52] recently showed that blowup solutions to the 3D Euler equations are hydrodynamically unstable. The instability results obtained in [35,52] require some strong regularity assumption on the initial data, which is not satisfied by the $C^{1,α}$ velocity field. In this paper, we generalize the analysis of [8,18,35,52] to show that the blowup solutions of the 3D Euler and 2D Boussinesq equations with $C^{1,α}$ velocity are unstable under the notion of stability introduced in [35,52]. These two seemingly contradictory results reflect the difference of the two approaches in studying the stability of 3D Euler blowup solutions. The stability analysis of the blowup solution obtained in [8,18] is based on the stability of a dynamically rescaled blowup profile in space and time, which is nonlinear in nature. The linear stability analysis in [35,52] is performed by directly linearizing the 3D Euler equations around a blowup solution in the original variables. It does not take into account the changes in the blowup time, the dynamic changes of the rescaling rate of the perturbed blowup profile and the blowup exponent of the original 3D Euler equations using a perturbed initial condition when there is an approximate self-similar blowup profile. Such information has been used in an essential way in establishing the nonlinear stability of the blowup profile in [8,18,19].

preprint2022arXiv

Potential singularity formation of incompressible axisymmetric Euler equations with degenerate viscosity coefficients

In this paper, we present strong numerical evidences that the incompressible axisymmetric Euler equations with degenerate viscosity coefficients and smooth initial data of finite energy develop a potential finite-time locally self-similar singularity at the origin. An important feature of this potential singularity is that the solution develops a two-scale traveling wave that travels towards the origin. The two-scale feature is characterized by the scaling property that the center of the traveling wave is located at a ring of radius $O((T-t)^{1/2})$ surrounding the symmetry axis while the thickness of the ring collapses at a rate $O(T-t)$. The driving mechanism for this potential singularity is due to an antisymmetric vortex dipole that generates a strong shearing layer in both the radial and axial velocity fields. Without the viscous regularization, the $3$D Euler equations develop a sharp front and some shearing instability in the far field. On the other hand, the Navier-Stokes equations with a constant viscosity coefficient regularize the two-scale solution structure and do not develop a finite-time singularity for the same initial data.

preprint2022arXiv

Potential singularity of the 3D Euler equations in the interior domain

Whether the 3D incompressible Euler equations can develop a finite time singularity from smooth initial data is one of the most challenging problems in nonlinear PDEs. In this paper, we present some new numerical evidence that the 3D axisymmetric incompressible Euler equations with smooth initial data of finite energy develop a potential finite time singularity at the origin. This potential singularity is different from the blow-up scenario revealed by Luo-Hou in \cite{luo2014potentially,luo2014toward}, which occurs on the boundary. Our initial condition has a simple form and shares several attractive features of a more sophisticated initial condition constructed by Hou-Huang in \cite{Hou-Huang-2021,Hou-Huang-2022}. One important difference between these two blow-up scenarios is that the solution for our initial data has a one-scale structure instead of a two-scale structure reported in \cite{Hou-Huang-2021,Hou-Huang-2022}. More importantly, the solution seems to develop nearly self-similar scaling properties that are compatible with those of the 3D Navier-Stokes equations. We will present numerical evidence that the 3D Euler equations seem to develop a potential finite time singularity. Moreover, the nearly self-similar profile seems to be very stable to the small perturbation of the initial data.

preprint2022arXiv

The potentially singular behavior of the 3D Navier-Stokes equations

Whether the 3D incompressible Navier-Stokes equations can develop a finite time singularity from smooth initial data is one of the most challenging problems in nonlinear PDEs. In this paper, we present some new numerical evidence that the incompressible axisymmetric Navier-Stokes equations with smooth initial data of finite energy seem to develop potentially singular behavior at the origin. This potentially singular behavior is induced by a potential finite time singularity of the 3D Euler equations that we reported in the companion paper (arXiv:2107.05870). We present numerical evidence that the 3D Navier--Stokes equations develop nearly self-similar singular scaling properties with maximum vorticity increased by a factor of 10^7. We have applied several blow-up criteria to study the potentially singular behavior of the Navier--Stokes equations. The Beale-Kato-Majda blow-up criterion and the blow-up criteria based on the growth of enstrophy and negative pressure seem to imply that the Navier--Stokes equations using our initial data develop a potential finite time singularity. We have also examined the Ladyzhenskaya-Prodi-Serrin regularity criteria. Our numerical results for the cases of (p,q) = (4,8), (6,4), (9,3) and (p,q)=(\infty,2) provide strong evidence for the potentially singular behavior of the Navier--Stokes equations. Our numerical study shows that while the global L^3 norm of the velocity grows very slowly, the localized version of the L^3 norm of the velocity experiences rapid dynamic growth relative to the localized L^3 norm of the initial velocity. This provides further evidence for the potentially singular behavior of the NavieStokes equations.

preprint2020arXiv

Analysis of Asymptotic Escape of Strict Saddle Sets in Manifold Optimization

In this paper, we provide some analysis on the asymptotic escape of strict saddles in manifold optimization using the projected gradient descent (PGD) algorithm. One of our main contributions is that we extend the current analysis to include non-isolated and possibly continuous saddle sets with complicated geometry. We prove that the PGD is able to escape strict critical submanifolds under certain conditions on the geometry and the distribution of the saddle point sets. We also show that the PGD may fail to escape strict saddles under weaker assumptions even if the saddle point set has zero measure and there is a uniform escape direction. We provide a counterexample to illustrate this important point. We apply this saddle analysis to the phase retrieval problem on the low-rank matrix manifold, prove that there are only a finite number of saddles, and they are strict saddles with high probability. We also show the potential application of our analysis for a broader range of manifold optimization problems.

preprint2020arXiv

Function Approximation via The Subsampled Poincar\' e Inequality

Function approximation and recovery via some sampled data have long been studied in a wide array of applied mathematics and statistics fields. Analytic tools, such as the Poincaré inequality, have been handy for estimating the approximation errors in different scales. The purpose of this paper is to study a generalized Poincar\' e inequality, where the measurement function is of subsampled type, with a small but non-zero lengthscale that will be made precise. Our analysis identifies this inequality as a basic tool for function recovery problems. We discuss and demonstrate the optimality of the inequality concerning the subsampled lengthscale, connecting it to existing results in the literature. In application to function approximation problems, the approximation accuracy using different basis functions and under different regularity assumptions is established by using the subsampled Poincaré inequality. We observe that the error bound blows up as the subsampled lengthscale approaches zero, due to the fact that the underlying function is not regular enough to have well-defined pointwise values. A weighted version of the Poincar\' e inequality is proposed to address this problem; its optimality is also discussed.

preprint2018arXiv

A Minimal Mechanosensing Model Predicts Keratocyte Evolution on Flexible Substrates

A mathematical model is proposed for shape evolution and locomotion of fish epidermal keratocytes on elastic substrates. The model is based on mechanosensing concepts: cells apply contractile forces onto the elastic substrate, while cell shape evolution depends locally on the substrate stress generated by themselves or external mechanical stimuli acting on the substrate. We use the level set method to study the behavior of the model numerically, and predict a number of distinct phenomena observed in experiments, such as (i) symmetry breaking from the stationary centrosymmetric to the well-known steadily propagating crescent shape, (ii) asymmetric bipedal oscillations and traveling waves in the lamellipodium leading edge (iii) response to mechanical stress externally applied to the substrate (tensotaxis), (iv) changing direction of motion towards an interface with a rigid substrate (durotaxis) and (v) the configuration of substrate wrinkles induced by contractile forces applied by the keratocyte.

preprint2016arXiv

A sparse decomposition of low rank symmetric positive semi-definite matrices

Suppose that $A \in \mathbb{R}^{N \times N}$ is symmetric positive semidefinite with rank $K \le N$. Our goal is to decompose $A$ into $K$ rank-one matrices $\sum_{k=1}^K g_k g_k^T$ where the modes $\{g_{k}\}_{k=1}^K$ are required to be as sparse as possible. In contrast to eigen decomposition, these sparse modes are not required to be orthogonal. Such a problem arises in random field parametrization where $A$ is the covariance function and is intractable to solve in general. In this paper, we partition the indices from 1 to $N$ into several patches and propose to quantify the sparseness of a vector by the number of patches on which it is nonzero, which is called patch-wise sparseness. Our aim is to find the decomposition which minimizes the total patch-wise sparseness of the decomposed modes. We propose a domain-decomposition type method, called intrinsic sparse mode decomposition (ISMD), which follows the "local-modes-construction + patching-up" procedure. The key step in the ISMD is to construct local pieces of the intrinsic sparse modes by a joint diagonalization problem. Thereafter a pivoted Cholesky decomposition is utilized to glue these local pieces together. Optimal sparse decomposition, consistency with different domain decomposition and robustness to small perturbation are proved under the so called regular-sparse assumption (see Definition 1.2). We provide simulation results to show the efficiency and robustness of the ISMD. We also compare the ISMD to other existing methods, e.g., eigen decomposition, pivoted Cholesky decomposition and convex relaxation of sparse principal component analysis [25] and [40].

preprint2016arXiv

Adaptive multiscale model reduction with Generalized Multiscale Finite Element Methods

In this paper, we discuss a general multiscale model reduction framework based on multiscale finite element methods. We give a brief overview of related multiscale methods. Due to page limitations, the overview focuses on a few related methods and is not intended to be comprehensive. We present a general adaptive multiscale model reduction framework, the Generalized Multiscale Finite Element Method. Besides the method's basic outline, we discuss some important ingredients needed for the method's success. We also discuss several applications. The proposed method allows performing local model reduction in the presence of high contrast and no scale separation.

preprint2016arXiv

Exploring the locally low dimensional structure in solving random elliptic PDEs

We propose a stochastic multiscale finite element method (StoMsFEM) to solve random elliptic partial differential equations with a high stochastic dimension. The key idea is to simultaneously upscale the stochastic solutions in the physical space for all random samples and explore the low stochastic dimensions of the stochastic solution within each local patch. We propose two effective methods to achieve this simultaneous local upscaling. The first method is a high order interpolation method in the stochastic space that explores the high regularity of the local upscaled quantities with respect to the random variables. The second method is a reduced-order method that explores the low rank property of the multiscale basis functions within each coarse grid patch. Our complexity analysis shows that compared with the standard FEM on a fine grid, the StoMsFEM can achieve computational saving in the order of $(H/h)^{d}/(\log(H/h))^k$, where $H/h$ is the ratio between the coarse and the fine gird sizes, $d$ is the physical dimension and $k$ is the local stochastic dimension. Several numerical examples are presented to demonstrate the accuracy and effectiveness of the proposed methods. In the high contrast example, we observe a factor of 2000 speed-up.

preprint2016arXiv

On the non-uniqueness of the instantaneous frequency

In this article, we investigate the debated Instantaneous Frequency (IF) topic. Here, we show that IF is non-unique inherently. We explain how this non-uniqueness can be quantified and explained from a mathematical perspective. The non-uniqueness of the IF can also be observed if different methods of adaptive signal processing are used. We will also show that even if we know the physical origin of an oscillatory signal, e.g. linear second order ordinary differential equation, the non-uniqueness is still present. All in all, we will end up with the conclusion that, without any a priori assumption about the relationship of the envelope and phase function of an oscillatory signal, there is not any preferred neither best representation of the IF of such oscillatory signal.

preprint2015arXiv

Extracting a shape function for a signal with intra-wave frequency modulation

In this paper, we consider signals with intra-wave frequency modulation. To handle this kind of signals effectively, we generalize our data-driven time-frequency analysis by using a shape function to describe the intra-wave frequency modulation. The idea of using a shape function in time-frequency analysis was first proposed by Wu. A shape function could be any periodic function. Based on this model, we propose to solve an optimization problem to extract the shape function. By exploring the fact that s is a periodic function, we can identify certain low rank structure of the signal. This structure enables us to extract the shape function from the signal. To test the robustness of our method, we apply our method on several synthetic and real signals. The results are very encouraging.

preprint2015arXiv

On the Uniqueness of Sparse Time-Frequency Representation of Multiscale Data

In this paper, we analyze the uniqueness of the sparse time frequency decomposition and investigate the efficiency of the nonlinear matching pursuit method. Under the assumption of scale separation, we show that the sparse time frequency decomposition is unique up to an error that is determined by the scale separation property of the signal. We further show that the unique decomposition can be obtained approximately by the sparse time frequency decomposition using nonlinear matching pursuit.

preprint2015arXiv

Optimal Local Multi-scale Basis Functions for Linear Elliptic Equations with Rough Coefficient

This paper addresses a multi-scale finite element method for second order linear elliptic equations with arbitrarily rough coefficient. We propose a local oversampling method to construct basis functions that have optimal local approximation property. Our methodology is based on the compactness of the solution operator restricted on local regions of the spatial domain, and does not depend on any scale-separation or periodicity assumption of the coefficient. We focus on a special type of basis functions that are harmonic on each element and have optimal approximation property. We first reduce our problem to approximating the trace of the solution space on each edge of the underlying mesh, and then achieve this goal through the singular value decomposition of an oversampling operator. Rigorous error estimates can be obtained through thresholding in constructing the basis functions. Numerical results for several problems with multiple spatial scales and high contrast inclusions are presented to demonstrate the compactness of the local solution space and the capacity of our method in identifying and exploiting this compact structure to achieve computational savings.

preprint2015arXiv

Self-similar Singularity of a 1D Model for the 3D Axisymmetric Euler Equations

We investigate the self-similar singularity of a 1D model for the 3D axisymmetric Euler equations, which is motivated by a particular singularity formation scenario observed in numerical computation. We prove the existence of a discrete family of self-similar profiles for this model and analyze their far-field properties. The self-similar profiles we find agree with direct simulation of the model and seem to have some stability.

preprint2015arXiv

Sparse Time-Frequency decomposition for multiple signals with same frequencies

In this paper, we consider multiple signals sharing same instantaneous frequencies. This kind of data is very common in scientific and engineering problems. To take advantage of this special structure, we modify our data-driven time-frequency analysis by updating the instantaneous frequencies simultaneously. Moreover, based on the simultaneously sparsity approximation and fast Fourier transform, some efficient algorithms is developed. Since the information of multiple signals is used, this method is very robust to the perturbation of noise. And it is applicable to the general nonperiodic signals even with missing samples or outliers. Several synthetic and real signals are used to test this method. The performances of this method are very promising.

preprint2014arXiv

A Heterogeneous Stochastic FEM Framework for Elliptic PDEs

We introduce a new concept of sparsity for the stochastic elliptic operator $-{\rm div}\left(a(x,ω)\nabla(\cdot)\right)$, which reflects the compactness of its inverse operator in the stochastic direction and allows for spatially heterogeneous stochastic structure. This new concept of sparsity motivates a heterogeneous stochastic finite element method ({\bf HSFEM}) framework for linear elliptic equations, which discretizes the equations using the heterogeneous coupling of spatial basis with local stochastic basis to exploit the local stochastic structure of the solution space. We also provide a sampling method to construct the local stochastic basis for this framework using the randomized range finding techniques. The resulting HSFEM involves two stages and suits the multi-query setting: in the offline stage, the local stochastic structure of the solution space is identified; in the online stage, the equation can be efficiently solved for multiple forcing functions. An online error estimation and correction procedure through Monte Carlo sampling is given. Numerical results for several problems with high dimensional stochastic input are presented to demonstrate the efficiency of the HSFEM in the online stage.

preprint2014arXiv

Sparse Time-Frequency decomposition by dictionary learning

In this paper, we propose a time-frequency analysis method to obtain instantaneous frequencies and the corresponding decomposition by solving an optimization problem. In this optimization problem, the basis to decompose the signal is not known. Instead, it is adapted to the signal and is determined as part of the optimization problem. In this sense, this optimization problem can be seen as a dictionary learning problem. This dictionary learning problem is solved by using the Augmented Lagrangian Multiplier method (ALM) iteratively. We further accelerate the convergence of the ALM method in each iteration by using the fast wavelet transform. We apply our method to decompose several signals, including signals with poor scale separation, signals with outliers and polluted by noise and a real signal. The results show that this method can give accurate recovery of both the instantaneous frequencies and the intrinsic mode functions.

preprint2013arXiv

Convergence of a data-driven time-frequency analysis method

In a recent paper, Hou and Shi introduced a new adaptive data analysis method to analyze nonlinear and non-stationary data. The main idea is to look for the sparsest representation of multiscale data within the largest possible dictionary consisting of intrinsic mode functions of the form $\{a(t) \cos(θ(t))\}$, where $a \in V(θ)$, $V(θ)$ consists of the functions smoother than $\cos(θ(t))$ and $θ'\ge 0$. This problem was formulated as a nonlinear $L^0$ optimization problem and an iterative nonlinear matching pursuit method was proposed to solve this nonlinear optimization problem. In this paper, we prove the convergence of this nonlinear matching pursuit method under some sparsity assumption on the signal. We consider both well-resolved and sparse sampled signals. In the case without noise, we prove that our method gives exact recovery of the original signal.

preprint2013arXiv

Generalized Multiscale Finite Element Methods (GMsFEM)

In this paper, we propose a general approach called Generalized Multiscale Finite Element Method (GMsFEM) for performing multiscale simulations for problems without scale separation over a complex input space. As in multiscale finite element methods (MsFEMs), the main idea of the proposed approach is to construct a small dimensional local solution space that can be used to generate an efficient and accurate approximation to the multiscale solution with a potentially high dimensional input parameter space. In the proposed approach, we present a general procedure to construct the offline space that is used for a systematic enrichment of the coarse solution space in the online stage. The enrichment in the online stage is performed based on a spectral decomposition of the offline space. In the online stage, for any input parameter, a multiscale space is constructed to solve the global problem on a coarse grid. The online space is constructed via a spectral decomposition of the offline space and by choosing the eigenvectors corresponding to the largest eigenvalues. The computational saving is due to the fact that the construction of the online multiscale space for any input parameter is fast and this space can be re-used for solving the forward problem with any forcing and boundary condition. Compared with the other approaches where global snapshots are used, the local approach that we present in this paper allows us to eliminate unnecessary degrees of freedom on a coarse-grid level. We present various examples in the paper and some numerical results to demonstrate the effectiveness of our method.

preprint2013arXiv

On Finite Time Singularity and Global Regularity of an Axisymmetric Model for the 3D Euler Equations

We investigate the large time behavior of an axisymmetric model for the 3D Euler equations. In \cite{HL09}, Hou and Lei proposed a 3D model for the axisymmetric incompressible Euler and Navier-Stokes equations with swirl. This model shares many properties of the 3D incompressible Euler and Navier-Stokes equations. The main difference between the 3D model of Hou and Lei and the reformulated 3D Euler and Navier-Stokes equations is that the convection term is neglected in the 3D model. In \cite{HSW09}, the authors proved that the 3D inviscid model can develop a finite time singularity starting from smooth initial data on a rectangular domain. A global well-posedness result was also proved for a class of smooth initial data under some smallness condition. The analysis in \cite{HSW09} does not apply to the case when the domain is axisymmetric and unbounded in the radial direction. In this paper, we prove that the 3D inviscid model with an appropriate Neumann-Robin boundary condition will develop a finite time singularity starting from smooth initial data in an axisymmetric domain. Moreover, we prove that the 3D inviscid model has globally smooth solutions for a class of large smooth initial data with some appropriate boundary condition.

preprint2013arXiv

On the Finite-Time Blowup of a 1D Model for the 3D Incompressible Euler Equations

We study a 1D model for the 3D incompressible Euler equations in axisymmetric geometries, which can be viewed as a local approximation to the Euler equations near the solid boundary of a cylindrical domain. We prove the local well-posedness of the model in spaces of zero-mean functions, and study the potential formation of a finite-time singularity under certain convexity conditions for the velocity field. It is hoped that the results obtained on the 1D model will be useful in the analysis of the full 3D problem, whose loss of regularity in finite time has been observed in a recent numerical study (Luo and Hou, 2013).

preprint2013arXiv

Potentially Singular Solutions of the 3D Incompressible Euler Equations

Whether the 3D 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 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 \times 10^{12})^{2}$ near the point of the singularity, we are able to advance the solution up to $τ_{2} = 0.003505$ and predict a singularity time of $t_{s} \approx 0.0035056$, while achieving a \emph{pointwise} relative error of $O(10^{-4})$ in the vorticity vector $ω$ and observing a $(3 \times 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.

preprint2012arXiv

Data-Driven Time-Frequency Analysis

In this paper, we introduce a new adaptive data analysis method to study trend and instantaneous frequency of nonlinear and non-stationary data. This method is inspired by the Empirical Mode Decomposition method (EMD) and the recently developed compressed (compressive) sensing theory. The main idea is to look for the sparsest representation of multiscale data within the largest possible dictionary consisting of intrinsic mode functions of the form $\{a(t) \cos(θ(t))\}$, where $a \in V(θ)$, $V(θ)$ consists of the functions smoother than $\cos(θ(t))$ and $θ'\ge 0$. This problem can be formulated as a nonlinear $L^0$ optimization problem. In order to solve this optimization problem, we propose a nonlinear matching pursuit method by generalizing the classical matching pursuit for the $L^0$ optimization problem. One important advantage of this nonlinear matching pursuit method is it can be implemented very efficiently and is very stable to noise. Further, we provide a convergence analysis of our nonlinear matching pursuit method under certain scale separation assumptions. Extensive numerical examples will be given to demonstrate the robustness of our method and comparison will be made with the EMD/EEMD method. We also apply our method to study data without scale separation, data with intra-wave frequency modulation, and data with incomplete or under-sampled data.

preprint2012arXiv

Finite Time Blow-up of a 3D Model for Incompressible Euler Equations

We investigate the role of convection on its large time behavior of 3D incompressible Euler equations. In \cite{HL09a}, we constructed a new 3D model by neglecting the convection term from the reformulated axisymmetric Navier-Stokes equations. This model preserves almost all the properties of the full Navier-Stokes equations, including an energy identity for smooth solutions. The numerical evidence presented in \cite{HL09a} seems to support that the 3D model may develop a finite time singularity. In this paper, we prove rigorously that the 3D inviscid model develops a finite time singularity for a family of smooth initial data whose energy is finite and conserved in time.

preprint2012arXiv

On Singularity Formation of a 3D Model for Incompressible Navier-Stokes Equations

We investigate the singularity formation of a 3D model that was recently proposed by Hou and Lei in [16] for axisymmetric 3D incompressible Navier-Stokes equations with swirl. The main difference between the 3D model of Hou and Lei and the reformulated 3D Navier-Stokes equations is that the convection term is neglected in the 3D model. This model shares many properties of the 3D incompressible Navier-Stokes equations. One of the main results of this paper is that we prove rigorously the finite time singularity formation of the 3D inviscid model for a class of initial boundary value problems with smooth initial data of finite energy. We also prove the global regularity of the 3D inviscid model for a class of small smooth initial data.

preprint2011arXiv

On the Local Well-posedness of a 3D Model for Incompressible Navier-Stokes Equations with Partial Viscosity

In this short note, we study the local well-posedness of a 3D model for incompressible Navier-Stokes equations with partial viscosity. This model was originally proposed by Hou-Lei in \cite{HouLei09a}. In a recent paper, we prove that this 3D model with partial viscosity will develop a finite time singularity for a class of initial condition using a mixed Dirichlet Robin boundary condition. The local well-posedness analysis of this initial boundary value problem is more subtle than the corresponding well-posedness analysis using a standard boundary condition because the Robin boundary condition we consider is non-dissipative. We establish the local well-posedness of this initial boundary value problem by designing a Picard iteration in a Banach space and proving the convergence of the Picard iteration by studying the well-posedness property of the heat equation with the same Dirichlet Robin boundary condition.

preprint2010arXiv

Dynamic Growth Estimates of Maximum Vorticity for 3D Incompressible Euler Equations and the SQG Model

By performing estimates on the integral of the absolute value of vorticity along a local vortex line segment, we establish a relatively sharp dynamic growth estimate of maximum vorticity under some assumptions on the local geometric regularity of the vorticity vector. Our analysis applies to both the 3D incompressible Euler equations and the surface quasi-geostrophic model (SQG). As an application of our vorticity growth estimate, we apply our result to the 3D Euler equation with the two anti-parallel vortex tubes initial data considered by Hou-Li. Under some additional assumption on the vorticity field, which seems to be consistent with the computational results of Hou-Li, we show that the maximum vorticity can not grow faster than double exponential in time. Our analysis extends the earlier results by Cordoba-Fefferman and Deng-Hou-Yu.

preprint2010arXiv

On Singularity Formation of a Nonlinear Nonlocal System

We investigate the singularity formation of a nonlinear nonlocal system. This nonlocal system is a simplified one-dimensional system of the 3D model that was recently proposed by Hou and Lei in [13] for axisymmetric 3D incompressible Navier-Stokes equations with swirl. The main difference between the 3D model of Hou and Lei and the reformulated 3D Navier-Stokes equations is that the convection term is neglected in the 3D model. In the nonlocal system we consider in this paper, we replace the Riesz operator in the 3D model by the Hilbert transform. One of the main results of this paper is that we prove rigorously the finite time singularity formation of the nonlocal system for a large class of smooth initial data with finite energy. We also prove the global regularity for a class of smooth initial data. Numerical results will be presented to demonstrate the asymptotically self-similar blow-up of the solution. The blowup rate of the self-similar singularity of the nonlocal system is similar to that of the 3D model.